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Chapter 1 
Introduction 



1.1 Summary 

In this Thesis we have investigated some effects appearing in top quark and Higgs bo- 
son decays with ffavor changing neutral currents (FCNC) in the framework of generic 
Two Higgs Doblet Models (2HDM) and the Minimal Supersymmetric Standard Model 
(MSSM). 

The Standard Model (SM) of Particle Physics interactions has had great success in 
describing the strong, weak and electromagnetic interactions and its validity has been 
tested up to the quantum level in past and present accelerators, such as the LEP at 
CERN or the Tevatron at Fermilab. The last great success of the SM was the discovery 
in 1994 of its last matter building block, namely the top quark, with a mass of = 
178.0 ± 2.7 ± 3.3 GeV. However, the mechanism by which all the SM particles get their 
masses is still unconfirmed, since no Higgs scalar has been found yet. The fermions couple 
to the Higgs bosons with a coupling proportional to their mass, so one expects that the 
large interactions between top quark and Higgs boson particles would give rise to large 
quantum effects. But the FCNC processes are very suppressed in the SM. So in some cases 
(specially the ones we are going to deal with in this Ph.D. Thesis) the sole observation of 
these FCNC processes would be instant evidence of new physics and could greatly help 
to unravel the type of underlying Higgs model. 

We have focused our work on two models: the generic Two Higgs Doblet Models 
(2HDM) and the Minimal Supersymmetric Standard Model (MSSM). The 2HDM is like 
the SM with an extra Higgs doblet, so it has more Higgs bosons and in particular some 
charged. There are two types, I and II, and they differ in the way they give mass to 
the fermions. The MSSM is an extension of the SM that incorporates Supersymmetry 
(SUSY). Supersymmetry is an additional transformation that can be added in the action 
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of Quantum Field Theory, leaving this action unchanged. The main phenomenological 
consequence is that for any SM particle (p) there should exist a partner for it, which we 
call sparticle (p) with different spin but with the same gauge quantum numbers. This 
extension of the SM provides elegant solutions to some theoretical problems of the SM, 
such as the hierarchy problem. 

We have applied these two extensions to the SM to see whether they can produce new 
FCNC effects. We have computed the following FCNC decays: (1) the branching ratios 
(B) of the top quark to Higgs bosons and charm quark in the 2HDM; (2) B and number of 
events at the LHC of the Higgs bosons to top and charm quarks in the 2HDM; (3) B of the 
Higgs bosons to bottom and strange quarks in the MSSM; (4) cross section and number 
of events at the LHC of the Higgs bosons to FCNC final states involving the heavy quarks 
like the top and bottom quark in the MSSM. We also studied the experimental signatures 
that would allow discover of the nature of these Higgs bosons in the LHC. In this study 
we have applied the severe restrictions from observed low-energy FCNC processes like 
b s'-f. 

Decays of the top quark induced by FCNC are known to be extremely rare events 
within the SM. This is so not only for the decay modes into gauge bosons, but most 
notably in the case of the Higgs channels, e.g. t — > H^'^^c, with a branching fraction of 
10~^^ at most. We have found that in the 2HDM the decays of the top quark to Higgs 
bosons, t H^, A^)c, can be the most favored FCNC modes - comparable or even 

more efficient than the gluon channel t — > gc. In both cases the optimal results are 
obtained for Type II models. However, only the Higgs channels can have rates reaching 
the detectable level (10~^), with a maximum of order 10~^. Compared with previous 
results obtained in the Higgs sector of the MSSM, the maximum branching ratios are 
similar but have different signatures. While in the 2HDM II there is only one Higgs boson 
that can reach the visible level or H^, but not both), in the MSSM all the channels 
can be competitive in some region. 

Similarly, Higgs boson decays mediated by FCNC are very much suppressed in the 
Standard Model, at the level of 10"^^ for Higgs boson masses of a few hundred GeV. 
We have computed the FCNC decays of Higgs bosons into a top quark, h ^ tc {h = 
h^,H^,A^), in a general 2HDM. The isolated top quark signature, unbalanced by any 
other heavy particle, is very clean (without background noise), and should help to identify 
the potential FCNC events much better than any other final state. We have computed 
the maximum branching ratios and the number of FCNC Higgs boson decay events at the 
LHC. The most favorable mode for production and subsequent FCNC decay is the lightest 
CP-even state (/i°) in the Type II 2HDM, followed by the other CP-even state (if°) if it is 
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not very heavy, whereas the CP-odd (A^) mode can never be sufficiently enhanced. Our 
calculation shows that the branching ratios of the CP-even states may reach 10~^, and 
that several hundreds of events could be collected in the highest luminosity runs of the 
LHC. We also point out some strategies in which to use these FCNC decays as a handle 
to discriminate between 2HDM and supersymmetric Higgs bosons. 

Furthermore, we analyzed the maximum branching ratios for the FCNC decays of the 
neutral Higgs bosons of the MSSM into bottom quarks, h ^ bs {h = h^,H^,A^), giving 
the maxima in the B{h ^ 6 s) ~ 10"'^ — 10""^ range. But this maximum could reach up 
to ~ 10~^ depending on whether or not it is allowed a fine-tunning in the B{b — > s'-f) 
restriction, which for naturalness reasons we do not allow. We consider that the bulk of 
the MSSM contribution to B{h —^bs) should originate from the strong supersymmetric 
sector, electroweak calculations are in progress. These calculations show that the FCNC 
modes h ^ bs can be competitive with other Higgs boson signatures and could play a 
helpful complementary role in identifying the supersymmetric Higgs bosons, particularly 
the lightest CP-even state in the critical LHC mass region m^jO ~ 90 — 130 GeV. 

Finally, we have also analyzed the production and subsequent FCNC decay of the 
neutral MSSM Higgs bosons to tc and bs in the LHC collider, pp h ^ tc, bs and h — > 
tc^bs {h = h'^,H^,A'^). Only the strongly-interacting FCNC sector has been computed 
because it expected to be the most important. We determined the maximum production 
rates for each of these modes and identified the relevant regions of the MSSM parameter 
space. The latter are different from those obtained by maximizing only the branching 
ratio, due to non-trivial correlations between the parameters that maximize/minimize 
each isolated factor. The production rates for the bs channel can be huge for a FCNC 
process (0.1 — 1 pb), but its detection can be problematic. The production rates for the tc 
channel are more modest (10~^ — 10~^pb), but its detection should be easier due to the 
clear-cut top quark signature. A few thousand tc events could be collected in the highest 
luminosity phase of the LHC, with no counterpart in the SM. 

Our general conclusion is that the physics of the processes with flavor changing neutral 
currents can be very important in seeing the physics beyond the Standard Model and to 
disentangle the nature of the most adequate model. Experiments at the LHC can be 
crucial to unravel signs of FCNC physics beyond the SM. 

1.2 Motivation 

The accepted model for the interactions between elementary particles is the Standard 
Model (SM) [1-5]. This model is composed of fundamental particles with spin 1/2 
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(fermions) , spin 1 (vector bosons) and one fundamental particle with spin 0, that is 
the Higgs boson, the only undiscovered particle of the SM. But there are good reasons 
to think that it is not the final model for the high energy physics. There are different 
models which try to explain the physics beyond the SM. The simplest model that extends 
the SM is the 2HDM [6]. The extra pieces of this model are some neutral and charged 
Higgs bosons. These neutral Higgs boson are those in which we are interested in order to 
compare with the SM Higgs boson. In this way we can distinguish which model describes 
better the future experimental results. Another not so simple model is the MSSM [7-10]^, 
and which has some characteristics of the 2HDM and, as we have said, has more or less 
twice the number of particles than the SM, as can be seen in chapter |21 

The top quark is the latest-discovered elementary particle of the SM. It almost com- 
pletes all of the building blocks of the SM. But for theoretical and aesthetic reasons we 
need to introduce a mechanism to give masses to the particles. One mechanism that 
allows this is the Higgs mechanism. It is based on introducing a new (undiscovered yet) 
particle of spin 0. So recently the Higgs boson has become the most wanted particle. The 
top quark physics is very important in the investigation of the high energy physics. It 
interacts with Higgs bosons with the highest possible strength (because it is proportional 
to the mass of the quark) and therefore this property may help to discover the Higgs bo- 
son. Moreover, it provides a big phase space, so it can decay in particles that the others 
cannot. 

The flavor changing neutral currents are a kind of processes, especially important to 
test the SM. They are characterised by one quark changing its flavor in the interaction 
(effective or not) with neutral particles (currents). Experimentally these processes are 
very suppressed, especially for the physics of the top quark, as we have seen. In the SM 
the branching ratios in the top quark case are so tiny that we cannot think of measuring 
them experimentally. But we could perhaps find models that, even being depressed, can 
give values nearer to our experimental possibilities. 

In the near and middle future, with the upgrades of the Tevatron (Run II, TeV33), the 
advent of the LHC, and the construction of an e+e~ Linear Collider (LC, nowadays called 
International Linear Collider ILC) [13-15], new results on top quark physics [16], and 
possibly on Higgs physics, will be obtained that may be extremely helpful complementing 
the precious information already collected at LEP I and II from Z and W physics. Both 
types of machines, the hadron colliders and the LC will work at high luminosities and 
produce large amounts of top quarks. In the LHC, for example, the production of top 
quark pairs will be a{tt) = 800 pb - roughly two orders of magnitude larger than in 

^For a review see [11,12] 
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the Tevatron Run II. In the so-called low-luminosity phase (10"^^ cm~^s~^) of the LHC, 
one expects about three tt-pairs per second (ten million tt-pairs per year!) [17]. And 
this number will be augmented by one order of magnitude in the high-luminosity phase 
(10^'^ cm~^s~^). As for a future LC running at e.g. = 500 GeV , one has a smaller 
cross-section a{tt) = 650 fb but a higher luminosity factor ranging from 5 x 10'^^ cm~'^s~^ 
to 5 X 10^^ cm^^s"^ and of course a much cleaner environment [18]. With datasets from 
LHC and LC increasing to several 100 fb'^/jeai in the high-luminosity phase, one should 
be able to pile up an enormous wealth of statistics on top quark decays. Therefore, not 
surprisingly, these machines should be very useful to analyze rare decays of the top quark 
and of the Higgs boson(s), viz. decays whose branching fractions are so small (< 10~^) 
that they could not be seen unless the number of collected decays is very large. 

The reason for the interest in these decays is at least twofold. First, the typical 
branching ratios for the rare top quark decays predicted within the Standard Model 
(SM) are so small that the observation of a single event of this kind should be "instant 
evidence", so to speak, of new physics; and second, due to its large mass (m^ = 178.0 ± 
2.7 ± 3.3 GeV [19]), the top quark could play a momentous role in the search for Higgs 
physics beyond the SM. While this has been shown to be the case for the top quark decay 
modes into charged Higgs bosons, both in the Minimal Supersymmetric Standard Model 
(MSSM) and in a general two-Higgs-doublet model (2HDM) [20,21]^, we expect that a 
similar situation would apply for top quark FCNC decays into non-SM neutral Higgs 
bosons and for Higgs boson FCNC decays. 

1.3 Today's situation 

The search for physics beyond the Standard Model (SM) is a very relevant, if not the most 
important, endeavor within the big experimental program scheduled in the forthcoming 
Large Hadron Collider (LHC) experiment at CERN [13,14]. There are several favorite 
searching lines on which to concentrate, but undoubtedly the most relevant one (due to 
its central role in most extensions of the SM) is the physics of the Higgs boson(s) with all 
its potential physical manifestations. 

As we mention above, experimentally, processes involving Flavor Changing Neutral 
Current (FCNC) have been shown to have rather low rates [19]. Letting aside the meson- 
meson oscillations, such as — and — B^, the decay processes mediated by FCNC 
are also of high interest and are strongly suppressed too. For instance, we have the 

^For a review of the main features of loop-induced supersymmetric effects on top quark production 
an decay, see e.g. Rcf. [22]. 
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radiative B-meson decays, with a typical branching ratio B{b — S7) ~ IQ—^. But we 
also have the FCNC decays with the participation of the top quark as a physical field, 
which are by far the most suppressed decay modes [23,24]. Indeed, the top quark decays 
into gauge bosons {t ^ cV; V = Z,g) are well known to be extremely rare events in 
the SM. The branching ratios are, according to Ref. [24]: ~ 5 x lO^^^ for the photon, 
slightly above 1 x IQ-^'^ for the Z-boson, and ~ 4 x IQ-^^ for the gluon channel, or even 
smaller according to other estimates [25]. Similarly, the top quark decay into the SM 
Higgs boson, H , is a very unusual decay, typically B{t cH^^) ~ IQ-^^ [26]. 





SM 


2HDM 


MSSM 


B{t^ C7) 


~ 5 X 10"i3 


< 1 X 10-^ 


< 1 X 10-6 


B{t^ cZ) 


> 1 X 10-13 


< 1 X 10-6 


< 1 X 10-^ 


B{t ^ eg) 


~ 4 X 10-11 


< 1 X 10-5 


< 1 X 10-5 


B{t^ cH) 


~ 10-13 - 10-1^ 


<io-4 


<io-4 


B{H ^ tc) 


~ 10-i3(mH < 2mw) 


<io-4 


<io-4 




< 10-i5(m^, > 2mw) 






B{H ^ bs) 


< 10-^(mH < 2Mw) 




<io-4 




< 10-iO(mH > mt) 







Table 1.1: Rare FCNC branching ratios of the top quark and the Higgs boson decays. 



The reason for this rareness is simple: for FCNC top quark decays in the SM, the loop 
amplitudes are controlled by down-type quarks, mainly by the bottom quark. Therefore, 
the scale of the loop amplitudes is set by ml and the partial widths are of order 
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and similarly for the FCNC Higgs boson decays 



r(i/-^ - tc) ^ (^)' run (Af-)^ ~ (M)' a.. (1.2) 



^tc) ~ ( aw Gprnl^lQ^^^ {ioiniH <2mw), (1.3) 



5(i75M^^-)^ (MV^^G^ I!^< 10-15 {ioimH>2mw). (1-4) 



B{H'^^ -.hs)^(^^}j ttvy (^^) < 10-^ (formH<2Mv^), (1.5) 

5(i7^A.^^-)^ I^M^a^G^ < 10-10 (formH>m,), (1-6) 

where a is a^m for ^ = 7, Z and ioiV = g,GF is Fermi's constant, a^y = g'^/4:7r and 
being the SU{2)i weak gauge coupling. Notice the presence of Af*'^ ~ nih/Mw, which is 
the SM Yukawa coupling of the bottom quark in units of g. The factor F ~ (1 — my/m^)^ 
results, upon neglecting m^, from phase space and polarization sums. Notice that the 
dimensionless fourth power mass ratio, in parenthesis in eq. (II. ip . and the fourth power 
of Af*^ stems from the GIM mechanism and is responsible for the ultra-large suppression 
beyond naive expectations based on pure dimensional analysis, power counting and CKM 
matrix elements. 

The GIM mechanism [27] is related to the unitarity of the mixing matrices between 
quarks. The Minimal Standard Model (SM) embeds the GIM mechanism naturally, due to 
the presence of only one Higgs doublet giving mass simultaneously to the down-type and 
the up-type quarks, and as a result no tree-level FCNCs interactions appear. FCNCs are 
radiatively induced, and are therefore automatically small. However, when considering 
physics beyond the SM, new horizons of possibilities open up which may radically change 
the pessimistic prospects for FCNC decays involving a Higgs boson and the top quark. 
Because of the loop-induced FCNCs effects, the SM and non-SM loops enter the FCNC 
observables at the same order of perturbation theory, and new physics competes on the 
same footing with SM physics to generate a non-vanishing value for these rare processes. 
It may well be that the non-SM effects are dominant and become manifest. Conversely, 
it may happen that they become highly constrained. 

The addition of further Higgs doublets to the SM in the most general way introduces 
potentially large tree-level FCNC interactions, which would predict significant FCNC rates 
in contradiction with observation. However, by introducing an ad-hoc discrete symmetry 
these interactions are forbidden. This gives rise to two classes of Two-Higgs-Doublet 
Models (2HDM) which avoid FCNCs at the tree-level, known conventionally as type I 
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and type II 2HDMs [6]. In Type I 2HDM (also denoted 2HDM I) one Higgs doublet, $i, 
does not couple to fermions at all and the other Higgs doublet, $2, couples to fermions 
in the same manner as in the SM. For more details see chapter |21 In contrast, in Type II 
2IIDM (also denoted 2HDM II) one Higgs doublet, $1, couples to down quarks (but not to 
up quarks) while $2 does the other way around. Such a coupling pattern is automatically 
realized in the framework of supersymmetry (SUSY), in particular in the MSSM, but it 
can also be arranged in non-supersymmetric extensions if we impose a discrete symmetry, 
e.g. $1 — $1 and $2 +^2 (or vice versa) plus a suitable transformation for the 
right-handed quark fields, this symmetry is only violated by soft terms of dimension two. 

Supersymmetry (SUSY) is certainly related to Higgs boson physics, and at the same 
time it may convey plenty of additional phenomenology. Ever since its inception, SUSY 
has been one of the most cherished candidates for physics beyond the SM, and as such 
it will be scrutinized in great detail at the LHC It is no exaggeration to affirm that the 
LHC will either prove or disprove the existence of SUSY, at least in its most beloved low- 
energy realization, namely the one which is needed to solve the longstanding naturalness 
problem in the Higgs sector of the SM [8] . On the other hand, SUSY provides an appealing 
extension of the SM, which unifies the fermionic and bosonic degrees of freedom of the 
fundamental particles and provides a natural solution to the hierarchy problem. The 
search for SUSY particles has been one of the main programs of the past experiments 
in high energy physics (LEP, SLD, Tevatron) and continues to play a central role in the 
present accelerator experiments (Tevatron II) and in the planning of future experimental 
facilities like the LHC and the LC. The Minimal Supersymmetric Standard Model (MSSM) 
is the simplest extension of the SM which includes SUSY, and for this reason its testing 
will be one of the most prominent aims of these powerful experiments. If SUSY is realized 
around the TeV scale, the LHC experiments shall be able to directly produce the SUSY 
particles for masses smaller than a few TeV [28,29]. On the other hand, the presence of 
SUSY may also be tested indirectly through the quantum effects of the supersymmetric 
particles. For one thing, it has been known since long ago that SUSY particles may 
produce large virtual effects on Higgs boson observables'^. 

In Ref. [42] it was shown that the vector boson modes can be highly enhanced within 
the context of the MSSM. This fact was also dealt with in great detail in Ref. [43] where 
a dedicated study was presented of the FCNC top quark decays into the various Higgs 
bosons of the MSSM (see also [44]) showing that these can be the most favored FCNC 
top quark decays - above the expectations on the gluon mode t ^ eg. For the 2HDM it 
was proven that while the maximum rates for t ^ eg were one order of magnitude more 

^See e.g. [20,30-41] and references therein. For a review see e.g. [11]. 
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favorable in the MSSM than in the 2HDM, the corresponding rates for t ^ ch^ were 
comparable both for the MSSM and the general 2HDM, namely up to the 10~^ level and 
should therefore be visible both at the LHC and the LC [25]. 

As in the 2HDM, in MSSM one has to impose some restrictions to avoid tree-level 
FCNCs among the extra predicted particles, which would induce one-loop FCNC interac- 
tions among the SM particles. But, in fact, MSSM requires their existence, because of the 
SU{2)l gauge symmetry. Of course, low energy measurements constrain the FCNC cou- 
plings (the most stringent being the B{b —>■ 57)). A potentially relevant FCNC interaction 
is the gluino with the squarks, not very much constrained experimentally. 

Concerning the FCNC interactions of Higgs bosons with third generation quarks, it was 
demonstrated long ago [43] that the leading term corresponds to a single particle insertion 
approximation. This produces a flavor change in the internal squark loop propagator, 
since in this case the chirality change can already take place at the squark-squark-Higgs 
boson interaction vertex. Adding this to the fact that the Higgs bosons (in contrast to 
gauge bosons) have a privileged coupling to third generation quarks, one might expect 
that the FCNC interactions of the type quark- quark- Higgs bosons in the MSSM become 
highly strengthened with respect to the SM prediction. This was already proven in the 
rare decay channels r(t ^ ch) [43] {h being any of the neutral Higgs bosons of the 
MSSM h = h^,H^,A^) where the maximum rate of the SUSY-QCD induced branching 
ratio was found to be B{t ch) ~ 10~^, eight orders of magnitude above the SM 
expectations B{t cH^^^) ~ 10"^'^. Similar enhancement factors have been found in 
the (top,bottom)-quark-Higgs boson interactions in other extensions of the SM, both in 
the MSSM (see chapters EE and Refs. [43,45-50]) and in the general two-Higgs-doublet 
model (2HDM) (see chapters |3]- El and Ref. [51]), and also in other extensions of the SM- 
see [25] for a review. 

The power of FCNC observables can be gauged e.g. by the implications of the bottom- 
quark rare decay h S7: the experimentally measured allowed range B{h S7) = 
(3.3 ± 0.4) X 10"^ [52-58] may impose tight constraints on extensions of the SM. For 
example, it implies a lower bound on the charged Higgs boson mass 'mH± > 350 GeV in 
general type II 2HDMs [59-62]. 

After this panoramic view of the FCNC processes in the SM and beyond, the work 
presented in this Thesis is as follows: in chapter |21 (resp. ^ we give the basic notations 
of the 2HDM (resp. MSSM); in chapter |3] we compute the FCNC top quark decay in the 
2HDM B(t — s> ch); in chapter^lwe compute the Higgs boson production and FCNC decay 
in the 2HDM a{pp ^ h ^ tc) at the LHC and B{h tc); in chapter IHl we compute the 
FCNC Higgs boson decay in the MSSM B{h — > bs); and in chapter [7| the Higgs boson 
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production and decay in the MSSM a{pp —^h-^ tc,bs) and B{h — > tc,bs). This PhD. 
work is based on the following articles [63-69]. 

In chapter |3] we show that within the simplest extension of the SM, namely the 
general two-Higgs-doublet model, the FCNC top quark decays into Higgs bosons, t — ^ 
[h^, H^, A^)c, can be the most favored FCNC modes - comparable or even more efficient 
than the gluon channel t — > gc. In both cases the optimal results are obtained for Type II 
models. However, only the Higgs channels can have rates reaching the detectable level 
10~^, with a maximum of order 10~^ which is compatible with the charged Higgs bounds 
from radiative B-meson decays. We compare with the previous results obtained in the 
Higgs sector of the MSSM. 

In chapter El we consider the FCNC decays of Higgs bosons into a top quark in a 
general two-Higgs-doublet model (2HDM). The isolated top quark signature, unbalanced 
by any other heavy particle, should help to identify the potential FCNC events much more 
than any other final state. We compute the maximum branching ratios and the number of 
FCNC Higgs boson decay events at the LHC collider at CERN. The most favorable mode 
for production and subsequent FCNC decay is the lightest CP-even state in the Type 
II 2HDM, followed by the other CP-even state, if it is not very heavy, whereas the CP- 
odd mode can never be sufficiently enhanced. Our calculation shows that the branching 
ratios of the CP-even states may reach 10~^, and that several hundred events could be 
collected in the highest luminosity runs of the LHC. We also point out some strategies to 
use these FCNC decays as a handle to discriminate between 2HDM and supersymmetric 
Higgs bosons. 

In chapterl^lwe analyze the maximum branching ratios for the FCNC decays of the neu- 
tral Higgs bosons of the MSSM into bottom and charm quarks, h ^ bs [h = h^, H^, A^). 
We consistently correlate these decays with the radiative B-meson decays {b — 57). A 
full-fledged combined numerical analysis is performed of these high-energy and low-energy 
FCNC decay modes in the MSSM parameter space. Our calculation shows that the avail- 
able data on B{b 57) severely restricts the allowed values of B{h bs). While the 
latter could reach a few percent level in flne-tuned scenarios, the requirement of natural- 
ness reduces these FCNC rates into the modest range B{h bs) ^ 10"^ — 10"'^. We 
expect that the bulk of the MSSM contribution to B{h bs) should originate from the 
strong supersymmetric sector. Our results are encouraging because they show that the 
FCNC modes h bs can be competitive with other Higgs boson signatures and could play 
a helpful complementary role to identify the supersymmetric Higgs bosons, particularly 
the lightest CP-even state in the critical LHC mass region m/jo ^ 90 — 130 GeV. 

In chapter [7| we analyze the production and subsequent decay of the neutral MSSM 



1.3 Today's situation 



11 



Higgs bosons {h = A^) mediated by FCNC in the LHC collider. We have 

computed the /i-pro duct ion cross-section times the FCNC branching ratio, a{pp ^ h ^ 
qq') = a{pp —^h)x B{h qq'), in the LHC focusing on the strongly-interacting FCNC 
sector. Here qq' is an electrically neutral pair of quarks of different flavors, the dominant 
modes being those containing a heavy quark: tc or bs. We determine the maximum 
production rates for each of these modes and identify the relevant regions of the MSSM 
parameter space, after taking into account the severe restrictions imposed by low energy 
FCNC processes. The analyses of a{pp h ^ qq') singles out regions of the MSSM 
parameter space different from those obtained by maximizing only the branching ratio, due 
to non-trivial correlations between the parameters that maximize/minimize each isolated 
factor. The production rates for the hs channel can be huge for a FCNC process (0.1 — 
Ipb), but its detection can be problematic. The production rates for the tc channel are 
more modest (10~^ — 10"^ pb), but its detection should be easier due to the clear-cut top 
quark signature. A few thousand tc events could be collected in the highest luminosity 
phase of the LHC, with no counterpart in the SM. 



Chapter 2 

Two Higgs Doublet Models (2HDM) 



2.1 Introduction 



The 2HDM are models that extend minimally the Higgs sector of the SM. They introduce 
one more doublet of complex scalar fields with hypercharge Y = +1. The most general 
Lagrangian with the SM gauge symmetry SU{2)l x U{1)y that contains these Higgs 
bosons can be divided in three terms: a kinetic term Ckin, the Yukawa couplings term 
(Higgs-fermions interactions) Cy and the potential for the two Higgs doublets V(0i,02): 



^Higgs — ^kin + — V(01, 02) , 



i=l,2 



a ^ Y 
D^ = df,- ig—W^ - ig'—B^ 



(2.1) 
(2.2) 

(2.3) 

where is the covariant derivative of SU{2)l ® U{1)y, and cij are the Pauli matrices^. 

The Higgs potential that spontaneously breaks the symmetry SU{2)l ® U{1)y to 
U{1)em is [6,70]: 

V(0i, 02) =Ai(0l0i - viy + A2(0|02 - vl)^+ 

+ A3 [(0l01 - Vl) + (0^02 - I'D 

+ A4 f(0i0l)(0|02) - (0l02)(0|0l) + (2.4) 



+ As 
+ Ae 



Re(0|02) - f if 2 cos ^ 
lm(0|02) - V1V2 sin^ 



+ 
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where all the Aj are real parameters, because the Lagrangian must be hermitic. This is the 
most general Lagrangian compatible with the gauge symmetry and the discrete symmetry 
01 — ^ —01 [6,70], this symmetry is only violated by soft terms of dimension two. We 
impose this last symmetry to forbid the FCNC at tree level. Moreover, this potential 
must be bounded from below, so the Aj must be non-negative. But in fact, the allowed 
range for the parameters Aj corresponding to this minimum is a range in the parameter 
space such that the square Higgs boson masses are positive and that \^(0, 0) > 0. 
In this context the minimum of the potential is: 



^2) 




(2.5) 
(2.6) 



which breaks the gauge symmetry giving U{1)em- 

We need two physical parameters in order to know their value, which are usually taken 
to be: 

1 „,2 „, _ 

(2.7) 



^ 2/2 

2^ 



g — , taxi(3 = — 
2 vi 



71 



' 2 



If we impose A5 = Ag (like in sypersymmetry) we can write the last two terms of 
eq. (El) as: 

0l02-fif2e*^r (2.8) 

The phase ^ can disappear with a redefinition of the fields without affecting the others 
terms of the potential (this phase will appear in other terms of the total Lagrangian). 
Then, the Higgs potential is CP conserving. 
At last, the final Higgs potential is: 



V(01, 02) =Ai(0l0i - Vlf + A2(0^02 - Vlf- 



+ A3 
+ A4 



402 



(0l0l)(0^02)-(0l02)(0^01 



+ 



(2.9) 



+ 



+ A5|0l02 - ViV2\^ 



There are different forms for the Yukawa terms of the Lagrangian to satisfy the 
Glashow and Weinberg [71] theorem. The Glashow and Weinberg theorem says that 
for a general SU{2) x U{1) gauge theory where we demand that the neutral-current inter- 
actions conserve all quark flavor naturally the necessary and sufficient conditions are: All 
quarks of fixed charge and helicity must (1) transform according to the same irreducible 
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representations of weak SU(2), (2) correspond to the same eigenvalue of weak T3, and 
(3) receive their contributions in the quark mass matrix form a single source (either from 
the vacuum expectations value of a single neutral Higgs boson or from a unique gauge- 
invariant bare mass term). In practice this implies that all fermions of a given electric 
charge couple to no more than one Higgs doublet. 

From the Lagrangian ()2.4j) with the potential ()2.9|) and the Yukawa terms we can 
obtain the full 2HDM spectrum, as well as the interactions, which contain the usual SM 
gauge interactions, the fermion-Higgs interactions, and the pure 2HDM interactions. A 
detailed treatment of this Lagrangian, and the process of derivation of the forthcoming 
results can be found in [72]. 



2.1.1 2HDM I 

In this model one of the Higgs doublets (^2) couple to all the fermions. The couplings 
with the quarks is of the form: 

^? = - E [^^- {itU^} qfn + (^f ^2) + h.c.] + leptons , (2.10) 



2J = 1 



where 



^ = ia2<p\ (2.11) 
- , (2.12) 

- (:) ' - (:) ' - ' 

and similarly for the leptonic doublets that contain the neutrinos and the leptons. 
This model is very related with the minimal model (SM), being the only difference a 
smaller VEV V2 < vsm {v ~ 174 GeV) and bigger Yukawa couplings. 



2.1.2 2HDM II 

Now one doublet {(pi) couples to the right-handed (RH) down fermions {qdR, Ida) and is 
responsible of the down masses; the other doublet (^2) couples to the RH up fermions 
{QuR^ Iur) and is responsible of their masses. Taking any flavor base, i.e. one in which /^*^ 
are isospin doublets the Lagrangian is: 

4''^ = - E (gIVi) qfn + U!, {<lf~^2) + h.c.] + leptons . (2.14) 
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The mass matrix will be proportional to the VEV of the Higgs as: 

M^"?'') = ws^^^"'^) , (2.15) 

M'f'^^ = wiD^"'^) , (2.16) 

This is basically the Higgs sector required in the MSSM. 
2.1.3 2HDM III 

This is the most general 2HDM without FCNC at tree level, being the other two important 
particular cases. The Yukawa interactions in this case, using any flavor base, is: 



4''^ = - E \dI,{<&^^& + (2.17) 

+ UlAlL^iWA + UlAlL^2)q^A + h.c.] + (2.18) 
+ [IHl terms] , (2.19) 

where the 3x3 matrices Di, D2, Ui, U2 are such that diagonalize simultaneously with 
the quark mass matrix. 



2.2 2HDM spectrum 

2.2.1 Higgs sector 

We will use this structure for the doublets: 

Re 0° + i Im ( 



1,2. (2.20) 



These fields are not physical fields, they do not have a well defined mass, as there are 
bilinear terms within the scalar fields with different fields. The next thing is to diagonalize 
the mass matrix. It can be seen that the mass matrix is a diagonal matrix in boxes for the 
fields a) (j)i,(j)2, b) lie (j)^, lie (j) 2 and c) Im^'j', Im</)2 (the real and imaginary part can be 
treated separately by CP invariance). So we have to separately diagonalize the different 
boxes. If we define the rotation angle as: 

, , / cos a; sincj\ , , 

R{u) = . , (2.21) 

\— smu; cos uo I 

the rotations (transformations) of the fields are: 
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AO 



= R{f3) 



a) 



Re0? 

Re 02 ^ f 2 y 

.0^ 



with their masses: 



1 



Mil + M22 ± \ {Mu - M22Y + 4Mi22 



where Mij are defined from the CP-even mass matrix 



M 



(4A3 + X^)viV2 Avl{\2 + A3) + vl\^ 



The mixing angles (3 and a are: 

tan l3 
sin 2q; 

cos 2a 



V2_ 



2M 



12 



V(Mii-M22)2+4Mf2 

Mil - M22 



V(Mii-M22)2+4M22' 

Now we can redefine the parameters of the theory as: 

4 masses : m/,o , m^/o , m^o , m/^± 
2 mixing angles : a, (3 



At tree level we get: 



Gf 9' 



V2 8M2,' 



(2.22) 
(2.23) 
(2.24) 



(2.25) 
(2.26) 

(2.27) 



(2.28) 



(2.29) 



(2.30) 



(2.31) 
(2.32) 



(2.33) 



where we find the value of v 



v{ + vi^ 



2^2-3/4^-1/2 ^174 6-61^. 



(2.34) 
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To obtain the magnitudes as functions of the physical magnitudes one has to invert 
the mass equations ()2.25|) and ()2.30|) : 

_ cos^am^o + sin^am^o vl 
_ sin^ a m|o + cos^ a vl 

AV2 AV2 

. m|o - mlo A5 . . 

A3 = cos asm a — , 2.37 

^ 4viV2 4 ' ^ ^ 



A4 = ^ , (2.38) 



A5 = ^ . (2.39) 

V 



2.3 Interactions in the mass-eigenstate basis 

We need to convert the interaction Lagrangian to a Lagrangian in the mass-eigenstate 
basis, which is the one used in the computation of the physical quantities. We quote only 
the interactions that we will need in our studies. 

• H^-Higgs: this interaction is obtained from the kinetic term of the Lagrangian: 




+ h. c. 

(2.40) 



• quarks-Higgs: they follow after replacing in ()2.10|) and ()2.14j) the mass-eigenstates 
Higgs fields ^TT^ : 



'-Htb 

nil 
'-Htb 

^hqq 

nil 
'-'hqq 



gVt 



tb 



H-b 



V2Mw 
-9 



nit cot P Pr + nih 



— cot P 
tan P 



2M, 



sin P 



cos a 



sm a 



w 



+ 



cos P 

1 grrih \ — cot P 

2 Mw \ tan/5 | " " " ' 2Mwtanp 
2 Mw sin p 



Pl 



t + h.c. (2.41) 



sma 



cos a 



t cos a + sin a\ t . 



(2.42) 



2.3 Interactions in the mass-eigenstate basis 



19 





+ (m^o ~ '^/fo) cos2/5sin(/3 — a)] 




MvKsm2/3 [("^^± m^o + 2"^'o)sin2/5sm(/3 a) + 
+ f?T2,?n — Tn?An^ COS 2/3 cosf/5 — au 




-llgS [(2m^o + m^o) sin 2a- 
— m^o (3 sin 2a — sin 2/3)] 




— „ , , ^ . \vn?jjn sin 2/3 cosf/5 — a^H- 

2Mvi/sm2/3 L -ff" ■^^^^ ''A^ ^^■^VA^ ^) i 

+2(m^o — "o^^o) cos 2/3 sin(/5 — a)] 




2M/sm2/3 Ko sin2/?sin(/5 a) + 
+2(m^o ~ "^^o) cos 2/3 cos(/5 — a)] 


H^-H^ - A^ 





H^-G^ - 


in \ 

{ ^a)(m^ .sin(/^-a) 


H^-G^ - 


5 r, cosfp — a) 


H^-G^ - A^ 


/" 2 2 N 


G^-G^ - 


, . . m^o cos(/3 - a) 
( 2M,. 


G^-G^ - 


^, m^o sin(/3 - a) 
^ 2Mw 


G^-G^ - A^ 






Table 2.1: Feynman rules for the trilinear couplings involving the Higgs self-interactions 
and the Higgs and Goldstone boson vertices in the Feynman gauge, with all momenta 
pointing inward. These rules are common to both Type I and Type II 2HDM under the 
conditions explained in the text. We have singled out some null entries associated to CP 
violation. 



where we have used the third quark family, the Vth is the corresponding element of 
the CKM matrix and Pl,r = T 75) are the chiral projectors. 



Trilinear Higgs couplings: they are summarized in 12. ![ and are valid for Type I and 
Type II models. Had we not imposed the restriction A5 = Ag, then the trilinear 
rules would be explicitly dependent on the A5 parameter. 
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2.4 Constraints 

There are multiple constraints that must be imposed, obviously one of such constraints 
is that they must reproduce the behaviour of the SM up to energy scales probed so far. 
Analysing the perturbativity of the theory one finds that the allowed range for tan (3 

is: 

0.1 < tan/3 < 60. (2.43) 

The custodial symmetry [73,74] {SU{2)) is a good symmetry at tree level, so the 
quadratic violations of this symmetry must be experimentally fixed. So, the one-loop 
corrections of the parameter p from the 2HDM sector can not be bigger than one per mil 
oftheSM [75]: 

|^^2HDM| ^ 00^_ (2.44) 

To be precise, the latter is the extra effect that 5p can accommodate at one standard 
deviation (1 a) from the 2HDM fields beyond the SM contribution [43]. This is a stringent 
restriction that affects the possible mass splittings among the Higgs fields of the 2HDM, 
and its implementation in our codes does severely prevent the possibility from playing 
with the Higgs boson masses to artificially enhance the FCNC contributions. 

Moreover, the charged Higgs bosons have an important indirect restriction from the 
radiative decays of the B meson, specially the ratio B{B — Xs7) - or B{h — * S7) at the 
quark level [52-58]: 

5(fe^ S7) = (3.3±0.4) X 10"^ (2.45) 

The Higgs contribution to B{h sj) (that have been computed at the NLO in QCD [62]) 
is positive: bigger experimental ratio means that the charged Higgs mass can be smaller. 
From the different analysis of the literature [59-62] we get: 

mH± > 350 GeV (2.46) 

for virtually any tan/5 > 1. This bound does not apply to Type I models because at 
large tan P the charged Higgs couplings are severely suppressed, whereas at low tan /3 we 
recover the previous unrestricted situation of Type II models. 

We can derive lower bounds for the neutral Higgs masses in these models [76,77]. 
For example, using the Bjorken process e"^e^ ^ Z + and the production of pairs of 
Higgs boson e~^e~ h^{H^) + we can get the following restrictions in almost all the 
parameter space [78,79]: 

f> 100 GeV Vtan/3 
rriho + rajifi < . (2.47) 

150 GeV tan/5 >1 
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In each of these cases there is a small region in the parameter space in the ranges of the 
CP-even Higgs masses and CP-odd masses around [78,79]: 

m;,,o^o = 20 - 30 GeV . (2.48) 

Although, as can be seen in the electroweak precision fits in Ref. [80], in the high tan/? 
range a light Higgs boson is statistically correlated with a light -ff^, so this situation is 
not favoured by the b s'-f restriction. Moreover, since our interest in Type II models is 
mainly focused in the large tan (3 regime, the corner in the light CP-even mass range is a 
bit contrived. At the end of the day one finds that, even in the worst situation, the strict 
experimental limits still allow generic 2HDM neutral scalar bosons as light as 70 GeV or 
so. As we said, most of these limits apply to Type II 2HDM's, but we will conservatively 
apply them to Type I models as well. 

Finally, the unitarity bound can be approximately formulated by imposing that the 
absolute value of the trilinear coupling of the 2HDM Higgs can not be bigger than the 
Trilinear coupling of the SM Higgs: 



Chapter 3 



The Minimal Supersymmetric 
Standard Model (MSSM) 

3.1 Introduction 

It goes beyond the scope of this Thesis to study the formal theory of Supersymmetry 
[81,82], however we would like, at least, to give a feeling on what is it. Supersymmetry 
(SUSY) can be introduced in many manners. 

Let us consider the symmetries of the scattering matrix S, that is, those transforma- 
tions that can be reduced to an interchange of asymptotic states. Before the discovery of 
Supersymmetry the only symmetries know were: the following: (1) the ones correspond- 
ing to the Poincare group; (2) the so called internal global symmetries, both of them ruled 
by a Lie algebra; and (3) discrete symmetries such as parity (P), charge conjugations (C) 
and the time reversal (T). A 1967 theorem due to Coleman and Mandula establishes rig- 
orously that, under very general conditions, these are the only symmetries allowed for S if 
we do not want to induce trivial scattering (fixed angles and speeds) in 2 — * 2 processes. 

This theorem may be eluded relaxing some of its hypotheses. The Supersymmetry 
appears precisely when assuming that the generators of the new symmetry we want to 
add have a spinorial character instead of a scalar one, therefore transforming under (|, 0) 
and (0, ^) representations of the Lorentz group. Fermionic spinorial generators necessarily 
have an anti-commutative algebra, generically known as a graded Lie algebra. The algebra 
is not closed with just the SUSY generators, thus it can not be understood as an internal 
symmetry, but it rather forms an extension of the space-time symmetries of the Poincare 
group. 

Following this line of thought, one could relax some other hypotheses of the Coleman- 
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Mandula theorem in order to introduce new theories. SUSY is the only known extension 
for the S matrix symmetries. Accepting as the only valid extensions of the Coleman- 
Mandula theorem conditions the presence of a graded Lie algebra, one can show (Haag, 
Lopuszahki and Sohnius theorem) that spinorial generators different from those of SUSY 
are forbidden. 

Next, we have to define the "superspace" , the supersymmetric space of the "super- 
fields" . We add to the space-time coordinates x other sets of spinorial coordinates (6*, 9) 
(as many sets as the dimension of the space-time) that are Grassmann variables, i.e. 
they anti-commute. In the case of adding just one set it is said that we have a = 1 
Supersymmetry and a = 1 superspace: 

space — time — ^ A^ = 1 Superspace 

^ (a;^,^",^°) (3.1) 

where a = 1,2. The supersymmetric transformations have the parameters (A,a,^,^), 
where A is the Lorentz matrix, a is the translation 4-vector, and the Weyl spinors ^, ^. 
The generators of SUSY transformations consist of are the ones of the Poincare group 
and the new spinorial generators Qa and Qa, satisfying the graded Lie algebra. The 
infinitesimal purely SUSY transformation of a superfield <I>(x, 6*, 9) is: 

$ ^ $ + (3.2) 

5s = -i{^''Qa + LQ") (3.3) 

The functions defined in the Superspace are polynomial functions of the (0, 9) variables 
(since 9^^ = 9\ = Q). Thus we can decompose the functions (superfields) of this Superspace 
in components of 9^, 9^, 9^, 9a9p, . . . each of these components will be a function of the 
space-time coordinates. Analogously to the space-time, we can define in the Superspace 
scalar superfields, vector superfields, . . . For example in a 4-dimensional space time with 
A^ = 1 supersymmetry a scalar superfield has 10 components. 

We can define fields with specific properties with respect to the 9 variables [81,82]. A 
scalar chiral field in a 4D A^ = 1 Superspace has 4 components: 

^L = A + V29ip + 99F = {A,ip,F) (3.4) 
^r = A* + V29iIj + 99F* = {A\iP,F*) (3.5) 

where A is a scalar field, ip and ip are Weyl spinors (left-handed and right handed 
Dirac fermions) and F is an auxiliary scalar field. This auxiliary field is not a dynamical 
field since its equations of motion do not involve time derivatives. To this end we are left 
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with a superfield, whose components represent an ordinary scalar field and an ordinary 
chiral spinor. So if nature is described by the dynamics of this field we would find a chiral 
fermion and a scalar with identical quantum numbers. That is Supersymmetry relates 
particles which differ by spin 1/2. When a SUSY transformation [Q) acts on a superfield 
it transform spin s particles into spin s ± 1/2 particles. 

Thus, for a = 1 SUSY, we find that to any chiral fermion there should be a scalar 
particle with exactly the same properties. This fact is on the basis of the absence of 
quadratic divergences in boson mass renormalization, since for any loop diagram involving 
a scalar particle there should be a fermionic loop diagram, which will cancel quadratic 
divergences between each other, though logarithmic divergences remain. 

Supersymmetric interactions can be introduced by means of generalized gauge trans- 
formations, and by means of a generalized potential function, the Superpotential, which 
give rise to masses, Yukawa-type interactions, and a scalar potential. 

As no scalar particles have been found at the electroweak scale we may infer that, if 
SUSY exists, it is broken. We can allow SUSY to be broken maintaining the property that 
no quadratic divergences are allowed: this is the so called Soft-SUSY-Breaking mechanism 
[83]. We can achieve this by only introducing a small set of SUSY-Breaking terms in the 
Lagrangian, to wit: masses for the components of lowest spin of a supermultiplet and triple 
scalar interactions. However, other terms like explicit fermion masses for the matter fields 
would violate the Soft-SUSY-Breaking condition. 

The MSSM is the minimal Supersymmetric extension of the Standard Model. It is 
introduced by means of a = 1 SUSY, with the minimum number of new particles. Thus 
for each fermion / of the SM there are two scalars related to its chiral components called 
"sfermions" for each gauge vector V there is also a chiral fermion: "gaugino" ({)), 

and for each Higgs scalar H another chiral fermion: "higgsino" {h). In the MSSM it turns 
out that, in order to be able of giving masses to up-type and down- type fermions, we 
must introduce two Higgs doublets with opposite hypercharge, and so the MSSM Higgs 
sector is of the so called Type II 2HDM (see chapter |21 section 13.4. II and Ref. [6]). 

To build the MSSM Lagrangian we must build a Lagrangian invariant under the 
gauge group SU(3)c; x SU(2)i x U(l)y, it must also include the superfields with the 
particle content of the Table 13.11 and in addition it must contain the terms that breake 
supersymmetry softly. But this Lagrangian violates the baryonic and leptonic number, so 
we have to introduce an additional symmetry. In the case of the MSSM this symmetry 
is the so-called i?-symmetry. In its discrete form it relates the spin (S*), the baryonic 
number {B) and the leptonic number (L) in the so-called -R-parity: 
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so that is 1 for the SM fields and —1 for its supersymmetric partners. In the way the MSSM 
is implemented i?-parity is conserved, this means that i?-odd particles (the superpartners 
of SM particles) can only be created in pairs, also that in the final product decay of an 
-R-odd particle at least one SUSY particle exists, and that the Lightest Supersymmetric 
Particle (LSP) is stable. 

3.2 Field content 

The field content of the MSSM consist of the fields of the SM plus all their supersymmetric 
partners, and an additional Higgs doublet. The Table ITTl shows all the correspondences 
and all the fields. All these fields suffer some mixing, so the physical (mass eigenstates) 
fields look much different from these ones, as shown in Table The gauge fields mix up 
to give the well known gauge bosons of the SM, W"^, A^, the gauginos and higgsinos 
mix up to give the chargino and neutralino fields, and finally the Left- and Right-chiral 
sfermions mix among themselves in sfermions of indefinite chirality. Letting aside the 
intergenerational mixing between fermions and sfermions that give rise to the well known 
Cabibbo-Kobayashi-Maskawa (CKM and superCKM) matrix. 



Super field 


SM particle 


Sparticle 


SU(3)c 


SU(2)i 


U(l)y 


Matter 
















L 
R 


leptons < 
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sleptons 
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U = u*n 
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1 
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Hi 

H2 


Higgs< 

[H2={H+,H^,) 
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'h, = {hIh^) 

H2 = {H^,Hl) 


1 
1 


2 
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-1 
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G 
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V 
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1 


3 





V' 




b 5° 




bino {B^ 


1 


1 






Table 3.1: Particle contents of the MSSM superfields 
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Name 


Mass eigenstates 


Gauge eigenstates 


Higgs bosons 


/^O ^0 ^0 JJ± 


Hi Hi H7 Ht 

12 12 


squarks 


ii h hi 62 


h in hi Ir 


sleptons 


n f2 J>r 


Tl Tr Ur 


neutralinos 


Ni N2 N3 Ni 


B° lyo Hi Hi 


charginos 


Ct Ct 


Hi Ht 



Table 3.2: Mass eigenstates of the MSSM particles. For notational simplicity only the 
third sfermion generation is presented. 



3.3 Lagrangian 



The MSSM interactions come from three different kinds of sources: 



Superpotential: 



fH\UR + hdHlQW + huHiQ'U - ^iH{Hi 



(3.7) 



The superpotential contributes to the interaction Lagrangian ()3.14|) with two differ- 
ent kind of interactions. The first one is the Yukawa interaction, which is obtained 
from (j3.7|) just replacing two of the superfields by its fermionic field content, whereas 
the third superfield is replaced by its scalar field content: 



Vy 



+eij 



fH\UR + hdH\QW + KHIQ'U - fiHlHi 
fH\UR + hdHlQW + KHiQ'U 
fH{UR + hdHlQW + KHiQ'U 



(3. 



+ h.c. 



The second kind of interactions are obtained by means of taking the derivative of 
the superpotential: 



w 



E 



dW {if) 



difi 



(3.9) 



ipi being the scalar components of superfields. 



Interactions related to the gauge symmetry, which contain: 



— the usual gauge interactions 
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the gaugino interactions: 

Vg^^ = ^V2ga^,y {T'^),, + h.c. (3.10) 

where {fjip) are the spin and spin 1/2 components of a chiral superfield 
respectively, is a generator of the gauge symmetry, is the gaugino field 
and g"' its coupling constant. 

and the D-terms, related to the gauge structure of the theory, but that do not 
contain neither gauge bosons nor gauginos: 



with 

D'^ = g'^^UT%^j , (3.12) 
(fi being the scalar components of the superfields. 

Soft-SUSY-Breaking interaction terms: 



miAiH\UR + m^A^HiQW - m^A^HiQU 



+ h.c. . 

(3.13) 

The trilinear Soft-SUSY-Breaking couplings Af can play an important role, specially 
for the third generation interactions and masses, and they are in the source of the 
large value of the bottom quark mass renormalization effects. 



The full MSSM Lagrangian is then: 

-^MSSM = -^Kinetic + -^Gaugc — ^Gff — — Vy — 



dW{ip) ^ 



difi 



-m\ L*L - m| R*R - m| Q*Q - ml U*U - m\ D*D , (3.14) 

where we have also included the Soft-SUSY-breaking masses. 

From the Lagrangian (j3.14j) we can obtain the full MSSM spectrum, as well as the 
interactions, which contain the usual SM gauge interactions, the fermion-Higgs interac- 
tions that correspond to a Type II Two-Higgs-Doublet Model [6], and the pure SUSY 
interactions. A very detailed treatment of this Lagrangian, and the process of derivation 
of the forthcoming results can be found in [84] . 
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3.4 MSSM spectrum 
3.4.1 Higgs boson sector 

As seen in 12.11 when a Higgs doublet is added to the SM there exist two possibihties for 
incorporating it, avoiding Flavour Changing Neutral Currents (FCNC) at tree level [6]. 
The first possibility is not to allow a coupling between the second doublet and the fermion 
fields, this is the so called Type I 2HDM. The second possibility is to allow both Higgs 
doublets to couple with fermions, the first doublet only coupling to the Right-handed 
down-type fermions, and the second one to Right-handed up-type fermions, this is the so 
called Type II 2HDM. 

The Higgs sector of the MSSM is that of a Type II 2HDM [6], with some SUSY 
restrictions. After expanding ()3.14|) the Higgs potential reads 

V = m\\Hi\^ + ml\H2\^ -ml^ [eijH\Hi + \i.c.) 

+ \{g' + g''){\H,\'-m''Y^\g'\HlH,\\ (3.15) 

This is equivalent to the 2HDM potential (j2.9|) with the following restrictions: 

Ai = A2 (3.16) 

\, = \{g' + g")-\i (3.17) 

A4 = 2Ai - (3.18) 

A5 = A6 = 2Ai-^(/ + (?'2). (3.19) 

The neutral Higgs bosons fields acquire a vacuum expectation value (VEV), 

(/^i)o=(';) Wo=(") ■ (3.20) 

We need two physical parameters in order to know their value, which are usually taken 
to be Mw and tan (3: 

Ml^ = \g\vl + vl)^g''^ (3.21) 

Ml = \{9^ + 9'V = cos' Ow (3.22) 
tanp = —, 0<p<^ (3.23) 

tangly = - (3.24) 
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These VEV's make the Higgs fields to mix up. There are five physical Higgs fields: a 
couple of charged Higgs bosons {H^); a "pseudoscalar" Higgs {CP = —1) ^4°; and two 
scalar Higgs bosons {CP = 1) H° (the heaviest) and h° (the lightest). There are also the 
Goldstone bosons and G^. The relation between the physical Higgs fields and that 
fields of is 





(3.25) 



cos /3 —sm/3 
sin f3 cos f3 

+ ± I I ■;„ I (3.26) 

v2 \sin/3 cos (3 

I /-(cos/5 -sin/5)\ /GO' 
^V^\ sin/3 cos/3 J \A° 

where a is given in (j3.28|) [6]. 

All the masses of the Higgs sector of the MSSM can be obtained with only two pa- 
rameters, the first one is tan/3 ()3.24p . and the second one is a mass; usually this second 
parameter is taken to be either the charged Higgs mass or the pseudoscalar Higgs 
mass rriAO. We will take the last option. From ()3.15|1 one can obtain the tree- level mass 
relations between the different Higgs particles, 

m^± = m\o + , 
"^HO,ho = I (rnlo + M| ± (m^o + M|)' - 4m\o M| cos2 2/3^ , (3.27) 

and the mixing angle between the two scalar Higgs is obtained by means of: 

cos 2a = - cos 2/3 ^^^^ f , sin 2a = - sin 2/3 ^ \ . (3.28) 

The immediate consequence of such a constrained Higgs sector, is the existence of absolute 
bounds (at tree level) for the Higgs masses: 

< niho < Mz < mno, Mw<mH±- (3.29) 

It must be taken into account, though, that the radiative corrections, mainly due to the 
top-stop supermultiplet, and also the bottom-sbottom one, are susceptible of relaxing the 
limits ()3.29|) in a significant manner. 

A good approximation for effective mixing angle aes including only the leading one- 
loop contributions of top, stop, bottom and sbottom follows from the diagonalization of 
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the one-loop Higgs mass matrix [85-88]: 



A/f2 

Higgs 



where 



sin 2(3 I cot {3 M| + tan {3 m\t) + at + uJb 



-Ml 



m\o + At + Aft 



-Ml 



m\o + At + Afe 



tan P Ml + cot /5 m^o + ^^t + K 



V2n^ sin^ /5 



log 



,mt^mt^^ ^ A(A-yUCot/5) ^^^ "^i^ 



2 2 
II 12 



At 



y27r2 sin^ p 



(m? — m? )2 
'/i(A-/icot/3) 

log 



2 2 



log 



t2 



t'Z 



2 2 

— m~ 



tl t2 t2 

2^At{At- II cot (3f 



tl t2' 
2 /' /I , , R\2 



il i2_ 

2 2 

TTlj — m- 

tl t2 



log- 



mr 



m 



t2 



0-t 



NcGprnf {At - jj, cot 13) 



y27r2 sin^ 



(m? — m? )2 

^ tl i2 ' 



2 2 



tl 



t2 



ni- 



t2 



with the following substitutions for the bottom/sbottom factors: 

t 



b 

cos P 

Ab — yutan/3 



{ub, Afo, ab) ^ (ujt, At, at) with <^ sin P 

At — fi cot P 

This approximate effective mixing angle a^ff is determined by 

{m\o + Ml) tan/5 + (At + Afe)(l + ta.ii^p)/2 



tana 



Ml + m^o tan2 p + (at + tOb - Ml,^^^^){l + tan^/?) 
where M^o^eff is the solution for the light Higgs mass: 



771% + Ml + LJt + at + LJb + ab 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



HO,hO,eff 



± 



(m^o + M|)2 + (cjt -at + ab- cUb)' 
4 



-(m^o-M|) 



-m^oM|cos^2/5+ 0^1^^i±^|^) eos2/5,^, 

(At + A.sin2/j)^^^^^^,^^(At + . 
2 4 



(At + Xb)' 



1/2 



(3.34) 



In the limit where m^o ^ M^, cos(/9 — a) = (9(M|/m^o), which means that the 
couplings to Standard Model particles approach values corresponding precisely to the 
couplings of the SM Higgs boson. There is a significant region of MSSM Higgs sector 
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0.1 
0.01 
-J 0.001 

I 

CO. 

S 0.0001 

le-05 
le-06 
le-07_. 

Figure 3.1: The value of cos^(/9 — a) is shown as a function of m^o for two choices of 
tan/9 = 5 and tan/5 = 50. When radiative-corrections are included, one can define an 
approximate loop-corrected angle a as a function of m^o, tan (3 and the MSSM parameters. 
In the figures above, we have incorporated radiative corrections, assuming that Msusy = 
rrtq = = rriu = ITeV. The decoupling effect, in which cos^(/3 — a) oc M|/m^o^ for 
m^o ^ mz, continues to hold even when radiative corrections are included. 



parameter space in which the decoupling limit applies, because cos(/3 — a) approaches zero 
quite rapidly once m^o is larger than about 200 GeV, as shown in Fig. 13.11 As a result, 
over a significant region of the MSSM parameter space, the search for the lightest CP-even 
Higgs boson of the MSSM is equivalent to the search for the SM Higgs boson. This result 
is more general; in many theories of non-minimal Higgs sectors, there is a significant 
portion of the parameter space that approximates the decoupling limit. Consequently, 
simulations of the SM Higgs signal are also relevant for exploring the more general Higgs 
sector. 

3.4.2 The SM sector 

In this section we give some expressions to obtain some MSSM parameters as a function 
of the SM parametrization. 

As stated above (sec. 13.4. Ij) the VEV's can be obtained by means of (|3.24j) . and the 
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Z mass can be obtained at tree-level by the relation: 

sm ew = l-^^. 

Fermion masses are obtained from the Yukawa potential ()3.8p letting the neutral Higgs 
fields acquire their VEV ()H.2()|1 . The up-type fermions get their masses from the H2 
whereas gives masses to down-type fermions, so 

/i„V2Miysin/3 , hdV2Mw cos (3 

niu = Kv2 = , rrid = hsJi = , 

9 9 

and the Yukawa coupling can be obtained as 

^ _hu _ rriu X _ _ rud . , 

9 V2MvKsin/3 9 V^Mw cos f3 

3.4.3 Sfermion sector (Flavor-diagonal case) 

The sfermion mass terms are obtained from the derivative of the superpotential the 
D-terms (jH.llj) and the Soft-SUSY-Breaking terms ()3.14|1 letting the neutral Higgs fields 
get their VEV ()3.20|) . and one obtain the following mass matrices: 



" ' rn.Mlj, +m2 + cos 2/? s^j, M|; ' 



1 \ ^ n^i 



being Q the corresponding fermion electric charge, the third component of weak 
isospin, Mq^ ^ the Soft-SUSY-Breaking squark masses [7-10] (by SU (2) ^-gauge invariance, 
we must have M^^ = M^^, whereas Mf^, M^^ are in general independent parameters), 
sg = sin 9w, and 

M^^ = A„-/icot/3 , 

Mlji = Ad- tan (3 . (3.37) 

We define the sfermion mixing matrix as [q'^ = {q'l = Ql, q'2 = ^-r} are the weak- 
eigenstate squarks, and qa = {qi,q2} are the mass-eigenstate squark fields) 

b 

i„\ I COS 9a — sin 6a \ , , 

= ' '] . (3.38) 

\sin^q cos6'q 

/?Wt_^2^(<7) = diag{m? , m? } {m^, > m,-J , (3.39) 
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tan 26 q 




(3.40) 




+ cos2/3(T| 




From eq. ()3.36p we can see that the sfermion mass is dominated by the Soft-SUSY- 
Breaking parameters (Mj ^ for / ^ top), and that the non-diagonal terms could 
be neglected, except in the case of the top squark (and bottom squark at large tan/5), 
however we will maintain those terms, the reason is that, although the A parameters do 
not play any role when computing the sfermion masses, they do play a role in the Higgs- 
sfermion-sfermion coupling, and thus it has an effect on the Higgs self-energies. Moreover 
these A parameters are constrained by the approximate (necessary) condition of absence 
of colour-breaking minima, 



where rriq is of the order of the average squark masses for q = t,b [89-92]. 

All the Soft-SUSY-Breaking parameters are free in the strict MSSM, however some 
simplifications must be done to be able of making a feasible numerical analysis. As the 
main subject of study are the third generation squarks we make a separation between 
them and the rest of sfermions. This separation is justified by the evolution of the squark 
masses from the (supposed) unification scale down to the electroweak scale [12] (see also 
section IB. 7. II for a more detailed discussion). 

So we will use the following approximations: 

• equality of the diagonal elements of eq. ()3.36|) 



for each charged slepton and each squark of the the first and second generation. 

• the first and second generation squarks share the same value of the A parame- 
ter ()3.37p and the mass parameter ()3.42|) 

• sleptons also share the same value for ()3.42|) and A parameters ()3.37|) . 

3.4.4 Sfermion sector (Non-flavor-diagonal case) 

Very important for our FCNC studies is when the squark mass matrix does not diagonalize 
with the same matrix as the one for the quarks. We introduce then intergenerational mass 
terms for the squarks, but in order to prevent the number of parameters from being too 




(3.41) 



(3.42) 



3.5 Interactions in the mass-eigenstate basis 



35 



large, we have allowed (symmetric) mixing mass terms only for the left-handed squarks. 
This simplification is often used in the MSSM, and is justified by RGE analysis [93]. 
The flavor mixing terms are introduced through the parameters 5ij defined as 

{Ml^)ij = m^j = SijUiiTTij , (3.43) 

where rrii is the mass of the left-handed i squark, and mfj is the mixing mass matrix 
element between the generations i and j. Thus we must diagonalize two 6x6 mass 
matrices in order to obtain the mass-eigenstates squark fields. Generalizing the notation 
in Sec. 13.4.31 we introduce the 6x6 mixing matrices as follows: 

Qa = (3.44) 

R('')^MlR = M\n = diag{m2^, . . . ,ml} , q = u, d , (3.45) 

where -A^^^j-, is the 6x6 square mass matrix for up-type (or down-type) squarks in 
the EW basis, with indices a = 1, 2, 3, . . . , 6 = ul, ur, ■ ■ ■ .t^ for up-type squarks, and 
similarly for down-type squarks. In this study we are only interested in the up-type quarks- 
squarks system, so we will drop out the (g) super-index in the forthcoming expressions. 
The rotation matrix R introduces gluino mediated tree-level FCNC between quarks and 
squarks. 

To analyze the contributions from these flavor and chiral mixed squarks we can use the 
so-called mass insertion approximation. This is based on the fact that the 5ij parameters 
are small, so instead of diagonalizing the 6x6 squared mass matrix we can treat them 
as an interaction (see Fig. 13. 2j) : 



^ (^^ (4 "?) fc) 

(3.46) 

3.5 Interactions in the mass-eigenstate basis 

We need to convert the MSSM interaction Lagrangian ()3.14|1 to a Lagrangian in the mass- 
eigenstate basis, which is the one used in the computation of the physical quantities. As 



TiL TjL 



' -mf \ f QiL 
p2 - J yqji 



and similarly for the left-right squark matrix. 
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Figure 3.2: Diagrams in the mass insertion approximation 



the expression for the full interaction Lagrangian in the MSSM is rather lengthy we quote 
only the interactions that we will need in our studies. 

• quark-quark-neutral Higgs: this is the usual Yukawa interaction from Type II 
2HDM, in the MSSM it follows after replacing in ()3.8|) the mass-eigenstate Higgs 
fields 



^Hqn = 7^ — - — t; \ (cos a — siu a dd — i sin j3 d'-f^d A'^] 

2Mw cosp '-^ ^ 

^ " [(sin a + cos a h^) uu — i cos f3 u'y^u A^~^ , (3.47) 



2Mw sin p 

where we have replaced the Yukawa couplings hi in favour of masses and tan (]. 
quark-gluon interactions: this is the usual QCD Lagrangian 

^QCD = |G^A?^.g,7^g, . (3.48) 

quark-squark-gluino: the supersymmetric version of the strong interaction is ob- 
tained from ()3.10|) : 

^gqq = [R*5a " R^a Pr] Q*a,i Kj tj 



- ^ ij! [Rl, Pl - i?L Pr] . (3.49) 

= -^g;, (A=),^. (4f Pz. - R^1*Pr) qj + h.c. , (3.50) 



V2 

where A*^ are the Gell-Mann matrices 
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squark-squark-neutral Higgs: the origin of this interaction is twofold, on one side the 
superpotential derivative ()3.9|) . and on the other the Soft-SUSY-Breaking trihnear 
interactions. It is convenient to express the resuhs in the — Qr basis, and use the 
transformation 13.441 or 13.381 

^hqg = [(T*^ - Q^Se^)q*LqiL + QiSB^TiRQiR\ 

i=u,d 



X 



cos(/3 + a) - /i° sin(/? + a)) 



2 

^^'^ -X^ldL + d*RdR) cos a - /i° sin a) 



Mw cos (3 

2 

— - — 7—s{uIul + u*rUr){H'^ COS a - /i°sina) (3.51) 
IVI]^ sm p 

^^^^idndL + dldR) 

X [(— /i sin a + Ad cos a)H^ + (— cos a — Ad sin 

X [(— /i cos a + Au sin a)H^ + sin a + Au cos a)/?."] 

^(/i + A,tan/3)(44 - 
igrriu 



+ Au cot (3) {u*jiUL-uluR) A 



3.6 Flavor changing neutral currents 

The most general MSSM includes tree-level FCNCs among the extra predicted particles, 
which induce one-loop FCNC interactions among the SM particles. Given the observed 
smallness of these interactions, tree-level SUSY FCNCs are usually avoided by including 
one of the two following assumptions: either the SUSY particle masses are very large, and 
their radiative effects are suppressed by the large SUSY mass scale; or the soft SUSY- 
breaking squark mass matrices are aligned with the SM quark mass matrix, so that both 
mass matrices are simultaneously diagonal. However, if one looks closely, one soon realizes 
that the MSSM does not only include the possibility of tree-level FCNCs, but it actually 
requires their existence [93]. Indeed, the requirement of SU{2)i gauge invariance means 
that the up-left-squark mass matrix can not be simultaneously diagonal to the down-left- 
squark mass matrix, and therefore these two matrices can not be simultaneously diagonal 
with the up-quark and the down-quark mass matrices, that is, unless both of them are 
proportional to the identity matrix. But even then we could not take such a possibility 
too seriously, for the radiative corrections would produce non-zero elements in the non- 
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Figure 3.3: Feynman graphs contributing to the squark mixing. 



qi = {b + s + d) 



Figure 3.4: Tree level FCNC strong interactions. 



diagonal part of the mass matrix (i.e. induced by H"^ and see Fig. 13. 3|) . All in all, 
we naturally expect tree-level FCNC interactions mediated by the SUSY partners of the 
SM particles. As an example in the MSSM that one can not set the FCNC Higgs bosons 
interactions to zero without being inconsistent notice that r(t — > cx°) is UV divergent in 
the absent of these couplings [94]. The potentially largest FCNC interactions are those 
originating from the strong supersymmetric (SUSY-QCD) sector of the model (viz. those 
interactions involving the squark-quark-gluino couplings, see Fig. 13. 4|) . and in chapters 
ini and [7| we mainly concentrate on those. These couplings induce FCNC loop effects on 
more conventional fermion-fermion interactions, like the gauge boson-quark vertices Vqq'. 

Of course, low energy meson physics puts tight constraints on the possible value of the 
FCNC couplings, especially for the first and second generation squarks which are sensitive 
to the data on K'^ — (see Fig. 13. 5|) and — [95,96]. The third generation system 
is, in principle, very loosely constrained since present data on — mixing still leaves 
much room for FCNCs and the most stringent constraints are given by the B[h ^ 57) 
measurement [52-58]. Therefore, the relevant FCNC gluino coupling ^23 [95,96] is not 
severely bound at present. The lack of tight FCNC constraints in the top-bottom quark 
doublet enables the aforementioned lower bound on the charged Higgs boson mass in the 
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MSSM to be easily avoided, to wit: by arranging that the SUSY-electroweak (SUSY- 
EW) contribution to B{h 57) from the top-squark/chargino loops screens partially 
the charged Higgs boson contribution. This situation can be naturally fulfilled if the 
higgsino mass parameter (/i) and the soft SUSY-breaking top-squark trilinear coupling 
[At) satisfy the relation /lA^ < [97-101]. The relevant Wilson operator in the effective 
theory involves a chirality fiip: 

^7 = TTT^ (^^ ^""^ ^r) ^m- (3-52) 

In Fig. 13. (il one can see the leading contributions to this operator. The most important 
contribution is expected to be SUSY-QCD one. Roughly speaking its amplitude goes like: 

j,svsy-^cn^, ^ ^ - ^an/j) ^ ^ 

This amplitude can be as large as minus twice the SM contribution, so the total MSSM 
amplitude can respect the experimental 6 — > S7 restriction. We would regard this choice as 
a fine-tuning of the parameters, hence unnatural, but in fact it is excluded experimentally 
[102]. 

On the other hand, if we assume that the squark mass matrices are diagonal at a certain 
energy scale (A) and write the Renormalization Group Equations [93] , then non-diagonal 
terms are generated in the left-left sector, because A„ and are not simultaneously 
diagonal. 

(M^^)2 = /io(4)l + /ii(4)AiAd 

{M^j^y = /io(«ij)l + yUi(S^)AtA„ 
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Figure 3.6: Leading contributions to the b ^ s'y Wilson operator O7. 
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^ \ ^ 1- „i- „ ^ 2 2 



diagonalize ^ = ± A (3.55) 



And then the FCNC terms are communicated from the up to the down sector by the 
CKM matrix, due to the SU{2) gauge invariant: 

The FCNC gluino interactions also induce large contributions to B{b s'-f), see 
Fig. 13. (il It should be noted however that the leading contributions to the Vqq' FCNC 
interactions from the third quark generation correspond to a double insertion term, in 
which the squarks propagating in the loop suffer a double mutation: a flavor conversion 
and a chirality transition. This fact has been demonstrated in the B{b ^ s'-f) observable 
itself [103], as well as in the FCNC rare decay width r(t eg) [43]. As a consequence, 
the loose limits on the third generation FCNC interactions derived under the assumption 
that the leading terms contributing to b ^ S'-f correspond to the single particle insertion 
approximation [95, 96] are not valid and more complex expressions must be taken into 
account [104]. 

The FCNC effects in the K^ — K^ system can be described with an effective Lagrangian 
with a mass matrix as: 

A mj^ 

so the signal of FCNC can determined by the mass difference of the eigenstates: 

ArriK = \mK^ — m/^jl (3.56) 

In contrast the experimental signal for b s'-f is directly its branching ratio. The exper- 
imental measurements are [19,53-55,57]: 

S(6^S7) = (3.3±0.4) X 10"^ (3.57) 
AmK = (3.483 ± 0.006) x 10"^^ MeV (ds ^ ds) (3.58) 

rriK = 497.648 ± 0.022 MeV 
Amo < 4.6 X 10"^^ MeV {cu ^ cu) (3.59) 

mo = 1864.6 ±0.5 MeV 
AniB = (3.304 ± 0.046) x 10"^° MeV {bd ^ Id) (3.60) 

rriB = 5279.4 ±0.5 MeV 

that translated to Sij restrictions they read: 

512 < .1 V7F^^/500GeV 

513 < .098^/m^/500GeV (3.61) 
^23 < 8.2m,-mi/(500GeV)2 
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3.7 MSSM parametrization 
3.7.1 MSSM parameters 

Since SUSY is a broken symmetry we have to deal with a plethora of Soft-SUSY-Breaking 
parameters, namely 

• tan/5, 

• masses for Left- and Right-chiral sfermions, 

• a mass for the Higgs sector, 

• gaugino masses, 

• triple scalar couplings for squarks and Higgs. 

This set of parameters is often simplified to allow a comprehensive study. Most of 
these simplifications are based on some universality assumption at the unification scale. 
In minimal supergravity (mSUGRA) all the parameters of the MSSM are computed from 
a restricted set of parameters at the Unification scale, to wit: tan/3; a common scalar 
mass mo; a common fermion mass for gauginos mi/2; a common trilinear coupling for all 
sfermions Aq\ and the higgsino mass parameter /i. Then one computes the running of 
each one of these parameters down to the EW scale, using the Renormalization Group 
Equations (RGE), and the full spectrum of the MSSM is found. 

We will not restrict ourselves to a such simplified model. As stated in the introduc- 
tion we treat the MSSM as an effective Lagrangian, to be embedded in a more general 
framework that we don't know about. This means that essentially all the parameters 
quoted above are free. However for the kind of studies we have performed there is an 
implicit asymmetry of the different particle generations. We are mostly interested in the 
phenomenology of the third generation, thus we will treat top and bottom supermultiplet 
as distinguished from the rest. This approach is well justified by the great difference of 
the Yukawa couplings of top and bottom with respect to the rest of fermions. We are 
mainly interested on effects on the Higgs sector, so the smallness of the Yukawa couplings 
of the first two generations will result on small effects in our final result. We include them, 
though, in the numerical analysis and the numerical dependence is tested. On the other 
hand, if we suppose that there is unification at some large scale, at which all sfermions 
have the same mass, and then evolve these masses to the EW scale, then the RGE have 
great differences [12]. Slepton RGE are dominated by EW gauge interactions, 1st and 2nd 
generation squarks RGE are dominated by QCD, and for the 3rd generation squarks there 
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is an interplay between QCD and Yukawa couplings. Also, as a general rule, the gauge 
contribution to the RGE equations of left- and right-handed squark masses are similar, 
so when Yukawa couplings are not important they should be similar at the EW scale. 

With the statement above in mind we can simplify the MSSM spectrum by taking 
an unified parametrization for 1st and 2nd generation squarks (same for sleptons). As 
we have seen in Sect. 13.4.31 we will use: a common mass^ for Ul and {rriu); an unified 
trilinear coupling for 1st and 2nd generation; a common mass for all z/^ and 1^ (m^); 
and a common trilinear coupling A^. 

For the third generation we will use different trilinear couplings At and Ab, as these 
can play an important role in the kind of processes we are studying (see chapter IHl and 
[7j). Stop masses can present a large gap (due to its Yukawa couplings), being the right- 
handed stop the lightest one. We will use a common mass for both chiral sbottoms, which 
we parametrize with the lightest sbottom mass {m-^), and the lightest stop quark mass 
(m^-), as the rest of mass inputs in this sector. This parametrization is useful in processes 
where squarks only appear as internal particles in the loops (such as the ones studied in 
chapters ini and IZj) , as one- loop corrections to these parameters would appear as two-loop 
effects in the process subject of study. Gluino mass {fng), on the other hand, is let free. 

For the Higgs sector two choices are available, we can use the pseudoscalar mass niAO, 
or the charged Higgs mass mH±. Both choices are on equal footing. As the neutral Higgs 
particles are the main element for most of our studies we shall use its mass as input 
parameter in most of our work. 

Standard model parameters are well known, we will use present determinations of EW 
observables [19] 

Mz = 91.1899 ± 0.0021 GeV 
Mw = 80.418 ±0.054 GeV 
Gf = (1.16639 ± 0.00001) X 10~^GeV~2 
a;^{Mz) = 128.896 ±0.090 

(3.62) 

QCD related observables are not so precise. On the other hand as the main results are 
not affected by specific value of these observables we will use the following ones 

mt = 175 GeV 
rUb = 5 GeV 
as{Mz) = 0.11 (3.63) 
^Notc that after diagonalization of the squark mass matrix the physical masses will differ slightly. 
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3.7.2 Constraints 

The MSSM reproduces the behaviour of the SM up to energy scales probed so far [105]. 
Obviously this is not for every point of the full parameter space! 

There exists direct limits on sparticle masses based on direct searches at the high 
energy colliders (LEP II, SLC, Tevatron). Although hadron colliders can achieve larger 
center of mass energies than e+e~ ones, their samples contain large backgrounds that 
make the analysis more difficult. This drawback can be avoided if the ratio signal-to- 
background is improved, in fact they can be used for precision measurements of "known" 
observables (see e.g. [17]). e+e~ colliders samples are more clean, and they allow to put 
absolute limits on particle masses in a model independent way. 

No significant evidence for a Higgs signal has been detected at LEP [106]. As a result, 
one can obtain bounds on the possible MSSM Higgs parameters. These limits are often 
displayed in the m^o-tan [3 plane, see Fig. 13. 7t although there is additional dependence 
on various MSSM parameters that effect the radiative corrections to the Higgs masses 
as discussed above. In representative scans of the MSSM parameters, the LEP Higgs 
Working Group [106] finds that m^Q > 91.0 GeV and m^o > 91.9 GeV at 95% CL. 
These limits actually correspond to the large tan/3 region in which Zh^ production is 
suppressed, as shown in Fig. 13.71 In this case, the quoted Higgs limits result 
of the non-observation of h^A^ and H^A^ production. As tan /3 is lowered, the limits on 
m/jO and m^o become more stringent. In this regime, the h^A^ production is suppressed 
while the Zh^ production rate approaches its SM value. Thus, in this case, the SM Higgs 
limit applies {niho > 114 GeV [19]) as shown in Fig. I3.7r a). The precise region of MSSM 
Higgs parameter space that is excluded depends on the values of the MSSM parameters 
that control the Higgs mass radiative corrections. For example, a conservative exclusion 
limit is obtained in the maximal mixing scenario, since in this case the predicted value 
of rriho as a function of m^o and tan f] is maximal (with respect to changes in the other 
MSSM parameters). The excluded regions of the MSSM Higgs parameter space based 
on the maximal mixing benchmark scenario of Ref. [107], are shown in Fig. 13.71 and 
correspond to the exclusion of the range 0.5 < tan/3 < 2.4 at the 95% CL. However, the 
tan/3 exclusion region can still be significantly reduced (even to the point of allowing all 
tan/3 values) by, e.g., taking Msusy = 2 TeV and rrit = 180 GeV, see Fig 13.81 (which 
still lies within the error bars of the experimentally measured value), and allowing for the 
theoretical uncertainty in the prediction of m/jO™^^ [108]. 

Similarly, no evidence for the charged Higgs boson has yet been found. The LEP 
Higgs Working Group quotes a limit of m//± > 78.6 GeV at 95% CL [109], which holds 
for a more general non-supersymmetric two-Higgs doublet model and assumes only that 
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Figure 3.7: LEP2 contours of the 95% CL exclusion limits for MSSM Higgs sector 
parameters as a function of tan/3 and (a) niho and (b) m^o (in GeV), taken from Ref. [106]. 
The contours shown have been obtained for MSSM Higgs parameters chosen according to 
the maximal mixing benchmark of Ref. [107]. 



CO. 



5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 
1.0 
0.5 



HHW, M3U3^ = 1000 GeV 






H = - 100 GeV 






m, = 174.3 GeV 






m, = 179.4 GeV 






m, = 184.5 GeV 










/ / 






' / - 






/ 

/ 






/ 

/ 

/ 

/ 

/ 




/ 

/ 





70 



80 



90 



100 



110 



120 130 



m, [GeV] 



Figure 3.8: Theoretically exclusion limits as in Fig. 13.71 for three different values of the 
top-quark mass, rrit = 174.3, 179.4, 184.5 GeV. Taken from Ref. [108] 
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the decays dominantly into r+z/,- and/or cs. Although the MSSM tree-level bound 
mH± > Mw can be relaxed somewhat by radiative corrections, the LEP bound quoted 
above provides no useful additional constraints on the MSSM Higgs sector. Actual fits 
to the MSSM parameter space project a preferred value for the charged Higgs mass of 
mH± ^ 120 GeV [110]. 

Hadron colliders bounds are not so restrictive as those from e~^e~ machines. Most 
bounds on squark and gluino masses are obtained by supposing squark mass unification 
in simple models, such as mSUGRA. At present approximately the limits, with certain 
conditions, on squarks (1st and 2nd generation) and gluino masses are [19] 



From the top quark events at the Tevatron a limit on the branching ratio B(t ^ h) 
can be extracted, and thus a limit on the tan/5 — mH± relation [20,32]. 

Finally indirect limits on sparticle masses are obtained from the EW precision data. 
We apply these limits through all our computations by computing the contribution of 
sparticles to these observables and requiring that they satisfy the bounds from EW mea- 
surements. We require new contributions to the p parameter to be smaller than present 
experimental error on it, namely 



We notice that as 5pncw is also the main contribution from sparticle contributions to 
Ar [111], new contributions to this parameter are also below experimental constrains. 
Also the corrections in the a- and G^-on-shell renormalization schemes will not differ 
significantly [112]. 

We recall that at the moment the MSSM may escape the indirect bound from h 
s'-f because the positive charged Higgs virtual contributions can be compensated for by 
negative stop and chargino loops, if they are not too heavy. Therefore, in the MSSM 
the charged Higgs can stay relatively light, m/^± > 120 GeV^, just to comply with the 
aforementioned LEP 200 bounds on m^o [113]. 

There exist also theoretical constrains to the parameters of the MSSM. As a matter 
of fact the MSSM has a definite prediction: there should exist a light neutral scalar Higgs 
boson h^. Tree-level analysis put this bound to the Z mass, however the existence of large 
radiative corrections to the Higgs bosons mass relations grow this limit up to ~ 130 GeV 
[87,88]. More recently the two- loop radiative corrections to Higgs mass relations in the 
MSSM have been performed [114-118], and the present upper limit on m^jO is 



Trig > 250 GeV , > 195 GeV . 



(3.64) 



|5pnew| < 0.003 (3fr). 



(3.65) 



ruho < 130 - 135 GeV . 



(3.66) 
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Of course if the MSSM is extended in some way this hmit can be evaded, though not to 
values larger of ~ 200 GeV [119, 120]. 

Another theoretical constraint is the necessary condition ()3.41|) on squark trilinear 
coupling (A) to avoid colour-breaking minima. This constraint is easily implemented 
when the A parameters are taken as inputs, but if we choose a different set of inputs 
(such as the mixing angle 6g) then it constrains the parameter space in a non-trivial way. 

Whatever the spectrum of the MSSM is, it should comply with the benefits that SUSY 
introduces into the SM which apply if the following condition is fulfilled: 

MsusY < Oil TeV) . (3.67) 

If supersymmetric particles have masses heavier than the TeV scale then problems with 
GUT's appear. This statement does not mean that SUSY would not exist, but that then 
the SM would not gain practical benefit from the inclusion of SUSY. 

A similar upper bound is obtained when making cosmological analyses, in these type 
of analyses one supposes the neutralino to be part of the cold dark matter of the universe, 
and requires its annihilation rate to be sufficiently small to account for the maximum 
of cold dark matter allowed for cosmological models, while at the same time sufficiently 
large so that its presence does not becomes overwhelming. Astronomical observations also 
restrict the parameters of SUSY models, usually in the lower range of the mass parameters 
(see e.g. [121-123]). 

For the various ROE analyses to hold the couplings of the MSSM should be pertur- 
bative all the way from the unification scale down to the EW scale. This implies, among 
other restrictions, that top and bottom quark Yukawa couplings should be below certain 
limits {h^fj/Air < 1). In terms of tan/5 this amounts it to be confined in the approximate 
interval 

l<tan/5< 70. (3.68) 

There are additional phenomenological restrictions that bring the lower bound on tan (3 
to roughly 2.4 for the so-called maximal m° scenario, and 10.5 for the so-called no mixing 
scenario [113], the upper bound being the same [11]. 

In the MSSM case we use the more restrictive limits ()3.68p . see Fig. 13.71 whereas in 
the 2HDM model the lower limit can be smaller as shown in ()2.43|) . Any deviation from 
this framework of restrictions will only be for demonstrational purpouses, and will be 
explicitly quoted in the text. 
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Loop Induced FCNC Decays of the 
Top Quark in a General 2HDM 

4.1 Introduction 

As we have said in the introduction, in the LHC the production of top quark pairs will 
be a{tt) = 860 pb - roughly two orders of magnitude larger than in the Tevatron Run 
II. In the so-called low-luminosity phase (10'^^ cm~^s~^) of the LHC one expects about 
three tt-pair per second (ten million tt-pairs per year!) [17]. And this number will be 
augmented by one order of magnitude in the high-luminosity phase (10^^ cm~^s~^). As for 
a future LC running at e.g. ^/s = 500 GeV, one has a smaller cross-section a{tt) = 650 fb 
but a higher luminosity factor ranging from 5 x 10"^^ cm"^s~^ to 5 x 10'^^ cm"^s~^ and 
of course a much cleaner environment [18]. With datasets from LHC and LC increasing 
to several 100/6~^/year in the high-luminosity phase, one should be able to pile up an 
enormous wealth of statistics on top quark decays. Therefore, not surprisingly, these 
machines should be very useful to analyze rare decays of the top quark, viz. decays whose 
branching fractions are so small (< 10"^) that they could not be seen unless the number 
of collected top decays is very large. 

The situation is dramatic with the top quark decay into the SM Higgs boson, B{t 
cHsm) ~ 10~^3 _ 10-15 (mt = 175 GeV; Mz < Mh < 2 Mw) [26]. This extremely tiny 
rate is far out of the range to be covered by any presently conceivable high luminosity 
machine. On the other hand, the highest FCNC top quark rate in the SM, namely that of 
the gluon channel t ^ eg, is still 6 orders of magnitude below the feasible experimental 
possibilities at the LHC. All in all the detection of FCNC decays of the top quark at 
visible levels (viz. B{t —>■ eX) > 10~^) by the high luminosity colliders round the corner 
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(especially LHC and LC) seems doomed to failure in the absence of new physics. Thus 
the possibility of large enhancements of some FCNC channels up to the visible threshold, 
particularly within the context of the general 2HDM and the MSSM, should be very 
welcome. Unfortunately, although the FCNC decay modes into electroweak gauge bosons 
Vew = W, Z may be enhanced a few orders of magnitude, it proves to be insufficient to raise 
the meager SM rates mentioned before up to detectable limits, and this is true both in the 
2HDM - where B{t V^^ c) < lO'^ [24] - and in the MSSM - where B{t V^u, c) < 10"^ 
except in highly unlikely regions of the MSSM parameter space [124]^. In this respect 
it is a lucky fact that these bad news need not to apply to the gluon channel, which 
could be barely visible {B{t —^gc)< 10~^) both in the MSSM [42,43] and in the general 
2HDM [24]. But, most significant of all, they may not apply to the non-SM Higgs boson 
channels t H^, A^) + c either. As we shall show in the sequel, these Higgs decay 

channels of the top quark could lie above the visible threshold for a parameter choice 
made in perfectly sound regions of parameter space! 

While a systematic discussion of these "gifted" Higgs channels was made in Ref. [43] 
for the MSSM case and in other models^, to the best of our knowledge there is no similar 
study in the general 2HDM. And we believe that this study is necessary, not only to assess 
what are the chances to see traces of new (renormalizable) physics in the new colliders 
round the corner but also to clear up the nature of the virtual effects; in particular to dis- 
entangle whether the origin of the hypothetically detected FCNC decays of the top quark 
is ultimately triggered by SUSY or by some alternative, more generic, renormalizable 
extension of the SM such as the 2HDM or generalizations thereof. Of course the alleged 
signs of new physics could be searched for directly through particle tagging, if the new 
particles were not too heavy. However, even if accessible, the corresponding signatures 
could be far from transparent. In contrast, the indirect approach based on the FCNC 
processes has the advantage that one deals all the time with the dynamics of the top 
quark. Thus by studying potentially new features beyond the well-known SM properties 
of this quark one can hopefully uncover the existence of the underlying new interactions. 



^Namely, regions in which there are wave- function renormaUzation thresholds due to (extremely for- 
tuitous!) sharp coincidences between the sum of the sparticle masses involved in the self-energy loops 
and the top quark mass. See e.g. Ref. [125] for similar situations already in the conventional t ^ W b 
decay within the MSSM. In our opinion these narrow regions should not be taken too seriously. 

^Preliminary SUSY analysis of the Higgs channels are given in [126], but they assume the MSSM Higgs 
mass relations at the tree-level. Therefore these are particular cases of the general MSSM approach given 
in [43]. Studies beyond the MSSM (e.g. including R-parity violation) and also in quite different contexts 
from the present one are available in the literature, see e.g. [127]. 
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4.2 Relevant fields and interactions in the 2HDM 

We will mainly focus our interest on the loop induced FCNC decays 

t^ch {h = h°,H\A^) , (4.1) 

in which any of the three possible neutral Higgs bosons from a general 2HDM can be 
in the final state. However, as a reference we shall compare throughout our analysis the 
Higgs channels with the more conventional gluon channel 

t^cg. (4.2) 

Although other quarks could participate in the final state of these processes, their contri- 
bution is obviously negligible - because it is further CKM-suppressed. The lowest order 
diagrams entering these decays are one-loop diagrams in which Higgs, quarks, gauge and 
Goldstone bosons - in the Feynman gauge - circulate around. While the diagrams for the 
decays (j4.ip are depicted in Fig. 14.11 the ones for the decay (j4.2j) are not shown [24] . 

We wish to concentrate here on Type I and Type II models of a sufficiently generic 
nature, to wit, those which are characterized by the following set of free parameters: 

{rriho, mjjo,mj^o,mH±, tan a, tan P) , (4.3) 

tan /3 is a key parameter in our analysis. The numerical analysis that we perform in the 
next section does not depend in any essential way on the simplification A5 = Afi. ()2.8|l . In 
essence we have just traded A5 for m^o in these rules and so by varying with respect to 
rriAo we do explore most of the quantitative potential of the general 2HDM. 

Since we shall perform our calculation in the on-shell scheme, we understand that the 
physical inputs are given by the electromagnetic coupling and the physical masses of all 
the particles: 

{e,Mw,Mz,mho,mHO,mAo,mH±,mf) . (4.4) 

The remaining parameters, except the Higgs mixing angles, are understood to be given in 
terms of the latter, e.g. the SU (2) gauge coupling appearing in the previous formulae and 
in Table ITTl is given hj g = e/s^ , where the sinus of the weak mixing angle is defined 
through = 1 — M^/M|. It should be clear that, as there are no tree-level FCNC decays 
of the top quark, there is no need to introduce counterterms for the physical inputs in 
this calculation. In fact, the calculation is carried out in lowest order ("tree level") with 
respect to the effective tch and teg couplings and so the sum of all the one-loop diagrams 
(as well as of certain subsets of them) should be finite in a renormalizable theory, and 
indeed it is. 
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Figure 4.1: One-loop vertex diagrams contributing to the FCNC top quark decays ()4.1|) . 
Shown are the vertices and mixed self-energies with all possible contributions from the SM 
fields and the Higgs bosons from the general 2HDM. The Goldstone boson contributions 
are computed in the Feynman gauge. 



4.3 Numerical analysis 

Prom the previous interaction Lagrangians and Feynman rules it is now straightforward 
to compute the loop induced FCNC rates for the decays ()4.H1 and ()4.2|) . We shall refrain 
from listing the lengthy analytical formulae as the computation is similar to the one 
reported in great detail in Ref. [43]. Therefore, we will limit ourselves to exhibit the final 
numerical results. The fiducial ratio on which we will apply our numerical computation 
is the following: 
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for each Type j = I, II of 2HDM and for each neutral Higgs boson h = h^, H^, A^. While 
this ratio is not the total branching fraction, it is enough for most practical purposes and 
it is useful in order to compare with previous results in the literature. Notice that for 
mfji > (the most probable situation for Type II 2HDM's, see below) the ratio (j4.5|) 
reduces to B^{t ^ h + c) = P{t ^ h + c)/T{t ^ W+ + b), which is the one that 
we used in Ref. [43]. It is understood that r^{t ^ h + c) above is computed from the 
one-loop diagrams in Fig. 14.11 with all quark families summed up in the loop. Therefore, 
consistency in perturbation theory requires to compute r{t —>■ + b) and r{t ^ + h) 
in the denominator of ()4.5|1 only at the tree- level (for explicit expressions see e.g. [20]). 
As mentioned in Sec. 2, we wish to compare our results for the Higgs channels ()4.1|) with 
those for the gluon channel ()4.2|) . so that we similarly define 

We have performed a fully-fledged independent analytical and numerical calculation of 
r'' (t — i> (yf + c) at one- loop in the context of 2HDM I and II. Where there is overlapping, 
we have checked the numerical results of Ref. [24] , but we point out that they agree with 
us only if r(t — s> II~^ + b) is included in the denominator of eq. ()4.6p . in contrast to what 
is asserted in that reference in which B{t ^ g + c) is defined without the charged Higgs 
channel contribution. 

We have performed part of the analytical calculation of the diagrams for both pro- 
cesses ()4.1|) and ()4.2|) by hand and we have cross-checked our results with the help of 
the numeric and algebraic programs FeynArts, FormCalc and Loop Tools [128-130], with 
which we have completed the rest of the calculation. In particular, the cancellation of 
UV divergences in the total amplitudes was also verified by hand. In addition we have 
checked explicitly the gauge invariance of the total analytical amplitude for the pro- 
cess ()4.2p . which is a powerful test. And we have confirmed that our code reproduces the 
SUSY Higgs contribution of Ref. [43] when we turn on the MSSM Higgs mass relations. 

As mentioned above, a highly relevant parameter is tan j3, which must be restricted 
to the approximate range ()2.43p in perturbation theory'^. It is to be expected from the 
various couplings involved in the processes under consideration that the low tan j3 region 
could be relevant for both the Type I and Type II 2HDM's. In contrast, the high tan/3 
region is only potentially important for the Type II. However, the eventually relevant 



■^Some authors [131] claim that perturbativity allows tan/3 to reach values of order 100 and beyond, 
and these are still used in the literature. We consider it unrealistic and we shall not choose tan /3 outside 



the interval (|2.43|l . Plots versus tan/3, however, will indulge larger values just to exhibit the dramatic 
enhancements of our FCNC top quark rates. 
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regions of parameter space are also determined by the value of the mixing angle a, as we 
shall see below. 

Of course there are several restrictions that must be respected by our numerical anal- 
ysis, see 12.41 as the p-parameter restriction, (j2.44j) . and the b —>■ s'-f restriction, (j2.45p . 
With the charge Higgs boson mass restrictions in section IT^ in principle the top quark 
decay t + 6 is still possible in 2HDM I; but also in 2HDM II, if m^^i lies near 

the lowest end of the previous bound, and in this case that decay can contribute to the 
denominator of eqs. ()4.5j) - ()4.(jj) . 

The combined set of independent conditions turns out to be quite effective in narrowing 
down the permitted region in the parameter space, as can be seen in Figs. I4.2H4.5I where 
we plot the fiducial FCNC rates ()4.5p - ()4.6|) versus the parameters ()4.3|) . The cuts in some 
of these curves just reflect the fact that at least one of these conditions is not fulfilled. 

After scanning the parameter space, we see in Figs. I4.2M.3I that the 2HDM I (resp. 
2HDM II) prefers low values (resp. high values) of tana and tan/9 for a given channel, 
e.g. t c. Therefore, the following choice of mixing angles will be made to optimize 
the presentation of our numerical results: 

2HDM I : tan a = tan /3 = 1/4 ; 

2HDM II : tan a = tan /? = 50 . (4.7) 

We point out that, for the same values of the masses, one obtains the same maximal 
FCNC rates for the alternative channel t c provided one just substitutes a — > 

Ti/2 — a. Equations (|4.7p define the eventually relevant regions of parameter space and, 
as mentioned in section 12.41 depend on the values of the mixing angles a and /3, namely 
/? ^ a ^ for Type I and /3 ~ a ^ 7r/2 for Type II. 

Due to the a — > 7r/2 — a symmetry of the maximal rates for the CP-even Higgs 
channels, it is enough to concentrate the numerical analysis on one of them, but one has 
to keep in mind that the other channel yields the same rate in another region of parameter 
space. Whenever a mass has to be fixed, we choose conservatively the following values 
for both models: 

niho = 100 GeV , mno = 150 GeV , m^o = mH± = 180 GeV . (4.8) 

Also for definiteness, we take the following values for some relevant SM parameters in our 
calculation: 

mt = 175 GeV , = 5 GeV , a,(mi) = 0.11 , Kb = 0.040 , (4.9) 

*At the time of this work this restriction was mfj± > (165 — 200) GeV for virtually any tan/3 > 1 
[60,80]. 
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Figure 4.2: Evolution of the FCNC top quark fiducial ratios (g^-dlSl) in Type I 2HDM 
versus: (a) the mixing angle a in the CP-even Higgs sector, in units of vr; (b) tan/3. The 
values of the fixed parameters are as in eqs. ()4.7|) and ()4.8|) . 



and the remaining ones are as in [75]. Notice that our choice of m^o prevents the decay 
t —>■ c from occurring, and this is the reason why it does not appear in Figs. 14.21 - 14.31 
The variation of the results with respect to the masses is studied in Figs. I4.4M.5I In 
particular, in Fig. 14.41 we can see the (scanty) rate of the channel t ^ c when it is 
kinematically allowed. In general the pseudoscalar channel is the one giving the skimpiest 
FCNC rate. This is easily understood as it is the only one that does not have trilinear 
couplings with the other Higgs particles (Cf. Table l^?T|) . While it does have trilinear 
couplings involving Goldstone bosons, these are not enhanced. The crucial role played by 
the trilinear Higgs self-couplings in our analysis cannot be underestimated as they can 
be enhanced by playing around with both (large or small) tan j3 and also with the mass 
splittings among Higgses. This feature is particularly clear in Fig. 14.4b where the rate of 
the channel t ^ c is dramatically increased at large m^o, for fixed values of the other 
parameters and preserving our list of constraints. Similarly would happen for t c 
in the corresponding region a ti /2 — a. 

From Figs. 14.2b and 14.2b it is pretty clear that the possibility to see FCNC decays 
of the top quark into Type I Higgs bosons is plainly hopeless even in the most favorable 
regions of parameter space - the lowest (allowed) tan/3 end. In fact, the highest rates 
remain neatly down 10~^, and therefore they are (at least) one order of magnitude below 
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Figure 4.3: As in Fig. 14.21 but for the 2HDM II. The plot in (b) continues above the 
bound in eq. ()2.43|) just to better show the general trend. 



the threshold sensibility of the best high luminosity top quark factory in the foreseeable 
future (see Section 4). We remark, in Fig. 14.21 that the rate for the reference decay 
t — (7 c in the 2HDM I is also too small but remains always above the Higgs boson rates. 
Moreover, for large tan (3 one has, as expected, {t g c) ^ j^sm ^ g ^-^ ^ 4 x 10"^^ 
because in this limit all of the charged Higgs couplings in the 2HDM I (the only Higgs 
couplings involved in this decay) drop off. Due to the petty numerical yield from Type I 
models we refrain from showing the dependence of the FCNC rates on the remaining 
parameters. 

Fortunately, the meager situation just described does not replicate for Type II Higgs 
bosons. For, as shown in Figs. 14.3k and l4.3b . the highest potential rates are of order 10~^, 
and so there is hope for being visible. In this case the most favorable region of parameter 
space is the high tan (3 end in eq. ()2.43p . Remarkably, there is no need of risking values over 
and around 100 (which, as mentioned above, are sometimes still claimed as perturbative!) 
to obtain the desired rates. But it certainly requires to resort to models whose hallmark 
is a large value of tan/5 of order or above mt/nih > 35. As for the dependence of the 
FCNC rates on the various Higgs boson masses (Cf. Figs. 14.4H4.5|) we see that for large 
rrij^o the decay t c can be greatly enhanced as compared to t g c; and of course, 
once again, the same happens with t c in the alternative region a ^ n /2 — a. We 

also note (from the combined use of Figs. 14.3b . 14.4b . and 14.4b ) that in the narrow range 
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Figure 4.4: Evolution of the FCNC top quark fiducial ratios (g^-diSl) in Type II 2HDM 
versus: (a) the CP-odd Higgs boson mass m^o; (b) the charged Higgs boson mass m//±. 
The values of the fixed parameters are as in eqs. ()4.7p and ()4.8|) . The plot in (b) starts 
below the bound m/^± > 165 GeV mentioned in the text to better show the general trend. 



where t b could still be open in the 2HDM II, the rate oi t h c becomes the 

more visible the larger and larger is tan (3 and m^o. Indeed, in this region one may even 
overshoot the 10"'^ level without exceeding the upper bound ()2.4Hj] while also keeping 
under control the remaining constraints, in particular eq. ()2.44j) . Finally, the evolution 
of the rates ()4.5p - ()4.6|) with respect to the two CP-even Higgs boson masses is shown in 
Figs. 14.5b and 14.5b . The neutral Higgs bosons themselves do not circulate in the loops 
(Cf. Fig. 14. ip but do participate in the trilinear couplings (Cf. Table 12. ip and so the 
evolution shown in some of the curves in Fig. 14.51 is due to both the trilinear couplings 
and to the phase space exhaustion. 

Turning now to the light scalar and pseudoscalar corners in parameter space mentioned 
above, it so happens that, after all, they prove to be of little practical interest in our case. 
Ultimately this is due to the quadratic Higgs boson mass differences entering 6p which 
make very difficult to satisfy the bound (j2.44j) . The reason being that for Type II models 
the limit ^ 165 GeV from 6 — > 57 implies that the constraint ()2.44j) cannot be 

preserved in the presence of light neutral Higgses. In actual fact the analysis shows that 
if e.g. one fixes rriho = 20 — 30 GeV, then the minimum m^o allowed by 6p is 100 GeV and 
the maximum rate (|4.5p is of order 10"^. Conversely, if one chooses rrij^o = 20 — 30 GeV, 
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Figure 4.5: As in Fig. 14.41 but plotting versus: (a) the lightest CP-even Higgs boson mass 
m/jo; (b) the heaviest CP-even Higgs boson mass rriffo. 



then the minimum m/jO allowed by 6p is 120 GeV and the maximum rate ()4.5|) is near 
10~^. Although in the last case the maximum rate is higher than in the first case, it is 
just of the order of the maximum rate already obtained outside the light mass corners of 
parameter space. On the other hand, these light mass regions do not help us in Type I 
models either. Even though for these models we do not have the b ^ S'-f bound on the 
charged Higgs, we still have the direct LEP 200 bound mfj± > 78.7 GeV [132] which is of 
course weaker than the CLEO bound. As a consequence the 6p constraint can be satisfied 
in the 2HDM I for neutral Higgs bosons lighter than in the corresponding 2HDM H case, 
and one does get some enhancement of the FCNC rates. Specifically, one may reach up 
to 10"^. However, the maximal rates ()4.5|) for the 2HDM I Higgs bosons are so small (see 
Figs. I4.2b - I4.2b ) that this order of magnitude enhancement is rendered immaterial. The 
upshot is that the top quark FCNC processes are not especially sensitive to the potential 
existence of a very light Higgs boson in either type of 2HDM. 



4-4 Discussion and conclusions 



59 



4.4 Discussion and conclusions 

The sensitivities to FCNC top quark decays for 100 of integrated luminosity in the 
relevant colliders are estimated to be [133]: 

LHC -.Bit ^ cX) > 5 X 10-^ 
LC -.Bit ^ cX) > 5 X 10"^ , (4.10) 
TEV33 ■.B{t cX) > 5 X 10"^ . 

This estimation has been confirmed by a full signal-background analysis for the hadron 
colliders in Ref. [25]. From these experimental expectations and our numerical results 
it becomes patent that whilst the Tevatron will remain essentially blind to this kind of 
physics, the LHC and the LC will have a significant potential to observe FCNC decays 
of the top quark beyond the SM. Above all there is a possibility to pin down top quark 
decays into neutral Higgs particles, eq. ()4.H) . within the framework of the general 2HDM 11 
provided tan (3 > mt/mh ~ 35. The maximum rates are of order 10~^ and correspond to 
the two CP-even scalars. This conclusion is remarkable from the practical (quantitative) 
point of view, and also qualitatively because the top quark decay into the SM Higgs 
particle is, in notorious contradistinction to the 2HDM 11 case, the less favorable top 
quark FCNC rate in the SM. On the other hand, we deem practically hopeless to see 
FCNC decays of the top quark in a general 2HDM I for which the maximum rates are of 
order 10"''. This order of magnitude cannot be enhanced unless one allows tan/? ^ 0.1, 
but the latter possibility is unrealistic because perturbation theory breaks down and 
therefore one cannot make any prediction within our approach. 

We have made a parallel numerical analysis of the gluon channel t —* c g in both types 
of 2HDM's. We confirm that this is another potentially important FCNC mode of the top 
quark in 2HDM extensions of the SM [24] but, unfortunately, it still falls a bit too short 
to be detectable. The maximum rates for this channel lie below 10~^ in the 2HDM 1 (for 
tan/3 > 0.1) and in the 2HDM 11 (for tan/5 < 60), and so it will be hard to deal with it 
even at the LHC. 

We are thus led to the conclusion that the Higgs channels ()4.1|) . more specifically the 
CP-even ones, give the highest potential rates for top quark FCNC decays in a general 
2HDM II. Most significant of all: they are the only FCNC decay modes of the top quark, 
within the simplest renormalizable extensions of the SM, that have a real chance to be 
seen in the next generation of high energy, high luminosity, colliders. 

The former conclusions are similar to the ones derived in Ref. [43] for the MSSM case, 
but there are some conspicuous differences on which we wish to elaborate a bit in what 
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follows [134]. First, in the general 2HDM II the two channels t {h^, H^) c give the same 
maximum rates, provided we look at different (disjoint) regions of the parameter space. 
The t ^ c channel is, as mentioned, negligible with respect to the CP-even modes. 
Hereafter we will discard this FCNC top quark decay mode from our discussions within the 
2HDM context. On the other hand, in the MSSM there is a most distinguished channel, 
viz. t — > c, which can be high-powered by the SUSY stuff all over the parameter 
space. In this framework the mixing angle a becomes stuck once tan /3 and the rest of the 
independent parameters are given, and so there is no possibility to reconvert the couplings 
between and as in the 2HDM. Still, we must emphasize that in the MSSM the other 
two decays t c and t c can be competitive with t c in certain portions 

of parameter space. For example, t — > H'^ c becomes competitive when the pseudoscalar 
mass is in the range 110 GeV < rrij^o < 170 GeV [43]. The possibility of having more 
than one FCNC decay (|4.H) near the visible level is a feature which is virtually impossible 
in the 2HDM II. Second, the reason why t ^ c in the MSSM is so especial is that it 
is the only FCNC top quark decay ()4.H) which is always kinematically open throughout 
the whole MSSM parameter space, while in the 2HDM all of the decays 1)4.11) could be, 
in the worse possible situation, dead closed. Nevertheless, this is not the most likely 
situation in view of the fact that all hints from high precision electroweak data seem to 
recommend the existence of (at least) one relatively light Higgs boson [132,134]. This 
is certainly an additional motivation for our work, as it leads us to believe that in all 
possible (renormalizable) frameworks beyond the SM, and not only in SUSY, we should 
expect that at least one FCNC decay channel ()4.H) could be accessible. Third, the main 
origin of the maximum FCNC rates in the MSSM traces back to the tree-level FCNC 
couplings of the gluino [43]. These are strong couplings, and moreover they are very 
weakly restrained by experiment. In the absence of such gluino couplings, or perhaps by 
further experimental constraining of them in the future, the FCNC rates in the MSSM 
would boil down to just the electroweak (EW) contributions, to wit, those induced by 
charginos, squarks and also from SUSY Higgses. The associated SUSY-EW rate is of 
order 10~^ at most, and therefore it is barely visible, most likely hopeless even for the 
LHC. In contrast, in the general 2HDM the origin of the contributions is purely EW and 
the maximum rates are two orders of magnitude higher than the full SUSY-EW effects in 
the MSSM. It means that we could find ourselves in the following situation. Suppose that 
the FCNC couplings of the gluino get severely restrained in the future and that we come 
to observe a few FCNC decays of the top quark into Higgs bosons, perhaps at the LHC 
and/or the LC. Then we would immediately conclude that these Higgs bosons could not 
be SUSY-MSSM, whilst they could perhaps be CP-even members of a 2HDM II. Fourth, 
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the gluino effects are basically insensitive to tan/3, implying that the maximum MSSM 
rates are achieved equally well for low, intermediate or high values of tan/5, whereas the 
maximum 2HDM II rates (comparable to the MSSM ones) are attained only for high 
tan j3. 

The last point brings about the following question: what could we possibly conclude if 
the gluino FCNC couplings were not further restricted by experiment and the tagging of 
certain FCNC decays of the top quark into Higgs bosons would come into effect? Would 
still be possibly to discern whether the Higgs bosons are supersymmetric or not? The 
answer is, most likely yes, provided certain additional conditions would be met. 

There are many possibilities and corresponding strategies, but we will limit ourselves 
to point out some of them. For example, let us consider the type of signatures involved in 
the tagging of the Higgs channels. In the favorite FCNC region fl4.7|l of the 2HDM II, the 
combined decay t ^ h c ^ ebb is possible only for or for H^, but not for both - Cf. 
Fig. 14.3b - whereas in the MSSM, together with H'^, are highlighted for 110 GeV < 
rrLj^o < rrit, with no preferred tan (3 value. And similarly, t — > A" c is also non-negligible for 
nij^o < 120 GeV [43]. Then the process t h c ebb gives rise to high px charm-quark 
jets and a recoiling hb pair with large invariant mass. It follows that if more than one 
distinctive signature of this kind would be observed, the origin of the hypothetical Higgs 
particles could not probably be traced back to a 2HDM II. 

One might worry that in the case of and they could also (in principle) decay 
into electroweak gauge boson pairs hP,H^ VewVew, which in some cases could be 
kinematically possible. But this is not so in practice for the 2HDM II if we stick to 
our favorite scenario, eq. (|4.7p . In fact, we recall that the decay — > VewVew is not 
depressed with respect to the SM Higgs boson case provided (3 — a = tx /2, and similarly 
for VewVew a (3 — a = 0. However, neither of these situations is really pinpointed 

by FCNC physics because we have found (3 ~ 7r/2 in the most favorable region of our 
numerical analysis, and moreover a was also seen there to be either a ~ 7r/2 (for h^) or 
(for H^), so both decays — > VewVew are suppressed in the regions where the FCNC 
rates of the parent decays t — > (/i°, H^) e are maximized. Again, at variance with this 
situation, in the MSSM case VewVew is perfectly possible - not so h'^ VewVew 

due to the aforementioned upper bound on niho - because tan (3 has no preferred value in 
the most favorable MSSM decay region of t ^ e. Therefore, detection of a high pj- 
charm-quark jet against a VewVew pah of large invariant mass could only be advantageous 
in the MSSM, not in the 2HDM. Similarly, for tan /5 > 1 the decay H° /i° h° (with real 
or virtual h^) is competitive in the MSSM [135, 136] in a region where the parent FCNC 
top quark decay is also sizeable. Again this is impossible in the 2HDM II and therefore 
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it can be used to distinguish the two (SUSY and non-SUSY) Higgs frames. 

Finally, even if we place ourselves in the high tan (3 region both for the MSSM and the 
2HDM II, then the two frameworks could still possibly be separated provided that two 
Higgs masses were known, perhaps one or both of them being determined from the tagged 
Higgs decays themselves, eq. ()4.H1 . Suppose that tan/9 is numerically known (from other 
processes or from some favorable fit to precision data), then the full spectrum of MSSM 
Higgs bosons would be approximately determined (at the tree level) by only knowing 
one Higgs mass, a fact that could be used to check whether the other measured Higgs 
mass becomes correctly predicted. Of course, the radiative corrections to the MSSM 
Higgs mass relations can be important at high tan/? [114], but these could be taken into 
account from the approximate knowledge of the relevant sparticle masses obtained from 
the best fits available to the precision measurements within the MSSM. If there were 
significant departures between the predicted mass for the other Higgs and the measured 
one, we would probably suspect that the tagged FCNC decays into Higgs bosons should 
correspond to a non-supersymmetric 2HDM II. 

At the end of the day we see that even though the maximum FCNC rates for the 
MSSM and the 2HDM II are both of order 10~^ - and therefore potentially visible - at 
some point on the road it should be possible to disentangle the nature of the Higgs model 
behind the FCNC decays of the top quark. Needless to say, if all the fuss at CERN [132] 
about the possible detection of a Higgs boson would eventually be confirmed, this could 
still be interpreted as the discovery of one neutral member of an extended Higgs model. 
Obviously the combined Higgs data from LEP 200 and the possible discovery of FCNC 
top quark decays into Higgs bosons at the LHC/LC would be an invaluable cross-check 
of the purportedly new phenomenology. 

We emphasize our most essential conclusions in a nutshell: i) Detection of FCNC top 
quark decay channels into a neutral Higgs boson would be a blazing signal of physics 
beyond the SM; ii) There is a real chance for seeing rare events of that sort both in 
generic Type II 2HDM's and in the MSSM. The maximum rates for the leading FCNC 
processes fl4.1|) and ()4.2j) in the 2HDM II (resp. in the MSSM) satisfy the relations 

B{t ^gc)< 10~^(10^^) <B{t^hc)^ 10"^ , (4.11) 

where it is understood that h is or but not both, in the 2HDM II; whereas h is 
most likely but it could also be and in the MSSM ; iii) Detection of more than 
one Higgs channel would greatly help to unravel the type of underlying Higgs model. 

The pathway to seeing new physics through FCNC decays of the top quark is thus 
potentially open. It is now an experimental challenge to accomplish this program using 
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the high luminosity super-colhders round the corner. 



Chapter 5 



Higgs Boson FCNC Decays into Top 
Quark in a General 2HDM 

5.1 Introduction 

When considering physics beyond the SM, new horizons of possibihties open up which 
may radically change the pessimistic prospects for FCNC decays involving a Higgs boson 
and the top quark. For example, in Ref. [42] it was shown that the vector boson modes 
can be highly enhanced within the context of the Minimal Supersymmetric Standard 
Model (MSSM) [7-10]. This fact was also dealt with in great detail in Ref. [43] where in 
addition a dedicated study was presented of the FCNC top quark decays into the various 
Higgs bosons of the MSSM (see also [44]), showing that these can be the most favored 
FCNC top quark decays - above the expectations on the gluon mode t ^ eg. A similar 
study is performed in the chapter |^ for the FCNC top quark decays into Higgs bosons in 
a general two-Higgs-doublet model (2HDM). 

In the previous chapter analysing the FCNC top quark decays in 2HDM extensions 
of the SM it was proven that while the maximum rates for t eg were one order of 
magnitude more favorable in the MSSM [43] than in the 2HDM, the corresponding rates 
for t eh^ were comparable both for the MSSM and the general 2HDM, namely up to 
the 10~^ level and should therefore be visible both at the LHC and the LC [25]. 

Similarly, one may wonder whether the FCNC decays of the Higgs bosons themselves 
can be of some relevance. Obviously the situation with the SM Higgs is essentially hope- 
less, so again we have to move to physics beyond the SM. Some work on these decays, 
performed in various contexts including the MSSM, shows that these effects can be im- 
portant [47,48,137,138], as seen in chapter IHl This could be expected, at least for 
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heavy quarks in the MSSM, from the general SUSY study (including both strong and 
electroweak supersymmetric effects) of the FCNC vertices htc {h = ,A^) made 

in Ref. [43]. However, other frameworks could perhaps be equally advantageous. Here 
we are particularly interested in the FCNC Higgs decay modes into top quark within a 
general 2HDM, which have not been studied anywhere in the literature to our knowledge. 
It means we restrict to Higgs bosons heavier than m^. From the above considerations, 
and most particularly on the basis of the detailed results obtained in the previous chapter 
m one may expect that some of the decays of the Higgs bosons 

h^tc, h^tc {h = h\H\A^) (5.1) 

in a general 2HDM can be substantially enhanced and perhaps can be pushed up to 
the visible level, particularly for which is the lightest CP-even spinless state in these 
models [6]. This possibility can be of great relevance on several grounds. On the one hand 
the severe degree of suppression of the FCNC Higgs decay in the SM obviously implies 
that any experimental sign of Higgs-like FCNC decay (jS.H) would be instant evidence 
of physics beyond the SM. On the other hand, the presence of an isolated top quark 
in the final state, unbalanced by any other heavy particle, is an unmistakable carrier of 
the FCNC signature. Finally, the study of the maximum FCNC rates for the top quark 
modes ()5.1|) within the 2HDM, which is the simplest non-trivial extension of the SM, 
should serve as a fiducial result from which more complicated extensions of the SM can 
be referred to. Therefore, we believe there are founded reasons to perform a thorough 
study of the FCNC Higgs decays in minimal extensions of the Higgs sector of the SM and 
see whether they can be of any help to discover new physics. 

The study is organized as follows. In Section E21 we summarize the 2HDM interactions 
most relevant for our study and estimate the expected FCNC rates of the Higgs decays 
in the SM and the general 2HDM. In Section 15.31 a detailed numerical analysis of the 
one-loop calculations of the FCNC decay widths and production rates of FCNC Higgs 
events is presented. Finally, in Section 15.41 we discuss the reach of our results and its 
phenomenological implications, and deliver our conclusions. 

5.2 Expected branching ratios in the SM and the 
2HDM 

Before presenting the detailed numerical results of our calculation, it may be instructive 
to estimate the typical expected widths and branching ratios (B) both for the SM decay 
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non-standard decays (jS-H) in a general 2HDM. This should be 
especially useful in this kind of rare processes, which in the strict context of the SM are 
many orders of magnitude out of the accessible range. Therefore, one expects to be able 
to grossly reproduce the order of magnitude from simple physical considerations based 
on dimensional analysis, power counting, CKM matrix elements and dynamical features. 
By the same token it should be possible to guess at the potential enhancement factors in 
the 2HDM extension of the SM. In fact, guided by the previous criteria the FCNC decay 
width of the SM Higgs of mass mn into top quark is expected to be of order 

r(i7^*^^tc)~ |l^;,Hcr«3v.mH(An'~ (l^y aH^G|mHm^ (5.2) 

where Gp is Fermi's constant and aw = (7^/47r, g being the SU {2)l weak gauge coupling. 
We have approximated the loop form factor by just a constant prefactor. Notice the 
presence of Af*^ ~ rrib/Mw, which is the SM Yukawa coupling of the bottom quark in 
units of g. The fourth power of Af *^ in ()5.2|) gives the non-trivial suppression factor rem- 
iniscent of the GIM mechanism after summing over flavors. Since we are maximizing our 
estimation, a missing function related to kinematics and polarization sums, F{mt/mH) ~ 
(1— m^/m|^)^, has been approximated to one. To obtain the (maximized!) branching ratio 
it suffices to divide the previous result by T{H^'^^ — > bb) ~ aw (Af*^)^ mn ~ GpTnuml 
to obtain 

B{H'^' ^tc)^(^^^ awGFml^lQ'^\ (5.3) 

with Vbc = 0.04, mb = 5GeV. In general this B will be even smaller, specially for higher 
Higgs boson masses {rriH > 2 Mw) for which the vector boson Higgs decay modes H^'^^ 

W~{Z Z) can be kinematically available and become dominant. In this case it is easy 
to see that B{H^^^ ^ tc) will be suppressed by an additional factor of ml/nijj, which 
amounts at the very least to two additional orders of magnitude suppression, bringing it to 
a level of less than 10^^^. Already the optimized branching ratio ()5.3p will remain invisible 
to all foreseeable accelerators in the future! To obtain the corresponding maximized 
estimation for the 2HDM we use the couplings from (j2.42j) and the trilinear couplings 12. 11 

Let us now first assume large tan (3 and restrict to Type II models. From the interaction 
Lagrangians above it is clear that we may replace Af — > A^^^ tan (3 in the previous 
formulae for the partial width. Moreover, the leading diagrams in the 2HDM contain the 
trilinear Higgs couplings Xn+H-h- Therefore, the maximum B associated to the FCNC 
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decays ()5.1|) in a general 2HDM II should be of order^ 



B^\h ^ tc) i^j^^ aw Gprnltw? [3 (5.4) 

where \h+ h- h is defined here in units of g and dimensionless as compared to Table 12.11 
Clearly a big enhancement factor tan^ jS appears, but this does not suffice. Fortunately, 
the trilinear couplings Xh+ h- h for h = h?,H^ (but not for h = A^) carry two additional 
sources of potential enhancement (Cf. Table 12.11) which are absent in the MSSM case. 
Take e.g. then we see that under appropriate conditions (for example, large tana and 
large tan /5) the trilinear coupling behaves as Xh+ h- ~ ("^^o — 'm\o) taxi(3/{Mw mjj±), 
and in this case 

B'\h' -^tc)^ f ' aw ml tan^ (3 i'^^^) ' • (5.5) 

So finally B^\h^ tc), and of course B^\h^ — * tc), can be augmented by a huge factor 
tan^ /5 times the square of the relative splitting among the CP-even Higgs decaying boson 
mass and the CP-odd Higgs mass. Since the neutral Higgs bosons do not participate in 
the loop form factors fsee 14.11) . it is clear that various scenarios can be envisaged where 
these mass splittings can be relevant. In the next section this behaviour will be borne out 
by explicit calculations showing that —* tc can be raised to the visible level in the case 
of the Type II model. As for the Type I model the Higgs trilinear coupling enhancement 
is the same, but in the charged Higgs Yukawa coupling all quarks go with a factor cot j3] 
hence when considering the leading terms in the loops that contribute one sees that in 
the corresponding expression (j5.4j) the term ml tan^ (3 is traded for m^ cot^ j3, which 
is negligible at high tan (3. Both sources of enhancement are needed, and this feature is 
only tenable in the 2HDM II. Of course one could resort to the range tan /9 -C 1 for the 
Type I models, but this is not theoretically appealing. For example, for tan/3 < 0.1 the 
top quark Yukawa couphng gt = g mt/{V2 Mw sin/3), which is present in the interaction 
Lagrangians above, is pushed into the non-perturbative region gf/lQir"^ > 1 and then 
our calculation would not be justified. And what is worse: for the 2HDM I we would 
actually need tan/3 < O{10~^) to get significant FCNC rates! In short, we consider that 
B^{h tc+ tc) is essentially small (for all h), and that these decays remain always 
invisible to speak of. Hereafter we abandon the study of the decays (|5.1|) for the 2HDM 
I and restrict ourselves to the general 2HDM II. 



""^Here we have normalized the B with respect to the h ^ bb channel only, because the gauge boson 
modes will be suppressed in the relevant FCNC region. Cf. Section 3. 
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5.3 Numerical analysis 

Let us now substantiate the previous claims and provide the precise numerical results of 
the full one-loop calculation of B^\h ^ tc + ic)"^ as well as of the LHC production rates 
of these FCNC events. We refer the reader to previous chapters for more details. The 
diagrams for the decays ()5.1|) are shown in Fig. 15.11 and the diagrams for the productions 
can be seen in Fig. 15.21 In what follows we present the final results of our numerical 
analysis together with a detailed discussion, interpretation and phenomenological appli- 
cation. We have performed the calculations with the help of the numeric and algebraic 
programs FeynArts, FormCalc and Loop Tools [128-130]. The calculation must obviously 
be finite without renormalization, and indeed the cancellation of UV divergences in the 
total amplitudes was verified explicitly. 

The input set for our numerical analysis is given by the data row 

[rriho , rriHO , m AO , rriHi , tan a, tan/5) (5.6) 

made out of six independent parameters in the general 2HDM. Remaining inputs as 
in [75]. In practice there are some phenomenological restrictions on the data ()5.6p which 
were already described in the chapter |2l particularly [79]. Again, a key parameter is 
tan/5, with the restriction (j2.43p . In practice, since Type I models are not considered, the 
effective range for our calculation will be the high tan/3 end of ()2.43|) . 

With these restrictions in mind we have computed the number of FCNC Higgs decay 
events into top quark at the LHC: 

pp^ h + X ^tc{ic)+X {h = h\H^,A°) . (5.7) 

The necessary cross-sections to compute the production of neutral Higgs bosons at this 
collider, including all known QCD corrections, have been computed by adapting the codes 
HIGLU 1.0 and HQQ 1.0 [139] - originally written for the MSSM case [140]- to the general 
case of the 2HDM'^. Folding the cross-sections with the one-loop branching ratios of the 
processes (j5.H) we have obtained the number of FCNC Higgs decay events at the LHC. Let 
us first consider the branching ratios themselves. In Fig. l5.3b .b we show B^^{hP t c+t c) 
for the lightest CP-even state. In particular. Fig. 15.3b shows B^\h^ tc + ic) versus 
the charged Higgs mass mH±. In this figure we fix the values of the parameters in ()5.6|) 

^Here and throughout we use the notation B" ^tc + tc) = B"{h^ ^tc) + B"{h° -^tc). 

^We have used the default parton distribution functions and renormahzation/factorization scales used 
in these programs, namely GRV94 with nn ~ jip — rrih for HIGLU, and CTEQ4L with nn = jip — \/I = 
^{Ph+PQ+PqY for HQQ. 
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Figure 5.1: One-loop vertex diagrams contributing to the FCNC Higgs decay (jS.lj) . Shown 
are the vertices and mixed self-energies with all possible contributions from the SM fields 
and the Higgs bosons from the general 2HDM. The Goldstone boson contributions are 
computed in the Feynman gauge. 
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Figure 5.2: Leading order Feynman diagrams for the Higgs boson production at the LHC 
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Figure 5.3: (a) B^\h^ tc + ic) versus m/^±; (b) Idem, versus m^o; (c) The production 
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Figure 5.4: Contour lines in the (m//±, m/io) -plane for the branching ratios (2HDM II 
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Figure 5.5: (a) Contour lines in the (m//±, m/^o) -plane for the maximum number of light 
CP-even Higgs FCNC events h° tc + tc produced at the LHC for 100 fb~^ of inte- 
grated luminosity within 5p^^°^ at 1 a; (b) Contour lines showing the value of m^o that 
maximizes the number of events. 
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Figure 5.6: As in Fig. O but within ^p^™^ at 3 a. 
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Figure 5.7: Contour lines a/n = const, [a is the mixing angle in the CP-even sector) 
corresponding to Figs. 15.5115.61 for 5p^^°^ at (a) 1 a and (b) 3 a . 
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Number of events 




Figure 5.8: (a) Contour lines in the (m//±, mj^/o) -plane for the maximum number of heavy 
CP-even Higgs FCNC events tc + ic (2HDM II case) produced at the LHC for 

100 fb~^ of integrated luminosity within Sp^^^^ at 1 a. 
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Figure 5.9: As in Fig.EHlbut within 5p™^^ at 3 a. 
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Figure 5.10: Contour lines a/vr = const, as in Fig l5.71 but for the heavy CP-even Higgs. 
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which are not varying as follows: 

{rriho = 350 GeV, tuho = 600 GeV, m^o = 550 GeV, mj^± = 375 GeV, , , 

(5.8) 

tana = 30, tan /3 = 60) 

After crossing a local maximum (associated to a pseudo-threshold of the one-loop vertex 
function involving the H~ coupling) the subsequently falling behavior of the B with 

mH± clearly shows that the previously discussed bounds on mH± are quite relevant. The 
branching ratio, however, stays within 10~^ — 10~^ for a wide range of heavy charged Higgs 
masses extending up to < 600 GeV^ in Fig. 15.3b . Hence, for mH± heavy enough to 
satisfy the indirect bounds from radiative 5-meson decays [62], the maximum B is still 
sizeable. In Fig. 15. 3h the production rate of bosons at the LHC is shown as a function 
of for fixed parameters ()5.6|) . The production cross-sections for the subprocesses 

gg^ h^ + X, gg.qq^ + QQ , (5.9) 

contributing to ()5.7p in the case of the light CP-even Higgs are explicitly separated in 
Fig. 15.3b . The gluon-gluon fusion process proceeds at one-loop and the h^QQ associated 
production proceeds at tree-level [136,141]. Similar subprocesses and results apply for 
and production. At large tan (3 and the larger the Higgs boson masses the particular 
associated production mechanism with the bottom quark, Q = b, i.e. bb, becomes 
dominant by far. All other mechanisms for Higgs boson production in Type II models 
[31,136,140,141], like vector-boson fusion (which contributes also to h^QQ when Q are 
light quarks), vector-boson bremsstrahlung {qq ^ hV) and associated ti production, are 
subdominant at large tan P and can be neglected for our purposes. Admittedly, some of 
these mechanisms can be relevant for Higgs boson production in the case of the Type I 
2HDM at low tan /3, but we have already warned that the corresponding FCNC branching 
ratios are never sufficiently high. 

The control over Sp"^™^ is displayed in Fig. I5.3l i. Recall that Sp is not sensitive 
to the mass splitting between m/jO and m^^o, because of CP-conservation in the gauge 
boson sector, but it does feel all the other mass splittings among Higgs bosons, charged 
and neutral. A more systematic search of B values in the parameter space is presented in 
Figs. I5.4b .b corresponding to B^\h^ —* tc+ic) and B^\H^ —* tc+ic) respectively. Here 
we have scanned independently on the parameters (j5.6|) while holding the bound 
at 1 a. The contour lines in these figures represent the locus of points in the {mH±,mha)- 
plane giving maximized values of the B in the 2HDM II. Let us remark that the highest 
value of tan/? is always preferred, and therefore all these contour lines correspond to 
tan/5 = 60. 
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In practice, to better assess the possibility of detection at the LHC, one has to study the 
production rates of the FCNC events. These are determined by combining the production 
cross-sections of neutral 2HDM II Higgs bosons at the LHC and the FCNC branching 
ratios. If we just adopt the mild LEP bound > SOGeV and let mfj± approach the 
maximum in Fig. 15.3b then the B can be as large as 10"^ and the number of FCNC 
events can be huge, at the level of ten thousand per 100 fb~^ of integrated luminosity. 
But of course the region near the maximum is too special. Moreover, if we switch on the 
above mentioned indirect bound from b s'-f [62], then the typical B is much smaller 
(of order 10~^) and the number of events is reduced dramatically, at a level of hundred 
or less for the same integrated luminosity. On the other hand it may well happen that 
there are regions of parameter space where B ~ 10~^ (see Fig. 15. 4j) but the production 
cross-section is too small because the decaying Higgs boson is too heavy. Therefore, it is 
the product of the two quantities that matters. 

The systematic search of the regions of parameter space with the maximum number of 
FCNC events for the light CP-even Higgs is presented in the form of contour lines in the 
multiple Figs. I5.5] and l5.6( For instance, each isoline in Figs. 15.5k and l5.6b corresponds to 
a fixed number of produced FCNC events at the LHC while keeping the value of 6p^^^^ 
within 1 (T or 3 0" respectively of its central experimental value. When scanning over the 
parameter space ()5.6|) we have found again that tan /3 is preferred at the highest allowed 
value (tan (3 = 60) - for Type II models. We have also determined (see Figs. 15. 51 b and 
I5.6l d) the corresponding contour lines for m^o associated to these events. The m^o-lines 
are important because the FCNC processes under consideration are sensitive to the mass 
splittings between m^o and the corresponding decaying Higgs boson, see e.g. eq. (j5.5|) 
and Table 12.11 The combined figures 15.51 and 15.61 are very useful because they give a 
panoramic view of the origin of our results in the parameter space. To complete the map 
of the numerical analysis we provide Fig. 15. 71 in which we have projected the contour lines 
of the CP-even mixing angle a associated to the previous plots. For a given contour line 
a/vr = const., the set of inner points have a value of a/vr smaller than the one defined 
by the line itself. In particular, the large domains in Figs. I5.7b .b without contour lines 
correspond to a/ir > 0.4 and so to relatively large (and positive) tana. There are a few 
and small neighborhoods where the FCNC rates for can be sizeable also for small tana. 

Knowing that high tan a is generally preferred by tc + tc, and noting from 

Figs. 15.51 and 15. 61 that large mass splittings between m/jO and m^o are allowed, we find that 
the trilinear coupling \h+ h- /i" can take the form \h+ h- ho ~ i^ho """^^o ) tan P / {Mw fnH± ) 
Hence it provides a substantial additional enhancement beyond the tan (3 factor. One can 
check from the approximate formula (|5.5p that the maximum FCNC branching ratios 
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B^\h^ ^ tc + ic) can eventually reach the 10~^ level even in regions where the charged 
Higgs boson mass preserve the stringent indirect bounds from radiative i?-meson de- 
cays [62]. These expectations are well in agreement with the exact numerical analysis 
presented in Fig. 15.41 thus showing that eq. (j5.5p provides a reasonable estimate, and there- 
fore a plausible explanation for the origin of the maximum contributions. As a matter of 
fact, we have checked that the single (finite) Feynman diagram giving rise to the estima- 
tion ()5.5|) - the one-loop vertex Feynman diagram with a couple of charged Higgs bosons 
and a bottom quark in the loop - reproduces the full result with an accuracy better than 
10% for tan/5 > 10 — 20. At lower tan/3 values large deviations are possible but, as warned 
before, eq. ()5.5p is expected to be valid only at large tan j3. Furthermore, for low values 
of tan/3 < 20 the FCNC Bs are too small to be of any phenomenological interest. The 
exact numerical analysis is of course based on the full expression for the branching ratio 

jj Q _ - r(/i° — > t c + t c) 

(/, ^tc + tc) = ^^^^ ^ ^-^ ^ ^^^^ ^ ^ ^^^^ -.VV) + r(/iO -^HH)' ^^-^^^ 

where all decay widths in the denominator of this formula have been computed at the 
tree-level in the 2HDM II, since this provides a consistent description of eq. ()5.10|) at 
leading order. Here we have defined 

r(/i° ^vv) = r(/i° ^ w^w) + t{h^ ^ z z), (5.11) 



T{h° -^HH)= r(/i° A° A°) + r(/i° ^ H+ H-) . (5.12) 

We disregard the loop induced decay channels, since they have branching ratios below the 
percent level all over the parameter space. The r-lepton decay channel is also neglected, 
since it is suppressed by a factor of (9(10^^) with respect the 66-channel in the whole 2HDM 
parameter space. In general the effect of the gauge boson channels h*^ — ^ W^W~ , ZZ 
in the B (I5.10|) is not so important as in the SM, actually for j3 = a they vanish in the 
case because they are proportional to sin^(/3 — a). This is approximately the case 
for large tana and large tan/3, the dominant FCNC region for /i° decay (Cf. Fig J5.7b 
and 15.7b ). In this region, the mode ^ ti is, when kinematically allowed, suppressed: 
B{h? ^ti) oc cos^ a/ sin^ (3^0 (Cf. Eq. ^T^ ). On the other hand there are domains in 
our plots where the decays H~ and are kinematically possible and 

non- (dynamically) suppressed. Indeed, this can be checked from the explicit structure 
of the trilinear couplings H~ and A^ A^ in Table 12.11 in the dominant region 

for the decays tc + tc both of these couplings are tan /3-enhanced. Nevertheless 

the decay — ^ A^ A^ is only possible for m^o < m^o/2 , and since the optimal FCNC 
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regions demand the largest possible values of m^o, this decay is kinematically blocked 
there. On the other hand the mode H~ is of course allowed if we just take the 

aforementioned direct limits on the 2HDM Higgs boson masses. But it is never available 
if we apply the indirect bound from 6 ^ S7 on the charged Higgs mass mentioned above, 
unless m/jo > 2mu± > 700 GeV, in which case is so heavy that its production cross- 
section is too small for FCNC studies to be further pursued. 

The corresponding results for the heavy CP-even Higgs boson are displayed in Figs. 
15. 8| 15.91 and 15.101 The exact formula for the B in this case reads 

(5.13) 

where we have defined 

T{H° ^VV)= T{H° ^ W+W^) + r(if° ^ Z Z), (5.14) 



r(if° ^HH) = r(i7° ^ /i°) + r(ij° ^ a° a°) + r(/7° ^ h~) . (5.15) 

From the contour lines in Figs. 15. 8k and 15.9b it is patent that the number of FCNC top 
quark events stemming from decays is comparable to the case of the lightest Higgs 
boson. However, Fig. l5.1(Jb .b clearly reveals that these events are localized in regions of 
the parameter space generally different from the h*^ case, namely they prefer tana ~ 0. 
Even so, there are some "islands" of events at large tana. This situation is complementary 
to the one observed for in Fig. 15. 71 However, in both cases these isolated regions are 
mainly concentrated in the segment m//± < 350 GeV^. Therefore, if the bound on m//± 
from 6 — > s 7 is strictly preserved, it is difficult to find regions of parameter space where 
the two CP-even states of a general 2HDM II may both undergo a FCNC decay of the 
type (jHH). 

In the dominant regions of the FCNC mode ^ tc (where tana is small and 
tan/3 is large), the decay of into the ti final state is suppressed: B{H ^ tt) oc 
sin^ a/ sin^ /5 — 0. In the same regions the gauge boson channels in ()5.13p are suppressed 
too because r{H^ — > W~^W~, ZZ) oc cos^(/5 — a). In principle the heavy CP-even Higgs 
boson also could (as h^) decay into and H~ . But there is a novelty here 

with respect to the decays, in that there could be regions where could decay into 
the final state . This contingency has been included explicitly in eq. ()5.15|) . However, 
in practice, neither one of these three last channels is relevant in the optimal FCNC 
domains of parameter space. First, the decay h^, although it is kinematically 

possible, is dynamically suppressed in the main FCNC region for H'^. This can be seen 
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from Table 12.11 where the triUnear couphng becomes vanishingly small at large 

tan/? and small tana. Second, the coupling H^A^A^ in Table ITT] is non-suppressed in 
the present region, but again the mode H'^ — A^ A^ is kinematically forbidden in the 
optimal FCNC domains because the latter favor large values of the CP-odd mass (see 
Figs. 15. 8b and 15.9b ). Third, although in these domains the decay — > H~ is also 
non-dynamically suppressed (see the corresponding trilinear coupling in Table 12.11) . it 
becomes kinematically shifted to the high mass range rriijo > 700 GeV if we switch on the 
indirect bound from b s'-f. Obviously, in this latter case the production cross section 
becomes too small and the FCNC study has no interest. All in all the contributions from 
()5.11|) . ()5.12|) . ()5.14|) and ()5.15|) are irrelevant for myfi^mHO < 700 GeV as their numerical 
impact on B^^{h^, — > tc + tc) is negligible. Our formulae ()5.10p and ()5.13p do contain 
all the decay channels and we have verified explicitly these features. 

As remarked before, in general the most favorable regions of parameter space for the 
FCNC decays of and do not overlap much. The trilinear Higgs boson self-couplings 
in Table im falso the fermionic ones) are interchanged when performing the simultaneous 
substitutions a ixjl — a and m/^o rrifjo (see last chapter Furthermore, the 
LHC production rates of the neutral Higgs bosons fall quite fast with the masses of these 
particles, as seen e.g. in Fig. 15.3b for the state. As a consequence that exchange 
symmetry on the branching ratios does not go over to the final event rates, so in practice 
the number of FCNC events from decays are smaller (for the same values of the other 
parameters) as compared to those for h^; thus requires e.g. lighter charged Higgs 
masses to achieve the same number of FCNC events as h^. As for the CP-odd state 
we have seen that it plays an important indirect dynamical role on the other decays 
through the trilinear couplings in Table 12.11 but its own FCNC decay rates never get a 
sufficient degree of enhancement due to the absence of the relevant trilinear couplings, so 
we may discard it from our analysis. 

We notice that this picture is consistent with the decoupling limit in the 2HDM: 
for a (3, the heaviest CP-even Higgs boson (H^) behaves as the SM Higgs boson, 
whereas decouples from the electroweak gauge bosons and may develop enhanced 
couplings to up and down-like quarks, depending on whether tan/3 is small or large 
respectively; in the opposite limit {a (3 — tt /2) , it is that behaves as H^^^ , while H*^ 
decouples from gauge bosons and may develop the same enhanced couplings to quarks 
as did in the previous case. Indeed these are the situations that we find concerning 
the FCNC decay rates. We recall that the numerical results presented in our figures 
correspond to an integrated luminosity of 100 However, the combined ATLAS and 
CMS detectors might eventually accumulate a few hundred inverse femtobarn [13,14]. 
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Therefore, hopefully, a few hundred FCNC events ()5.1|) could eventually be collected in 
the most optimistic scenario. Actually, the extreme rareness of these events in the SM 
suggests that if only a few of them could be clearly disentangled, it should suffice to claim 
physics beyond the SM. 

5.4 Discussion and conclusions 

Detection strategies at the CERN-LHC collider for the search of the SM Higgs boson, 
and also for the three spinless fields of the MSSM Higgs sector, have been described in 
great detail in many places of the literature [6,11,13,14,17,142,143], but not so well 
for the corresponding charged and neutral Higgs bosons of the general 2HDM. The result 
is that the discovery of the SM Higgs boson is guaranteed at the LHC in the whole 
presumed range 100 GeV^ < mn ^ ITeV. However, the discovery channels are different 
in each kinematical region and sometimes the most obvious ones are rendered useless. For 
example, due to the huge irreducible QCD background from bb dijets, the decay mode 
jjSM ^ jg ciifgcult and one has to complement the search with many other channels, 
particularly H^^^ 77 [13,14]. We have shown in this work that there are scenarios in 
the 2HDM parameter space where alternative decays, like the FCNC modes ^ tc + ic 
and H° ^ tc + tc, can also be useful. For instance, in the case, this situation occurs 
when tan/3 and tana are both large and the CP-odd state is much heavier than the 
CP-even ones. The potential enhancement is then spectacular and it may reach up to 
ten billion times the SM value B{H^^^ ^ tc) ^ 10^^^, thereby bringing the maximum 
value of the FCNC branching ratio B{h^ — tc) to the level of ~ 10^^. As a matter of 
fact, the enhancement would be much larger were it not because we eventually apply the 
severe (indirect) lower bound on the charged Higgs mass from 6^57 [62]. Although 
these decays have maximal ratios below B[h '-/'-/) ~ 10"'^, they should be essentially 
free of QCD background ^. 

While in the MSSM almost the full (m^o, tan/5)-parameter space is covered, with 
better efficiency at high tan /3 though, we should insist that within the general 2HDM the 
tagging strategies are not so well studied and one would like to have further information 
to disentangle the MSSM scenarios from the 2HDM ones. Here again the study of the 
FCNC Higgs decays can play a role. Of course the statistics for the FCNC Higgs decays 
is poor due to the weakness of the couplings and the large masses of the Higgs bosons 

^Misidentification of 6-quarks as c-quarks in tb production might be a source of background to our 
FCNC events. However, to rate the actual impact of that misidentification one would need a dedicated 
simulation of the signal versus background, which is beyond the scope of this work. 
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to be produced. However, in the favorable regions, which are generally characterized 
by large values of tan/5 and of tana, one may collect a few hundred events of the type 
()5.H) - mainly from - in the high luminosity phase of the LHC As we have said, this 
is basically due to the enormous enhancement that may undergo the FCNC decay rates, 
but also because in the same regions of parameter space where the Bs are enhanced, also 
the LHC production rates of the Higgs bosons can be significantly larger (one order of 
magnitude) in the 2HDM II as compared to the SM. 

Interestingly enough, in many cases one can easily distinguish whether the enhanced 
FCNC events ()5.H) stem from the dynamics of a general, unrestricted, 2HDM model, or 
rather from some supersymmetric mechanisms within the MSSM. This is already obvious 
from the fact that the ranges of neutral and charged Higgs boson masses in the 2HDM 
case can be totally incompatible with the corresponding ones in the MSSM. But there 
are many other ways to discriminate these rare events. For instance, in the 2HDM case 
the CP-odd modes ^ tc + ic are completely hopeless whereas in the MSSM they can 
be enhanced see chapter [43,47,48]. Using this information in combination with the 
masses of potentially detected Higgs bosons could be extremely useful to pinpoint the 
supersymmetric or non-supersymmetric nature of them. We may describe a few specific 
strategies. As it was first shown in Ref. [43], the leading SUSY-FCNC effects associated 
to the htc vertices {h = , A^) come from the FCNC gluino interactions which 

are induced by potentially large misalignments of the quark and squark mass matrices 
[93, 144]. These effects are not particularly sensitive to tan (3 and they can be very sizeable 
for both high and moderately low values of this parameter. This sole fact can be another 
distinguishing feature between FCNC events (|5.ip of MSSM or 2HDM origin. If, for 
example, a few of these events were observed and at the same time the best MSSM fits to 
the electroweak precision data would favor moderate values of tan say in the range 10 — 
20, then it is clear that those events could originate in the FCNC gluino interactions but 
in no way within the context of the general 2HDM. In this respect it should be mentioned 
that the FCNC gluino couplings became more restricted from the low-energy meson data 
[104] , and will presumably become further restricted in the near future. The reason being 
that the same couplings are related, via SU (2) gauge invariance and CKM rotation, to 
those affecting the down-like quark sector, which will most likely become constrained by 
the increasingly more precise low-energy meson physics [95,96,104]. In that circumstance 
the only source of FCNC Higgs decays in the MSSM will stem purely from the electroweak 
interactions within the super-CKM basis. Then, in the absence of these SUSY-QCD 
FCNC effects, we could judiciously conclude from the work of Ref. [43] - in which both 
the SUSY-QCD and the SUSY electroweak contributions were computed for the htc 
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vertices - that the FCNC rates in the MSSM should diminish dramatically (two to three 
orders of magnitude). In such case we can imagine the following "provocative" scenario. 
Suppose that the LHC finds a light neutral Higgs boson of mass < 140 GeV (suggestively 
enough, in a mass range near the MSSM upper bound for m/^o!) and subsequently, or about 
simultaneously, a charged Higgs boson and another neutral Higgs boson both with masses 
around 400 GeV or more. At this point one could naively suspect that a MSSM picture 
out these findings is getting somehow confirmed. If, however, later on a few FCNC events 
()5.H1 are reported and potentially ascribed to the previously discovered heavy neutral 
Higgs boson (presumably H^), then the overall situation could not correspond at all 
to the MSSM, while it could be perfectly compatible with the 2HDM II. Alternatively, 
suppose that the FCNC gluino couplings were not yet sufficiently restricted, but (still 
following the remaining hypotheses of the previous example) a third neutral Higgs boson 
(presumably A^) is found, also accompanied with a few FCNC events. Then this situation 
would be incompatible with the 2HDM II, and in actual fact it would put forward strong 
(indirect) evidence of the MSSM!! 

We should also mention that there are other FCNC Higgs decay modes, as for example 
h bs + bs, which could be, in principle, competitive with the top quark modes ()5.H1 . 
In some cases these bottom modes can be highly enhanced in the MSSM case [47,48]. 
Actually, a more complete assessment of the FCNC bottom modes in the MSSM case is 
studied in the chapter IHl - namely one which takes also into account the supersymmet- 
ric contributions to the highly restrictive radiative B-meson decays - shows that they 
are eventually rendered at a similar level of the top modes under study in most of the 
parameter space. 

To summarize, the FCNC decays of the Higgs bosons into top quark final states can 
be a helpful complementary strategy to search for signals of physics beyond the SM in the 
LHC. Our comprehensive numerical analysis shows that the FCNC studies are feasible for 
CP-even Higgs masses up to about 500 GeV. While the statistics of these FCNC decays 
is of course poor, the advantage is that a few tagged and well discriminated events of this 
sort could not be attributed by any means to the SM, and therefore should call for various 
kinds of new physics. In this work we have shown that a general 2HDM II is potentially 
competitive to be ultimately responsible for these FCNC decays, if they are ever found, 
and we have exemplified how to discriminate this possibility from the more restricted one 
associated to the MSSM. 



Chapter 6 



Higgs Boson FCNC Decays into 
Bottom Quarks in the MSSM 

6.1 Introduction 

The most general MSSM includes tree-level FCNCs among the extra predicted particles, 
which induce one-loop FCNC interactions among the SM particles, as discussed in section 

Eini 

Concerning the FCNC interactions of Higgs bosons with third generation quarks, it was 
demonstrated long ago [43] that the leading term corresponds to a single particle insertion 
approximation, which produces a flavor change in the internal squark loop propagator, 
since in this case the chirality change can already take place at the squark-squark-Higgs 
boson interaction vertex. Adding this to the fact that the Higgs bosons (in contrast to 
gauge bosons) have a privileged coupling to third generation quarks, one might expect 
that the FCNC interactions of the type quark- quark- Higgs bosons in the MSSM become 
highly strengthened with respect to the SM prediction. This was already proven in the 
rare decay channels r(t ch) [43] {h being any of the neutral Higgs bosons of the MSSM 
h = if", A^), where the maximum rate of the SUSY-QCD induced branching ratio 
was found to be B(t ch) ^ 10~^, eight orders of magnitude above the SM expectations 
B(t — > cH^^^) ~ 10~^^. Similar enhancement factors have been found in the top-quark- 
Higgs boson interactions in other extensions of the SM, as seen in chapters |31|7| 

From the experience of the previous calculations with the top quark, we expect similar 
enhancements in the FCNC interactions of the MSSM Higgs bosons with the bottom 
quark. Indeed, the purpose of this work is to quantify, in a reliable way, the MSSM 
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Figure 6.1: One-loop SUSY-QCD vertex diagram contributing to the decay h ^ qq' and 
diagrams contributing to mixed h — s self-energy. d{a,i3} represent mass-eigenstate down 
type squarks of any generation. 



expectations on the FCNC Higgs boson decay modes 

h^hs, h^bs {h = h'^ , H\ A^) . (6.1) 

The Feynman diagrams for these decays are depicted in Fig. 16.11 There are other FCNC 
decay modes involving light quarks. However, only these bottom quark channels are 
relevant, as the remaining FCNC decays into light quarks are negligible in the MSSM. 
Moreover, the FCNC decays of Higgs bosons into bottom quarks are specially interesting 
as they can provide an invaluable tool to discriminate among different extended Higgs 
boson scenarios in the difficult LHC range 90 < < 130 GeV [13,14]. 

In this work we present what we believe is the first realistic estimate of the SUSY- 
QCD contributions to the FCNC branching ratios of the MSSM Higgs bosons into bottom 
quark. Specifically, we compute 

for the three Higgs bosons of the MSSM, h = h^,H^,A^, where T{h X) is the - 
consistently computed - total width in each case. The maximization process of the above 
branching ratios in the MSSM parameter space is performed on the basis of a simultaneous 
analysis of the relevant partial decay widths and of the branching ratio of the low-energy 
FCNC process b ^ s'-f, whose value is severely restricted by experiment [52-58]. It turns 
out that the maximum FCNC rates that we find disagree quite significantly with Ref. [47] 
as they do not impose the restriction oi B{b ^ s'-y)} 



^Scc also Rcf. [48] for a combined analysis of flavor- violating and CP- violating MSSM couplings. 
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The structure of the study is as follows. In Section 16.21 we estimate the expected 
branching ratios and describe the structure of Eq. ()6.2|) in the MSSM in more detail; in 
Section lUTSl we present the numerical analysis, and in Section lU^ we deliver our conclusions. 

6.2 Partial widths and branching ratios 

As seen in eqs. (|1.5p and the branching ratios into bottom quark are much larger 

than the Higgs boson FCNC branching ratio into top quark in the SM as seen in chapter 
El However, even in the case p.5|) it is still too small to have a chance for detection in 
the LHC. It is clear that unless new physics comes to play the process H^'^^ ^ bs (and 
of course if-^*^ —*bs) will remain virtually invisible. Nonetheless the result is not 
too far from being potentially detectable, and one might hope that it should not be too 
difficult for the new physics to boost it up to the observable level. 

Consider how to estimate the potentially augmented rates for the MSSM processes 
if only within a similarly crude approximation as above. Because of the strong 
FCNC gluino couplings mentioned in Section 16.11 and the tan /3-enhancement inherent to 
the MSSM Yukawa couplings (see Ref. [43] for details), we may expect several orders 
of magnitude increase of the branching ratios ()6.2j) as compared to the previous SM 
result. A naive approach might however go too far. For instance, one could look at the 
general structure of the couplings and venture an enhancement factor typically of order 
(as/aw)'^ tan^/? which for ^23 ^ 1 and tan/? > 30 could easily rocket the SM 

result some 5 — 6 orders of magnitude higher, bringing perhaps one of the MSSM rates 
()6.2j) to the "scandalous" level of 10% or more. But of course only a more elaborated 
calculation, assisted by a judicious consideration of the various experimental restrictions, 
can provide a reliable result. As we shall see, a thorough analysis generally disproves the 
latter overestimate. 

The detailed computation of the SUSY-QCD one-loop partial decay widths T{h q q') 
in (Q within the MSSM follows closely that of T{t ch) (see Ref. [43]). The rather 
cumbersome analytical expressions will not be listed here as they are an straightforward 
adaptation of those presented in the aforementioned references. However, there are a 
few subtleties that need to be pointed out. One of them is related to the calculation 
of the total widths T{h X) for the three Higgs bosons h = in the MSSM. 

As long as r(/i ^ qq') in the numerator of Eq. ()6.2|) is computed at leading order, the 
denominator has to be computed also at leading order, otherwise an artificial enhancement 
of B{h qq') can be generated. For example, including the next-to-leading (NLO) order 
QCD corrections to T{h — > bb) reduces the decay width by a significant amount [145-149]. 
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Then, to be consistent, the NLO (two-loop) contributions to T{h qq') should also be 
included. Similarly, the one- loop SUSY-QCD corrections to r{h ^ bb) can be very large 
and negative [31], which would enhance B{h — qq'). At the same time these corrections 
also contribute to r{h qq'), such that contributions to the numerator and denominator 
of Eq. (|6.2j) compensate (at least partially) each other. Therefore the same order of 
perturbation theory must be used in both partial decay widths entering the observable 
B{h qq') to obtain a consistent result. By the same token, using running masses in the 
numerator of ()().2|1 is mandatory, if they are used in the denominator. Last, but not least, 
consistency with the experimental bounds on related observables should also be taken 
into account. In this respect an essential role is played by the constraints on the FCNC 
couplings from the measured value of B{b — > 57). They must be included in this kind of 
analysis, if we aim at a realistic estimate of the maximal rates expected for the FCNC 
processes in the MSSM. In our calculation we have used the full one-loop MSSM 
contributions to 5(6 — > 57) as given in [150] ^. 

Let us now summarize the conditions under which we have performed the computation 
and the approximations and assumptions made in the present analysis: 

• We include the full one-loop SUSY-QCD contributions to the partial decay widths 
r{h qq') in (Q. 

• We assume that FCNC mixing terms appear only in the LH-chiral sector of the 
squark mixing matrix. This is the most natural assumption, and, moreover, it was 
proven in Ref. [43] that the presence of FCNC terms in the RH-chiral sector enhances 
the partial widths by a factor two at most - not an order of magnitude. 

• The Higgs bosons total decay widths T{h X) are computed at leading order, 
including all the relevant channels: r{h ff, ZZ, W^W~ , gg). The off-shell decays 
r{h ZZ*, W^W^*) have also been included. The one-loop decay rate r{h — > gg) 
has been taken from [151] and the off-shell decay partial widths have been computed 
explicitly and found perfect agreement with the old literature on the subject [152]. 
We have verified that some of the aforementioned higher order decays are essential 
to consistently compute the total decay width of r{h^ X) in certain regions of the 
parameter space where the maximization procedure probes domains in which some 
(usually leading) two-body processes become greatly diminished. We have checked 
that our implementation of the various Higgs boson decay rates is consistent with 

^Ref. [150] contains a partial two-loop computation of B{b 57), but only the one-loop contributions 
have been used for the present work. 
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the results of HDECAY [153]. However, care must be exercised if using the full- 
fledged result from HDECAY. For example, it would be inconsistent, and numerically 
significant, to compute the total widths T{h X) with this program and at the 
same time to compute the SUSY-QCD one-loop partial widths V{h — > qq') without 
including the leading conventional QCD effects through e.g. the running quark 
masses. 

• The Higgs sector parameters (masses and CP-even mixing angle a) have been 
treated using the leading mt and mbtan/5 approximation to the one-loop result [85, 
154-156]. For comparison, we also perform the analysis using the tree-level approx- 
imation. 

• We include the constraints on the MSSM parameter space from B{h — S7). We 
adopt B{b —>■ 57) = (2.1 — 4.5) x 10~^ as the experimentally allowed range within 
three standard deviations [58]. Only the SUSY-QCD contributions induced from 
tree-level FCNCs are considered in the present work. 

Running quark masses {mg{Q)) and strong coupling constants (^^(Q)) are used through- 
out. More details are given below, as necessary. 



6.3 Full one-loop SUSY-QCD calculation: Numerical 
analysis 

Given the setup described in Section 16.21 we have performed a systematic scan of the 
MSSM parameter space with the following restrictions: 

623 < 10-°-°9 ~ 0.81 

Ab = -1500- ■■ 1500 GeV , , 

6.3 

fi = -1000- ■ -1000 GeV 
rrig = 150 ■■■ 1000 GeV 

and the following fixed parameters: 

tan/3 = 50 

= = ^in =^g = (6-4) 

At = -300 GeV . 

Here rrij^^ ^ are the left-chiral and right-chiral bottom-squark soft-SUSY-breaking mass 
parameters, and rriq is a common mass for the strange- and down-squark left- and right- 
chiral soft SUSY-breaking mass parameters. Following the same notation as in [43], the 
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(a) (b) 

Figure 6.2: Maximum SUSY-QCD contributions to B{h — > qq'), Eq. (j6.2j) . as a function 
of a) m^o and b) ^23 for "^ao = 200 GeV. 



parameter represents the mixing between the second and third generation squarks. 
Let us recall its definition: 

2 

mi _ 

523 = , (6.5) 

m- . being the non-diagonal term in the squark mass matrix squared mixing the second 
and third generation left-chiral squarks. The parameter ^23 is a fundamental parameter in 
our analysis as it determines the strength of the tree-level FCNC interactions induced by 
the supersymmetric strong interactions, which are then transferred to the loop diagrams 
of the Higgs boson FCNC decays 

The result of the scan is depicted in Fig. 16.21 To be specific: Fig. 16.2b shows the 
maximum value B'^^^iji — > qq') of the FCNC decay rate ()6.2|) under study as a function 
of m^o; Fig. 16.2b displays B^^^ih qq') as a function of the mixing parameter ^23 
for m^o = 200 GeV. Looking at Fig. 16.21 three facts strike the eye immediately : i) the 
maximum is huge (13%!) for a FCNC rate, actually it is as big as initially guessed from 
the rough estimates made in Section ii) very large values of ^23 are allowed; iii) the 
maximum rate is independent of the pseudo-scalar Higgs boson mass m^o. We will now 
analyze facts ii) and iii) in turn, and will establish their incidence on fact i). For further 
reference, in Table lUTTl we show the numerical values of B"^^^{h — > qq') together with the 
parameters which maximize the rates for m^o = 200 GeV. 

One would expect that a large value of ^23 should induce a large gluino contribution 
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Particle 








B{h an') 


3.3 X 10~2 


1.3 X 10^1 


3.3 X 10-2 


Y[h X) 


11.0 GeV 


1.6 X 10"^ GeV 


11.3 GeV 




10-0.09 


10-0.1 


10-0.09 


nig 


975 GeV 


975 GeV 


975 GeV 


Ab 


1500 GeV 


730 GeV 


1290 GeV 


fi 


980 GeV 


1000 GeV 


980 GeV 


B{b 57) 


4.42 X 10"^ 


4.23 X 10-^ 


4.50 X 10"^ 



Table 6.1: Maximum values of B{h qq') and corresponding SUSY parameters for 
niAo = 200 GeV. 



Particle 






A^ 


B{h qq') 


9.1 X 10-"^ 


3.1 X 10-3 


9.1 X 10-4 


r{h ^ X) 


11.2 GeV 


1.4 X 10-3 GeV 


11.3 GeV 


^23 


10-0.43 


10-0.8 


10-0.43 


nig 


1000 GeV 


975 GeV 


1000 GeV 


Ab 


-1500 GeV 


-1500 GeV 


-1500 GeV 




-460 GeV 


-1000 GeV 


-460 GeV 


B{b S7) 


4.49 X 10-^ 


4.48 X 10-^ 


4.49 X 10-4 



Table 6.2: Maximum values of B{h qq') and corresponding SUSY parameters for 
niAo = 200 GeV excluding the window region. 



to B{b 57). In fact it does! However our automatic scanning process picks up the 
corners of parameter space where the gluino contribution alone is much larger than the 
SM contribution, but opposite in sign, such that both contributions destroy themselves 
partially leaving a result in accordance with the experimental constraints. We exam- 
ine this behaviour in Fig. 16.31 where we show the values of B{h qq') together with 
B[b 57) as a function of ^23 for the parameters which maximize the FCNC rate of the 
lightest CP-even state /i" in Table 16.11 We see that, for small values of ^23; the gluino 
contribution to B{b 57) is small, and the total B{b 57) prediction is close to the 
SM expectation. In contrast, as ^23 steadily grows, B{b ^ 57) decreases fast (mean- 
ing a dramatic cancellation between the two contributions) until reaching a point where 
B{b — > S7) = 0. From there on it starts to grow with a large slope, and in its race 
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Figure 6.3: B{h — > qq') and B[h — > 37) as a function of ^23 for the parameters that 
maximize B{h^ qq') in Table lUTTl The shaded region is excluded experimentally. 



eventually crosses the allowed B{b ^ 57) region. The crossing is very fast, and so rather 
ephemeral in the ^23 variable, and it leads to the appearance of a narrow allowed window 
at large ^23 values, see Fig. Ifj.Hk . We would regard the choice of this window as a fine- 
tuning of parameters, hence unnatural ^. For this reason we reexamine the B{h — > qq') 
ratio by performing a new scan of the MSSM parameter space in which we exclude the 
fine-tuned (or window) region. The result for m^o = 200 GeV can be seen in Table 16.21 
and Fig. 16.41 This time we see that the maximum values of B{h qq') are obtained for 
much lower values of ^23, and the maximum rates have decreased more than one order of 
magnitude with respect to Table lU^Tl reaching the level of few per mil. These FCNC rates 
can still be regarded as fantastically large. Had we included the SUSY-EW contributions 
to B{b 57), further cancellations might have occurred between the SUSY-EW and the 
SUSY-QCD amplitudes. Even more: since each contribution depends on a separate set 
of parameters, one would be able to find a set of parameters in the SUSY-EW sector 
which creates an amplitude that compensates the SUSY-QCD contributions for almost 
any point of the SUSY-QCD parameter space. But of course this would be only at the 
price of performing some fine tuning, which is not the approach we want to follow here. 

On the other hand further contributions to 5(6 57) might exist. In the most general 
MSSM, flavor- changing interactions for the right-chiral squarks (^23^7?), and mixing left- 
and right-chiral squarks (523Lij) can be introduced. The latter can produce significant 

■^At the time of writing this Thesis this fine-tuning is excluded experimentaUy [102]. 



6.3 Full one-loop SUSY-QCD calculation: Numerical analysis 



95 



0.004 
0.0035 

0.003 
0.0025 

0.002 
0.0015 

0.001 
0.0005 




B(H->qq') 
B(h°->qq') 
B(A°->qq') 



(a) 




(b) 



Figure 6.4: Maximum value of the SUSY-QCD contributions to B{h qq') as a function 
of a) m^o and b) for f^A° = 200 GeV, for the scenario excluding the window regions. 



contributions to B{b — > 57), changing the allowed parameter space. The introduction 
of S23LR can produce two possible outcomes: First, in certain regions of the parameter 
space, the contributions of S23LR and are of the same sign, enhancing each other. In 
this situation, the maximum allowed value of ^23 is obtained for 623LR = 0. Second, in 
other regions of the parameter space the two contributions would compensate each other, 
producing an overall value oi B{b ^ 37) in accordance with experimental constraints, 
even though each contribution would be much larger. Again, we would regard these 
compensations as unnatural, and would discard that region of the parameter space. In 
the following we will require that the SUSY-QCD contributions induced by ^23 do not 
compensate the SM ones to give an acceptable value oi B{h S7); this is equivalent to 
the condition that the SUSY-QCD amplitude represents a small contribution to the total 
B(h ^ 57) value, and is therefore independent of the inclusion of the other contributions 
(SUSY-EW, 623lr).^ 

We turn now our view to the second fact, namely the independence of the maximum 
rates with respect to m^o. We will show that it also plays a central role as to the 
enhancement of B{h qq')- Actually, a good hint is given by the small values of the 
lightest Higgs boson decay width in Tables O and EH r(/i° ^ X) ~ 2 x lO-^GeV. 
The maximization process of B{h^ —* qq') does not only find the parameters for which 

^The analysis of Ref. [48] follows the opposite approach, that is: to find the fine-tuning conditions 
imposed by low energy data that allow for the largest possible value of the FCNC parameters. 



96 



Higgs Boson FCNC Decays into Bottom Quarks in the MSSM 









- 































)/ 

/\ 




•3 Or(h"->X) [GeV] \ 








B(H°/A°-> q q') 








- B(h°->qq') 


1 1 







-1500 -1000 



-500 500 

H [GeV] 



1000 1500 




-1500 -1000 



-500 

H [GeV] 



500 1000 1500 



(a) (b) 

Figure 6.5: B{h — * qq') and r{h^ X) (in GeV) as a function of fi for a) one-loop 
a angle; b) tree-level a angle, and for the parameters that maximize B{h^ qq') in 
Table lUTTl The if" and curves coincide. The B{b s'-f) constraint is not shown. 



r(/i° — * qq') is maximum, but also the parameters for which r{h^ X) is minimum. 
Specifically, since r{h^ — > bb) is the dominant decay decay channel of for large tan/3, 
the maximum of B{h^ qq') is produced in the parameter range of the so-called small 
Ueff scenario [157], that is, a parameter range where the radiative corrections make the 
CP-even Higgs boson mixing angle a vanish (or very small) , such that the leading partial 
decay width T{h^ bb) is strongly suppressed. The consequences of this scenario have 
been extensively studied in Ref. [158]. As advertised in Section [6.21 the possibility that 
the maximization process explores these regions of the parameter space is the reason why 
the leading higher order decay channels, and also the leading three-body decay modes 
have to be taken into account in the computation of the total width. 

In Fig. 16.51 we plot the value of the various branching ratios B{h qq') and of 
the total width of the lightest CP-even Higgs boson, T{h^ — > X), as a function of the 
higgsino mass parameter /i, the rest of the parameters being those of the third column of 
Table IHTH i.e. the ones that maximize the branching ratio B{h^ — > qq'). Fig. 16.5b shows 
that T{h^ X) has a deep minimum in the range of yU corresponding to the maximum of 
B{h^ qq'), which reaches the level of a few percent. If, instead of using the radiatively 
corrected a value we use the tree-level expression, we obtain the result shown in Fig. 16.5b . 
Here the total decay width of the Higgs boson is independent of and B{h qq') does 
not show any peak. Actually in this case the branching ratio for becomes smaller than 
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Particle 








B{h an') 


9.0 X 10""^ 


1.3 X 10"^ 


9.0 X 10"^ 


Y[h X) 


11.3 GeV 


5.4 X 10-3 GeV 


11.3 GeV 




10-0.43 


10-0.28 


10-0.43 


nig 


1000 GeV 


1000 GeV 


1000 GeV 


Ab 


-1500 GeV 


-1500 GeV 


-1500 GeV 


fi 


-460 GeV 


-310 GeV 


-460 GeV 


B{b 57) 


4.49 X 10"^ 


4.50 X 10"^ 


4.49 X 10"^ 



Table 6.3: Maximum values of B{h qq') and corresponding SUSY parameters for 
m^o = 200 GeV, using the tree-level expressions for the Higgs sector, and excluding the 
window region. 



that of and A^ for all /i. The maximization procedure in Fig. 16. 21 selects for each value 
of m^o the MSSM parameters corresponding to the small «£// scenario for that specific 
value of niAo. Of course, this discussion regarding the channels for large values of m^o 
has a correspondence with the channel for low values^ of m^o. 

As indicated in Section 16. 2[ we have used a one-loop approximation for the Higgs 
sector [85,154-156], instead of the more sophisticated complete two-loop result present 
in the literature [114,115]. However, we should stress that the exact MSSM parameters 
at which the small Oe// scenario is realized are not important for the sake of the present 
analysis. All that matters is that some portion of the parameter space exists, for which 
r(/i° bb) is strongly suppressed, but r{h^ qq') is not. 

To compare the maximum value of B{h^ qq') obtained with and without the 
small a^ff scenario, we have performed the maximization procedure using the tree-level 
expressions for the Higgs sector parameters. The result is shown in Table and Fig. 16.61 
In this case B{h^ — ^ qq') is reduced by a sizeable factor of > 20 with respect to TableEIS 
whereby the rate descends about an order of magnitude below that of the 
channels which remain basically unchanged. Notice also that r{h^ X) is larger than 
in previous tables. In spite of the reduction, achieving a FCNC ratio B{h^ ~^ QQ') ^ 
1.3 X 10"^ is a remarkable result, three orders of magnitude larger than the maximum SM 
rate (jl.5|) . and only one order of magnitude below the rare decay B{h^ — > 77) ~ 10"^. Also 

^Largc or low values here means m^o > m^^^ or m^o < m™^'^, i.e. above or below the maximum 
possible value for the mass of the lightest Higgs boson h^, respectively. 
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Figure 6.6: Maximum value of the SUSY-QCD contributions to B{h qq') as a function 
of a) m^o and b) ^23) for f^A" = 200 GeV and for the scenario excluding the window 
region and using the tree-level expressions for the Higgs sector parameters. 



worth noticing in Fig. 16.6b (and Fig. 16.4b ) is the fact that B^^^{h qq') is essentially flat 
in (523 in the upper range down to ^23 ~ 10~*^'^ ^ 0.16. The reason lies in the correlation 
between B{h qq') and B{h — ^ 57). In order to comply with the (non-fine-tuned) value 
of 5(6 57) for large 523, the absolute value of the 11 parameter must be small. When ^23 
decreases, can grow to larger values, leaving the overall maximum rates B^^^{h —* qq') 
effectively unchanged (see Eq. ()6.7|) below). 

The maximization process selects a squark mass scale in the vicinity of the maximum 
values used in the scanning procedure. We should point out, however, that the same order 
of magnitude for B{h — > qq') could be obtained with a much lower squark mass scale. 
In this case the lighter squark masses induce a much larger B{b ^ S7) value, and ^23 is 
much more constrained. For example, if we perform a scan in the parameter space ()6.3|) . 
but fixing the squark mass scale to be < 500 GeV, we obtain the following values for 
the maximal branching ratios for m^o = 200 GeV: 

^max^^^o ^qq') = 1.4 X 10"^ , B'^'^'^iHyA^ ^ qq') = 9.2 x 10"^ , (6.6) 

with ^23 ~ 10"°'^, /i ~ —110 GeV, and we have limited ourselves to the scenario avoiding 
the window regions and using the tree-level expression for the Higgs sector parameters. 
These numbers have to be compared with Table 

The reason behind this scale independence admits an explanation in terms of an ef- 
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Figure 6.7: One-loop SUSY-QCD vertex diagram contributing to the decay h ^ qq' and 
diagrams contributing to mixed h — s self-energy in the mass approximation 



fective Lagrangian approach, in which one can estimate the leading effective coupling to 
behave approximately as 16.71 [48] : 



0^ 
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Aside from ensuring (at least) a partial SUSY scale independence of the leading terms, this 
expression also shows that B{h — > qq') has a weak dependence on the soft-SUSY-breaking 
trilinear coupling A^. The observed situation is similar to the flavor-conserving hhh inter- 
actions, where the cancellation of the Ai, terms at leading order has been proven [41]. It 
also shows that the leading non- decoupling SUSY contributions to T{h^ q q') eventually 
fade out as the decoupling limit of the Higgs sector is approached: cos(/3 — a) ^ 0. We 
have found (using the tree-level expression for a) that the non-leading (SUSY-decoupling) 
contributions to r(/i° qq') dominate for m^o > 450 GeV, inducing a value r™'^^(/i° — > 
qq') ~ 1.2 X 10-^ with 623 ~ lO'S /i ~ 1000 GeV. 

We further investigate the role of the scale of SUSY masses, and the fine-tuning be- 
haviour in Figs. 16.81 and |6.9I In these figures we give up the equality rrig = rrig (j6.4j) . the 
squark masses are fixed at the values stated in Tables 16.21 and 16.11 respectively. Fig. 16.81 
shows the values of B{h — ^ qq') for the three Higgs decays and of B{b 57) as a func- 
tion of the gluino mass for the parameters that maximize B{h^ qq') when the window 
regions are excluded (third column of Table . Here we see that, while the gluino con- 
tribution to B{b — > 57) decouples fast as a function of nig, its contribution to B{h qq') 
is fairly sustained. Indeed, between rrig = 1 TeV and rrig = 5TeV B{h'^ qq') decreases 
only by a factor ~ 1/4, while the gluino contribution to B{b 57) becomes negligible 
at rrig = 5 TeV and we recover the SM prediction. As a consequence, the maximum rates 
B[h qq') that we have found are robust, in the sense that further theoretical refine- 
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Figure 6.8: B{h — > qq') and B{h 57) as a function of rrtg for the parameters that 
maximize B{h^ — > hs) excluding the window region (see third column of Table IHT^ . The 
shaded region is excluded experimentally. 
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Figure 6.9: As in Fig. 16. 81 but including the window region. The remaining parameters 
are fixed as in the third column of Table 16.11 
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Table 6.4: Maximum values of B{h — > qq') and corresponding r{h — > X) for the different 
scenarios studied in this work. 



ments and experimental results that change the allowed range oi B{b ^ s'-f) can easily 
be compensated for by a slight increase of the gluino mass {mg), which would leave the 
prediction for B{h qq') essentially unchanged. We note in Fig. l6.9l the correspond- 
ing behaviour of B{h^ qq') and B{b s'-f) in the presence of fine-tuning, i.e. as 
in Table lUm In contrast to the previous case, here we observe the presence of two tiny 
windows in the regions nig = 25 — 75 GeV and rrig = 950 — 1125 GeV. In the middle 
region rrig = 75 — 950 GeV, B{b s'-f) is one order of magnitude larger than the allowed 
experimental range, and in the region above rrig = 1125 GeV it only enters the allowed 
region for rrig > 4500 GeV. In this region B{h^ — > qq') is still large, but at the price of 
having a gluino five times heavier than the rest of the SUSY spectrum. This is another 
manifestation of the large fine-tuning that governs this region of the parameter space. 

Up to this point we have used the high tan/5 value quoted in Eq. (j6.4j) . But we have 
also looked at the impact of varying tan/3 on B™^^{h qq')- Since the latest LEP 
data restricts tan/? > 2.5, we have used a moderate value of tan/3 = 5. Note that, at low 
tan /3, the small scenario does not arise. As a consequence similar results are obtained 
using either the tree-level or one-loop expressions for the Higgs sector parameters. We 
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find that the three branching ratios B^^^{h qq') at tan/3 = 5 stay in the same order 
of magnitude as in the scenarios with tan /5 = 50 (default case) with the tree-level Higgs 
sector and no-window (Cf. Table lU^ . 

6.4 Remarks and conclusions 

The main numbers of our analysis are put in a nutshell in Table 16.4^ where we show the 
results presented previously, together with some other scenarios and the low tan/5 case. 
The computed maximum values of B{h qq') must not be taken as exact numbers 
in practice, but order of magnitude results. The implications that can be derived from 
Table 16.41 can be synthesized as follows: 

1. The SUSY-QCD contributions can enhance the maximum expectation for the FCNC 
decay rates B{h — qq') enormously. This is seen by comparing the results of 
Table with the maximum value of B{H^^'^ — > hs) considered in Eq. (jl.5j) . The 
optimized MSSM branching ratios are at the very least 3 orders of magnitude bigger 
than the SM result. 

2. If no special circumstances apply, that is, if no fine-tuning occurs between the param- 
eters contributing to B{h — S7) in the MSSM, and if T{h^ — > hh) is not suppressed, 
the maximum rates are B^'^'^Qi^ qq') ~ 1.3 x 10""', B'"''''^{H^ / ~^ Qq') — 
9 X 10~^. This corresponds to the tree- Higgs / no-window scenario in Table IfT^ 

3. If, however, r{h^ hh) is suppressed by the radiative corrections to the CP-even 
mixing angle a, then B(hP qq') can be an order of magnitude larger: B^^^{h^ — > 
qq') ~ 3 X 10"'^. This corresponds to the small a^ff scenario, and is indicated by 
small-ag/// no-window in Table The FCNC branching ratio that we find for 

in this case should be considered as the largest possible one within the conditions 
of naturalness (no fine-tuning). 

4. On the other hand, if fine-tuning between the gluino and the SM contributions to 
B{h ^ 57) is allowed, but the small-ae// scenario is not realized, then B^^^{h^ 
qq') grows one order of magnitude up to B"^^^{h^ ^ QQ') ~ 4 x 10^'^, whereas 
^max^^o^^o _^ gg/^ ^ 2 ^ 10~^. This corresponds to the case labelled tree- 
Higgs/ window in Table 

5. When both special conditions take place simultaneously, viz. fine-tuning in 5(6 ^ 
57) (triggered by a very special choice of the ^23 parameter in a narrow window 
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range) and small scenario (independent of assumptions on ^23), we reach an 
over-optimistic situation where B™^^{h? qq') could reach the ~ 10% level. This 
is the case referred to as small-agj-j/ window in Table 

6. If tan/5 is low/moderate, then B'^^^iji qq') lie in the lower range ~ 10~^, which 
can grow an order or magnitude for B^^^^{H^ / qq') in fine-tuned scenarios 
(last three rows in Table . 

Although the large FCNC rates mentioned in points |^ and El above seem to offer 
a rather tempting perspective, we will not elaborate on them any further since in our 
opinion the fine-tuning requirement inherent in them is too contrived. On the other 
hand, points |2l and El offer a moderate, but certainly much more realistic scenario, which 
in no way frustrates our hopes to potentially detect the FCNC Higgs boson decays (jfi.lj) . 
Indeed, in the case described in point El B{h — > qq') can be at most of order 10~^. But 
this is still a fairly respectable FCNC branching ratio (comparable to that of 6 — 57) and 
it may lead to a large number of events at a high luminosity collider as seen in chapter [71 
for the MSSM and the chapter jH for the general 2HDM. Moreover, if T{h X) becomes 
suppressed (e.g. by realizing the small ag// scenario, point then B{h^ — > qq') can be 
enhanced by an additional order of magnitude. 

Our analysis correlates the values of B{h qq') with that of B{b 57), taking into 
account only the SUSY-QCD contributions due to flavor mixing parameters among the 
left-chiral squarks. The presence of several other competing contributions to B{b ^ 57) 
alters the borders of the allowed parameter space: 

• For the fine-tuned scenarios, the presence and position of the allowed window regions 
in the parameter space depends significantly on all the contributions, and therefore 
also does the maximum value of B{h qq')- Outside the window regions, the 
computed value oi B[h ^ S7) can only be made consistent with the experimental 
range, by means of a large splitting between the squark and gluino masses. 

• For the non- fine-tuned scenarios, the inclusion of further contributions to B{b ^ S7) 
also alters the allowed parameter space, but the condition of non-fine-tuning ensures 
precisely that the change in the allowed range of is smooth, and the corresponding 
change in B^^{h — > qq') is not dramatic. 

Of course, the question immediately arises on what will happen if the data from present 
S-meson factories further constrains the ^23 parameter. In that case, we should take into 
account the (charged-current induced) SUSY-EW contributions to B{h qq'), see [49]). 
However, we can advance that the SUSY-EW effects on B^^^{h — * q q') that we find are in 
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the ballpark of ^ g g') ~ 3 x IQ-^ and fi™^^(i/7A° ^qq')^lx IQ-^ for a non- 

fine-tuned scenario, while ^ g g') ~ 2 x 10"^ and 5'^^^(i/7A° ^ g g') ~ 8 x 10"^ 

for a fine-tuned scenario. From the analysis of the chapter El we expect that even with 
these impoverished MSSM rates the number of FCNC events of that sort should be non- 
negligible at the LHC. 

Even though we have detected the existence of corners of the MSSM parameter space 
where a Higgs boson FCNC branching ratio can barely reach the 10% level (cf. the narrow 
windows in Fig. Ki.HI and Fig. l(i.9j) . we insist once more that they should be considered 
rather unlikely as they are associated to fine tuning of the parameters and excluded 
experimentally [102]. Moreover, in contrast to Ref [47], we find that it is the lightest 
CP-even state, the one that could have the largest FCNC branching ratio. 

To conclude, we have presented a first realistic estimate of the branching ratios of 
the Higgs boson FCNC decays (jHIH) within the MSSM, assuming that the SUSY-QCD 
corrections can be as large as permitted by the experimental constraints on B{h ^ 57). 
We have carried out a systematic and self-consistent maximization of the branching ratios 
()(i.2|) taking into account this crucial experimental constraint. At the end of the day the 
results that we obtain, especially for the lightest CP-even Higgs boson of the MSSM, 
are fairly large: gg') ~ 10"^ - lO'^. These MSSM rates turn out to be 

between three to four orders of magnitude larger than the maximum SM rate ()1.5|) . but 
not five or six orders as naive expectations indicated. Whether this branching ratio is 
measurable at the LHC [13,14] or at a high energy e^e~ Linear Collider [15] can only be 
established by means of specific experimental analyses. However, on the basis of related 
studies in the general 2HDM (chapter EI) and from the MSSM (chapter [7j), we can foresee 
that an important number of FCNC events ()6.H) can be potentially collected at the LHC. 
They could play a complementary, if not decisive, role in the identification of low-energy 
Supersymmetry. In this work we have dealt only with the maximum rates induced by the 
SUSY-QCD sector of the model. 
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Production and FCNC decay of 
MSSM Higgs bosons into heavy 
quarks in the LHC 



7.1 Introduction 



The rareness of the FCNC Higgs boson decay modes, as seen in the previous chapter El 
are an ideal laboratory to look for non-standard interactions superimposed onto the SM 
ones. Similar considerations apply to the FCNC processes associated to the Hbs vertex, 
but in this case it is more difficult to pin down the phenomenological signatures. Some 
work along these lines has already been done, both in the MSSM (see chapter [HI and 
Refs. [43,45-50]) and in the general two-Higgs-doublet model (2HDM) (see chapters |3] 
and El and Refs. [51,64]), and also in other extensions of the SM- see [25] for a review. 
Up to now, the main effort has been concentrated in computing the FCNC decay modes at 
one-loop within the new physics, and also in getting a realistic estimate of the maximum 
branching ratios expected. It is not enough to compute the FCNC branching ratios in, say 
the MSSM, and then evaluate them in some favorable region of the parameter space, for 
one has to preserve at the same time the stringent bounds on other observables in which 
the same physics can be applied, like the aforementioned low-energy 6 — > 57 decay. This 
kind of correlated study was done very carefully in chapter El for the specific Higgs boson 
FCNC decays into bottom quarks within the MSSM, h bs {h = h°,H°,A°). In this 
work we extend the latter work by computing also the top quark Higgs boson FCNC decay 
modes of the heavy MSSM Higgs bosons, h tc [h = H^, A^) under the same restrictions 
(recall that cannot participate in this decay because m/jO < in the MSSM [11]). 
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Furthermore, in this work we carry out an additional step absolutely necessary to make 
contact with experiment, namely we combine the FCNC decay branching ratios of the 
MSSM Higgs bosons (into both top and bottom quarks) with their MSSM production 
cross-sections in order to estimate the maximum number of FCNC events expected at 
LHC energies and luminosities. Only in this way one can assess in a practical way the 
probability of detecting such processes at the LHC. This computation was done for the 
general 2HDM in chapter^! but to the best of our knowledge the corresponding calculation 
in the MSSM case is not available. In this work we perform this calculation and compare 
the FCNC results obtained for the MSSM and 2HDM scenarios. 

The relevant observable quantity on which we shall focus hereafter is the cross-section 
for the production of electrically neutral pairs of heavy quarks of different flavors at the 
LHC, whose origin stems from the FCNC decays of the neutral Higgs bosons of the MSSM, 
h = h^, H^, [6,11]. Thus, we aim at the quantity 

T{h^qq') 



a(j)p h ^ qq) = a{pp — > hX)B{h ^ qq) = a{pp — hX)' 



,^S{h^qq' + qq') , ,^ 
= a{pp hX) ^ ^ ^ ^ {qq=bsoTtc). (7.1) 

Here r{h X) is the - consistently computed - total width in each case. In order to 
asses the possibility to measure these processes at the LHC, we have performed a scan 
of the MSSM parameter space to find the maximum possible value of the production 
rates ()7.1|) under study. The computation of the combined production rate is necessary, 
since the correlations among the different factors are important. For example, in chap- 
teriniit was shown that the maximum branching ratio for the lightest MSSM Higgs boson, 
B{hP ^ bs), is obtained in the regions of the parameter space where the coupling bb 
is strongly suppressed by quantum effects. On the other hand, the associated production 
a{pp — * h^bb) is one of the leading processes for the production of the lightest MSSM Higgs 
boson. It is clear then that, in the regions where B{h^ bs) is largest, a{pp — > h^bb) will 
be suppressed. Therefore, the maximum FCNC production rate at the LHC can only be 
obtained by the combined analysis of the two relevant factors in fj7.1|) (viz. the branching 
ratio and the Higgs boson production cross-section). We will see that the effects from 
each factor are different in different regions of the parameter space. Moreover, the real- 
istic production FCNC rates (j7.ip in the MSSM parameter space can be obtained only 
by including the restrictions imposed by the simultaneous analysis of the branching ratio 
of the low-energy process b —* s'y, whose range of values is severely limited by experi- 
ment [52-58]. As in chapter IHl in this work we limit ourselves to supersymmetric FCNC 
interactions mediated by the strongly interacting sector of the MSSM, i.e. the SUSY-QCD 
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flavor-violating interactions induced by the gluinos. The corresponding analysis for the 
electroweak supersymmetric FCNC effects requires a lengthy separate presentation, and 
will be reported elsewhere [159]. 

The work is organized as follows. In Section 17^ we describe the general setting for our 
numerical analysis. In Section 17.31 we present the LHC production rates of Higgs boson 
decaying into bottom quarks through supersymmetric FCNC interactions. In Section [73] 
we present the corresponding FCNC rates for the top quark channel. Finally, in Section l731 
we compare the MSSM results with the 2HDM results, and deliver our conclusions. 

7.2 General setting for the numerical analysis 

We have performed the calculations with the help of the numeric programs HIGLU, PPHTT [136, 
139,151,160] and LoopTools [129,161,162]. The calculation must obviously be finite with- 
out renormalization, and indeed the cancellation of UV divergences using either dimen- 
sional regularization or dimensional reduction - the two methods giving the same results 
here - in the total amplitudes was verified explicitly. In the following we will detail the 
approximations used in our computation: 

• We include the full one-loop SUSY-QCD contributions to the FCNC partial decay 
widths r(/i — qq') in the observable ()7.1|) . 

• We assume that FCNC mixing terms appear only in the LH-chiral sector of the 
6x6 squark mixing matrix. Therefore, this matrix has only non-diagonal blocks in 
the LH-LH sector. This is the most natural assumption from the theoretical point 
of view [93], and, moreover, it was proven in Ref. [43] that the presence of FCNC 
terms in the RH-chiral sector would enhance the partial widths by a factor two at 
most - not an order of magnitude. 

• The Higgs sector parameters (masses and CP-even mixing angle a) have been 
treated using the leading rrit and mfetan/5 approximation to the one-loop result [85, 
154-156]. For comparison, we also perform the analysis using the tree-level approx- 
imation. 

• The Higgs bosons total decay widths r{h —>■ X) are computed at leading order, 
including all the relevant channels: r{h ff, ZZ, W~^W~ , gg). The off-shell decays 
r{h ZZ*,W^W^*) have also been included. This is necessary to consistently 
compute the total decay width of r{h^ X) in regions of the parameter space 
where the maximization of the cross-section (j7.ip is obtained at the expense of 
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greatly diminishing the partial decay widths of the two-body process — > bb (due 
to dramatic quantum effects that may reduce the CP-even mixing angle a to small 
values [157]). The one-loop decay rate r{h — > gg) has been taken from [151] and 
the off-shell decay partial widths have been recomputed explicitly and found perfect 
agreement with the old literature on the subject [152]. 

• The MSSM Higgs boson production cross-sections at the LHC have been computed 
using the programs HIGLU 2.101 and PPHTT 1.1 [136,139,151,160]. These pro- 
grams include the following channels: gluon-gluon fusion (T{pp{gg) — * h), associated 
production with top-quarks a{pp htt) and associated production with bottom- 
quarks a{pp hbb). In order to have a consistent description, we have used the 
leading order approximation for all channels. The corresponding Feynamn diagrams 
are depicted in 15.21 The QCD renormalization scale is set to the default values for 
each program, namely /iq = rrih for HIGLU and fiQ = {m^ + 2Mq)/2 for PPHTT. We 
have used the set of CTEQ4L Parton Distribution Functions [163] . 

Running quark masses {mq{Q)) and strong coupling constants (a^,(Q)) are used through- 
out, with the renormalization scale set to the decaying Higgs boson mass in the decay 
processes. More details are given below, as necessary. 

Using this setup, we have performed a maximization of the FCNC cross-section, 
Eq. ()7.1|1 . in the MSSM parameter space with the following restrictions on the parameters: 
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Here nig^ are the LH-chiral and RH-chiral squark soft-SUSY-breaking mass parameters, 
rrig^ H) common for the three generations; Mg. are the physical masses of the squarks, and 
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rrih is the mass of the decaying Higgs boson h = h^, H^, A^. These masses are fixed at the 
tree-level by the values of (tan j3, m^o) and the SM gauge boson masses and couplings [6]. 
Due to the structure of the Yukawa couplings in the MSSM, the value of tan/? is fixed at 
a high (small) value for the bottom (top) quark channel as indicated. The parameter m^o 
(the mass of the CP-odd Higgs boson) is assumed to vary in the range indicated in (|7.2j) . 
At one loop these masses receive corrections from the various SUSY fields, and therefore 
depend on the values of the remaining parameters in Eq. ()7.2|) . The characteristic SUSY 
mass scale Msusy defines the typical mass of the squark and gluino masses ^ . The rest 
of the parameters of the squark sector are determined by this setup. For instance, by 
SU{2) gauge invariance we have rriu^ = m^^. Following the same notation as in [43], the 
parameter represents the mixing between the second and third generation of LH-chiral 
squarks. Let us recall its definition: 

2 

fe^^^, (7.3) 

m? _ being the non-diagonal term in the squark mass matrix squared mixing the second 
and third generation of LH-chiral squarks - and an equivalent definition for the up-type 
quarks.^ The parameter ^23 is a fundamental quantity in our analysis as it determines 
the strength of the tree-level FCNC interactions induced by the supersymmetric strong 
interactions, which are then transferred to the loop diagrams of the Higgs boson FCNC 
decays in Eq. ()7.1|) . The last two restrictions in Eq. ()7.2j) ensure that the (heavy) Higgs 
boson decay channels into a pair of squarks are kinematically forbidden. We have checked 
explicitly for some of the heavy Higgs boson channels {h^ can never do it in practice) 
that we obtain the same results if we remove these conditions and include the partial 
widths r{h — > qq*) in the denominator of !^7.1\i . Strictly speaking this condition could be 
implemented, in the case of the boson, by just requiring rriHo < 2Mg-, but we have 
made it stronger by including an additive term. This term is arbitrary (provided it is not 
very small) and acts as a buffer, namely it impedes that by an appropriate choice of the 
squark masses we can approach arbitrarily close the threshold from above, and therefore 
avoids artificial enhancement effects in our loop calculations (see below). Similarly, for the 
CP-odd Higgs boson, the condition expressed in (j7.2j) ensures that we avoid a similar kind 
of enhancement. In this case, however, the condition is a bit different because the A° qq* 

^Our programs are able to deal with completely arbitrary masses for each squark, but we are forced 
to make some simplifications in order to provide a reasonable analysis within a manageable total CPU 
time, see below. 

^Recall that the Sij parameters in the up-sector are related to the corresponding parameters in the 
down-sector by the Cabibbo-Kobayashi-Maskawa matrix, see e.g. [95,96]. 
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vertex can only exist with squarks of different chirality types {A^ qiqn*) or, equivalently, 
with different mass eigenstates (A^qiqj*). We have used fixed values for the soft-SUSY- 
breaking trilinear couplings At and A}, for the bs and to channels respectively. Our results 
are essentially independent of these values and their signs. 

The task of scanning the MSSM parameter space in order to maximize a{pp — s> — > 
qq') for the various Higgs bosons is quite demanding and highly CPU-time consuming, 
even under the conditions imposed in Eq. ()7.2j) . As stated in the introduction, our code 
includes also the restrictions on the MSSM parameter space due to the experimental 
constraint on B{b 57), and therefore contains the full one-loop SUSY-QCD amplitude 
for b s'-f constructed from the FCNC interactions induced by the gluinos. The scan 
was carried out with the help of two entirely different methods. In the first method 
we used a systematic procedure based on dividing the parameter subspace fl7.2|) into a 
lattice which we filled with points distributed in a completely homogeneous way. The 
second is a Monte-Carlo based method, first proposed in [164]. We have adapted the well- 
known Vegas integration program [165] to generate a sufficient number of "interesting" 
points in our parameter subspace. The total number of points used in this case was far 
smaller than in the first method. Obviously the lattice procedure gives more accurate 
results by increasing arbitrarily the total number of points, but the CPU time becomes 
prohibitively long for the whole analysis. This is so even after factoring out in a suitable 
way the phase-space integrals of the Higgs boson production processes, so that these 
integrals are computed only once for every fixed Higgs boson mass and for all the MSSM 
points of our scan in the parameter subspace ()7.2j) . The second method is comparatively 
much faster, but it still involves a quite respectable amount of CPU time for the whole 
analysis. We found that the partial results obtained by the two methods are compatible 
at the level of 10 — 20%. For the study of our FCNC processes we consider that this level 
of accuracy should be acceptable, and for this reason all of the plots that we present in 
this work have been finally computed with the Vegas-based procedure. This also explains 
the wiggling appearance observed in the profiles of the curves presented in Sections 3-4. 
For any given abscissa point in each one of these curves, the corresponding value on the 
vertical (ordinate) axis is somewhere within a band whose width lies around 10 — 20% of 
the central value. 

A few words on the effects of the 5(6 —> 57) constraint in our analysis are now in order. 
The SUSY-QCD contribution to B{b S7) can be quite large, in fact as large as the SM 
one, and with any sign. This raises the possibility of "fine-tuning" between the two type 
of contributions in certain (narrow) regions of the parameter space. As a consequence we 
could highly optimize our FCNC rates in these regions without being in conflict with the 
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experimentally measured B{h ^ 57) band. We have checked that in these regions the 
number of FCNC events can be artificially augmented by one or two orders of magnitude. 
Our scanning procedure indeed finds automatically these fine-tuning domains. However, 
we have systematically avoided them in the presentation of our analysis (for more details 
see chapter EI). In all of our plots, therefore, we show the results obtained for the non- 
fine-tuned case only. We adopt B{h — ^ 57) = (2.1 — 4.5) x 10~^ as the experimentally 
allowed range within three standard deviations [19]. 

7.3 Analysis of the bottom-strange channel 

The main result of the numerical scan for the bottom channel is shown in Fig. 17.11 To 
be specific: Fig. 17.1b shows the maximum values of the production cross-sections a{pp — > 
h ^ qq') for the three MSSM Higgs bosons h = h^, H^, at the LHC, as a function of 
m^o ; Fig. 17.1b displays the cross-section as a function of tan /3. In this plot we indicate the 
value of a{pp ^ h ^ q q') (in pb) in the left- vertical axis, and at the same time we track 
number of FCNC events (per 100 fb~^ of integrated luminosity) on the right vertical axis. 
Looking at Fig. 17. II one can see immediately that at large tan/3: i) the maximum number 
of events is remarkably high (10^ events!) for a FCNC process; ii) there is a sustained 
region in the channel, comprising the interval 300 GeV < rrij^o < 900 GeV, with a fiat 
value of 5 X 10^ events; iii) the chosen value of tan /3 = 50 is not critical for and A^ as 
long it is larger than 10. In Fig. 17.1b we see that the dependence on tan/5 is essentially 
the same for and A^, but for it is quite different: in the region 10 < tan/5 < 30 the 
cross-section remains below 10~^ pb, but for tan /5 > 30 it starts climbing fast up to 0.3 pb 
at tan/5 = 50. The number of events here reaches a few times 10^ for all channels (for 
fixed m^o = 200 GeV) ^. For further reference, in Table ITHI we show the numerical values 
of a{pp ^ h ^ qq') together with the parameters which maximize the production for 
tan/5 = 50 and m^o = 200 GeV. We include the value of B{h — > hs) at the maximization 
point of the FCNC cross-section. We notice that at this point the lightest Higgs boson 
is the one having the smallest branching ratio. This is in contrast to the situation when 
one maximizes the branching ratios independently of the number of events as in chapterlHl 
In Fig. 17.2b we show the effect on a{pp ~* bs), and on the total decay width 

r(/i° —>■ X), of using the Higgs boson sector at the tree- level or at one-loop in our com- 
putation. It is well-known that the couplings to quarks are particularly sensitive to 
this issue, and for this reason we focus on the lightest CP-even Higgs boson for these 

■^As already advertised, in reading the plots and tables in this work, one must keep in mind that they 
are the result of a Monte-Carlo sampling near the region of the maximal values. 
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Figure 7.1: Maximum SUSY-QCD contributions to a{pp — > /i — > 6s), Eg. ()7.1|) . as a 
function of (a) mA_a (at fixed tan/5) and (b) tan/5 (at fixed m^o). In each plot the left- 
vertical axis provides the cross-section in pb and the right-vertical axis tracks the number 
of events per 100 fb~^ of integrated luminosity. 



considerations. We can appreciate the correlations among the different factors that enter 
the production rate (j7.1|) . The plotted values for a{pp ^ bs) and r(/i° X) 

in Fig. 17.2b correspond precisely to the parameters that maximize fl7.1|) at the tree-level 
or at one loop in each case. Fig. 17.2b shows a comparison of the various /i°-pro duct ion 
mechanisms for the values that maximize a{pp bs) at one-loop. Remarkably, the 

effect of the radiative corrections in the Higgs sector amounts to an enhancement of our 
maximal FCNC rates of up to three orders of magnitude. As we mentioned in chapter IB] 
the maximum of B{h^ bs) is attained under the conditions of the so-called "small a^s 
scenario" [157,158], where the two-body decay bb is strongly suppressed due to a 

corresponding suppression of the bb coupling. Since r(/i° bb) usually dominates the 
total width r{h^ — *■ X), the latter also diminishes drastically (at the level of the partial 
width of a three-body decay, as mentioned above). In this scenario the production 
cross-section a{pp h^bb) is also suppressed, so the final result is a compromise between 
the suppression of r(/i° bb) and the possible enhancement (or at least sustenance) of 
a{pp — s> X) by other mechanisms other than the associated production with bottom 
quarks (like the mechanism of gluon-gluon fusion, see Fig. 17.2b ). 
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AO 


a{pp h ^ q q') 


0.45 pb 


0.34 pb 


0.37pb 


events/ 100 fb~^ 


4.5 X 10^ 


3.4 X 10^ 


3.7 X 10^ 


B{h bs) 


9.3 X 10"^ 


2.1 X 10"^ 


8.9 X 10-4 


V{h X) 


10.9 GeV 


1.00 GeV 


11.3 GeV 




10-0.62 


10-1.32 


10-0.44 




990 GeV 


670 GeV 


990 GeV 


A, 


-2750 GeV 


-1960 GeV 


-2860 GeV 




-720 GeV 


-990 GeV 


-460 GeV 


B{b S7) 


4.50 X 10-^ 


4.47 X 10-4 


4.39 X 10-4 



Table 7.1: Maximum value of a{pp ^ h ^ qq') (and of the number of 6s events per 
100/6"^) in the LHC, for rriAo = 200 GeV and tan/3 = 50. Shown are also the corre- 
sponding values of the relevant branching ratio B{h ^ bs) and of the total width of the 
Higgs bosons, together with the values of the SUSY parameters. The last row includes 
B{b S7). 
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Figure 7.2: (a) production cross-section and decay width as a function of m^o with 
the Higgs mass relations at tree-level and at one-loop, (b) Different contributions to the 
production cross-section as a function of m^o corresponding to the maximization of 
a{pp ^ ^ bs) using the one-loop Higgs mass relations. 
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Indeed, for m^o < 300 GeV we can check in Fig. 17.2b that the most relevant factor 
for maximizing the FCNC cross-section is the enhancement of the production channel 
in association with bottom quarks, a{pp — h%h). This channel operates through the hb- 
fusion vertex hh and is highly enhanced at large tan (3. In the region m^o < 300 GeV 
stays as the dominant mechanism for production, although one can see that the alter- 
nate gg-iViSioYi mechanism remains all the way non- negligible. The corresponding effect 
on our FCNC cross-section ()7.1|) is nevertheless not obvious because this same parame- 
ter choice does also maximize the total width of /i^, mainly through the enhancement of 
r(/i'^ ^ hh). There is a delicate interplay of various factors here. In particular, in the 
region m^o < 300 GeV the maximized partial width of the FCNC process —>■ hs (which 
is a function of all the parameters in ()7.2j) ) is larger than in the region m^o > 300 GeV, 
but at the same time the total width becomes smaller in the latter region. Overall, the 
result is that a{pp ^ ^ hs) is larger in the former m^o range than in the latter. 
Furthermore, when we cross ahead the limit m^o ~ 300 GeV the dominant /i°-production 
channel changes turn: the associated /i°-production with bottom quarks falls abruptly 
down (see explanations below) and the gluon-gluon fusion mechanism takes over, so that 
in this range the /i°-production cross-section becomes completely dominated by gg h^. 
We note that while this mechanism is not so efficient at its maximum as the associated 
production, it has a virtue: it is non-suppressed in the entire range of m^o. As a re- 
sult, in the region above m^o > 300 GeV a small value of the h^hb coupling enhances 
B(h'^ hs) while it does not dramatically suppress the total cross-section a{pp — hX). 
For m^o > 300 GeV our sampling procedure finds the maximum by selecting the points 
in parameter space corresponding to the small acs scenario, where B{h^ — > hs) takes 
the highest values and r{h^ — * X) and a{pp — * h^hb) are strongly suppressed - the latter 
staying below the associated Higgs boson production with top quarks! In this way the net 
Higgs boson cross-section a{pp —>■ h^) is not drastically reduced thanks to the sustenance 
provided by the gg-iusion channel. Most of the loop contributions to it come from the top 
quark because the bottom quark contribution is suppressed and the squarks are rather 
heavy. One can clearly see this sustenance feature in Fig. 17.21 in the form of a long cross- 
section plateau up to around m^o = 1 TeV, beyond which the small a^s scenario cannot 
be maintained and r{h^ X) starts increasing and at the same time a{pp ^ ^ hs) 
starts decreasing. It is remarkable that this behavior is only feasible thanks to the large 
radiative corrections in the MSSM Higgs sector. When we, instead, perform the computa- 
tion using the tree-level relations for the Higgs sector, the small scenario is obviously 
not possible and the enhancements/suppressions of a{pp — > h^hb) /T{h^ X) cannot 
take place. As a result the FCNC rate is some 3 orders of magnitude smaller than in the 
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Figure 7.3: Maximum SUSY-QCD contributions to a{pp 
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previous case. 

Next, we turn our attention to the FCNC mixing parameter ^23 in Eq. ()6.5|1 . The 
value of ^23 at the cross-section maximum is not necessarily the maximum allowed value 
of ^23 in ()7.2|) . This is because it is a conditioned maximum, namely a maximum obtained 
under the restrictions imposed hj b 57, as illustrated in Fig. 17.31 For further reference 
in our discussion, and to better grasp some qualitative features of our results, let us write 
the general form of the SUSY-QCD contribution to b ^ 57. If we emphasize only the 
relevant supersymmetric terms under consideration (obviating the powers of the gauge 
couplings and other factors) we have 

B{b ^ S7) ~ Sl,ml{At - /itan/5)VM|usY • (7-4) 

Fig. 17.3b shows the maximum of a{pp ^ h ^ qq') as a function of ^23 for a fixed value 
of the CP-odd Higgs boson mass rriAo = 200 GeV, whereas Fig. 17.3b shows the computed 
value of B{b 57) corresponding to the parameter space points where each maximum 
is attained. At small ^23 the SUSY-QCD contribution to B{b S7) is negligible, and 
the experimental restriction B{b 57) = (2.1 — 4.5) x 10"^ does not place constraints 
on the other MSSM parameters ()7.2|) : in other words, in this region the dependence is 
a(j)p ^ h ^ qq') oc (^23)^ ^ the naively expected one. Here ^23 ^ lO"^'^ ~ 0.03 and 
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Figure 7.4: Value of (a) down-type squark physical masses (M^), four of them are degen- 
erate; (b) the parameters ()7.2|1 from the maximization of the channel in Fig. U.'Ah . 



the computed B{b 57) value lies well within the experimental limit. For larger ^23, 
B{h s'-f) can be saturated at its uppermost experimentally allowed limit, and the 
rest of the parameters in (j7.2j) must change accordingly in order not to cross that limit. 
This can be appreciated in Fig. 17.41 where we show the range of values taken by the 
physical down-type squark masses^ (Fig. 17.4b ) and the lagrangian parameters (Fig. 17.4b ) 
from Eq. ()7.2j) that provide the maximum values of a{pp — > — > 6 s) in Fig. 17.31 In 
the small ^23 region (^23 ^ 10"^'^ ~ 0.03) the parameters and masses that maximize 
a{pp ^ h q q') are constant - except for the statistical noise unavoidable in a Monte- 
Carlo procedure. Note also that there are two possible values for the parameters /i and 
Ah, due to the fact that (the leading contribution of) a{pp ^ h ^ qq') is independent 
of the sign of these parameters and the Monte-Carlo procedure picks either sign for each 
point with equal probability. At ^23 — 10~^'^ the value of B{b 57) becomes saturated 
and a{pp ^ bs) reaches its maximum (cf. Fig. 17.3b ): however ^23 can keep growing, 
yet without overshooting the B{b ^ 57) limits, because the increasing value of ^23 is 
compensated by the growing squark masses (cf. Fig. 17.4b ). But this is not all that simple. 



^There are six difFerent down-type squarks, but four of them are nearly degenerate in mass in our 
approximation. For the 6s-channel, the down squarks are so heavy that the conditions required in the 
last two rows of (|7.2(l are automatically satisfied by them in practically all the allowed range for m^o . 
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the higher range of ^23 can be further divided in two more segments where different 
dynamical features occur. In the first range, namely 10~^'^ < ^ 10~°'^^, the heavy 
Higgs boson channels keep on increasing their FCNC rates, but not so the lightest Higgs 
boson channel a{pp ^ bs), the reason being that for higher squark masses we reach 

the region where the small aes scenario is feasible and hence the couplings become 
weakened. The relevant terms of the cross-section can roughly be written as follows (see 
chapter ^ : 

a{pp ^ /i° ^ 6 s) ~ (t{pp X cos\p - a,s) ml /i'/MlusY , (7-5) 

so that for large tan/3 and small it becomes reduced. In the second high range of ^23, 
i.e. for ^23 ^ 10~°'^^ ~ 0.18, the SUSY mass parameter Msusy has already reached its 
allowed maximum value specified in ()7.2j) . therefore other parameters have to change to 
compensate for the larger 623- This is confirmed in Fig. 17.4b . where for 623 ^ 10~°'^^ the 
absolute value of /i decreases to preserve the B{b 57) upper bound. Correspondingly, in 
this region a{pp ^ bs) further falls down, as it is patent in FigESK- This additional 
feature can also be understood from the approximate expression of the cross-section given 
above. At the same time the FCNC rates for h = H^, keep further growing, but at a 
much lower pace. This is because their (approximate) contribution goes like (see (j6.7|) : 

a{pp ^ (if°, ^ bs) ~ a{pp ^ H^, A^) x 6l, {sm\p - a,s), 1) /iVMgV , (7-6) 
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similar to the case but with angular dependences on acs and /3 which are non- 
suppressing in this region (see chapter E)). Again, any further increase of is now 
partially cancelled by the 6 — > 57 constraint, which demands smaller values of fi. This 
explains the stabilization of the FCNC rates of the and channels in the highest ^23 
range (cf. Fig. 17.3b ). The profile of the squark mass curves in Fig. 17.4b . implies a mixing 
mass matrix with constant diagonal terms (with value Msusy) and growing mixing terms 

(523). 

We finish our analysis of a{pp ^ h ^ qq') by looking at its behavior as a function of 
the SUSY mass scale Msusy, viz. the overall scale for the squark and gluino masses - cf. 
Eq. ()7.2j) . Fig. 17.51 shows the maximum of a{pp ^ h qq') as a function of Msusy for 
fixed m^o = 200 GeV and tan/9 = 50. The interpretation of this figure follows closely the 
results of the previous ones. At small values of Msusy the potentially large contribution 
to B{h s'-f) has to be compensated - see Eq. (|7.4j) - by small values of ^23 and/or 
resulting in a (relatively) small value of a{pp — /i — > qq'). As Msusy grows, ^23 and 
can take larger values without disturbing the restrictions from B{b 57). The leading 
contribution to our FCNC cross-sections is actually independent of the overall SUSY mass 
scale Msusy, because for all Higgs boson channels we have found the general behavior 
(leaving aside other terms mentioned above) 

a(pp ^h^qq')^ a{pp ^ h) x 5l, /xV^Iusy (h = h\ H', A') . (7.7) 

The last factor effectively behaves as (J^a A^^/^Iusy' grows as for increasing Msusy 
at fixed ratio /i/MsusY- Under the same conditions B{b 57) causes no problem because 
it is additionally suppressed by m^/MgusY ^ cf. Eq. (j7.4|) . Therefore, we are led to a sort of 
"non-decoupling behavior" of the FCNC rates with increasing Msusy- In other words, we 
find that for the heavy neutral Higgs bosons {H^, A^) the interesting region is (contrary to 
naive expectations) the high Msusy range! The lightest Higgs boson (/i°) channel shows 
a similar overall behavior, but it presents additional features because it is more tied to 
the evolution of the CP-even mixing angle a. The most interesting region for this channel 
is (cf. Fig. 17. 5p the central squark mass scale Msusy ~ 600 — 800 GeV, where the small 
aeflf scenario can take place. 

From the combined analysis of Figs. 17.1117.51 we arrive at the following conclusions 
concerning the bs final state: 

• A significant event rate of FCNC Higgs boson decays a{pp ^ h ^ qq') is expected 
at the LHC, even after taking into account the limits on B{b S7); 

• Lightest Higgs boson case. 
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— For rriAo < 300 GeV the rate a{pp ^ ^ bs) decreases with 771^0 but it is the 
largest in this interval, being produced by the combination of a large production 
cross-section a{pp — h^bb) and a moderate B{h^ — » bs). It amounts to a 
number of events between ~ 5 x 10^ and ~ 12 x 10^ for every 100 fb~^ of 
integrated luminosity at the LHC; 

— For 300 GeV < m^o < 850 GeV we expect amaximumof ~ 6x10^ events/100 fb 
in the small aes scenario, provided by a large B{h^ ^ bs) and (y{pp{gg) h^) 
as the dominant production cross-section; 

— For m^^o > 850 GeV the number of events starts to decrease slowly; 

— In all cases, this maximum is attained for a large value of tan/3 ~ 50, a 
moderate value of the SUSY mass scale (Msusy ~ 600 — 800 GeV) and a low 
value of ^23 ~ 10-^-^ ~ 0.05; 

• Heavy Higgs bosons, H^, 

— Although not shown in our plots, we have checked that their production rate 
a{pp A^) decreases fast with the Higgs boson mass (due to the decreasing 
of the production cross-section). We find a maximum FCNC rate of ~ 5 x 10^ 
events for m^o ~ 200 GeV, and 20 events for m^io ^ 1 TeV. 

— The maximum is produced at i) large tan/5 > 30, ii) at the highest allowed 
values of the SUSY mass scale, Msusy ~ 1 TeV, and iii) at a relatively large 
value of the FCNC mixing parameter, ^23 ~ 10"°'^^ ~ 0.18, but not at the 
largest allowed value. The small aes scenario plays no role in the heavy Higgs 
boson channels. 

Altogether one should expect a total maximum of some 120,000 events/100 fb~^. 

7.4 Analysis of the top-charm channel 

The results of the numerical scan for this channel are similar to the bs channel, so we will 
focus mainly on the differences. Fig. 17.6b shows the maximum value of the production 
cross-section a{pp h tc), Eq. (j7.1|) . under study as a function of rrij^o; Fig. 17.6b 
displays the cross-section as a function of tan /5. Obviously the lightest Higgs boson (/i^) 
channel does not appear in these plots, since in the MSSM this boson is always lighter 
than the top quark. Looking at Fig. 17.61 one can see immediately the following: i) the 
dominant channel in this case is the heavy scalar Higgs boson, H^; ii) it varies between 1 
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(a) (b) 

Figure 7.6: Maximum SUSY-QCD contributions to cr^pp ^ h ^ tc), Eg. ()7.H) . as a 
function of (a) rriAO (at fixed tan/5) and (b) tan/5 (at fixed m^o). 



and 300 events/100 fb~ ; iii) the tan/5 value is critical with preference for low values. In 
Table 17^ we show the numerical values of a{pp ^ h ^ to) together with the parameters 
which maximize the production for tan /5 = 5 and niAo = 300 GeV. We have included the 
value of B{h ^ to) at the maximization point of the FCNC cross-section. It is remarkable 
that for the heavy CP-even Higgs boson one can reach B{H^ to) ~ 10"'^ compatible 
with the b —>■ S'-f constraint. 

Let us remark that for the heavy Higgs boson channels the features of the small a^s 
scenario play no significant role because the partial widths into bb are proportional either 
to cos^Ocfr (in the case) or to sin^/5 (A^ case). Moreover, in the low tan/? > 1 region 
(the relevant allowed one for the to channel) the small a^s scenario does not even have a 
chance to take place. 

We turn now our view to the role of the B{b s'-f) restriction in Figs. 17.71 and 17.81 
Fig. 17. 71 shows the maximum value of a{pp /i ^ t c) as a function of ^23 for a fixed value 
of m^o = 300 GeV together with the corresponding computed value ol B{b ^ s'-f), while 
Fig. 17. 81 shows the values of the parameters, Eq. ()7.2|) . that realize this maximum, together 
with the physical up-type squark masses. In this case the B{b — 57) restriction is not as 
critical as in the bs channel, in part due to the fact that the SUSY-QCD contribution to 
B{b S'-f) is not enhanced at low tan/5. For ^ 10^^'^ the value of B{b s'-f) is well 
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a[pp — > h — *• t c) 


2 A X 10"3pb 


5.8 X 10"^ pb 


events/100 fb"^ 


240 


58 


B(h tc) 


1.9 X 10"^ 


5.7 X 10"^ 


Y(h X) 


0.41 GeV 


0.39 GeV 


59 


10-0.10 


10-0.13 




880 GeV 


850 GeV 


At 


-2590 GeV 


2410 GeV 


fi 


-700 GeV 


-930 GeV 


B{b 57) 


4.13 X 10"^ 


4.47 X 10"^ 



Table 7.2: Maximum value of a{pp ^ h ^ tc) (and of the number of tc events per 
100/6"^) in the LHC, for m^o = 300 GeV and tan/3 = 5. Shown are also the corre- 
sponding values of the relevant branching ratio — > tc) and of the total width of the 
Higgs bosons, together with the values of the SUSY parameters. The last row includes 



inside the experimental limits fFig. 17.7b ). there is no restriction on the rest parameters 
(Fig. 17.8b ). the squark masses remain constant (Fig. 17.8b ). and the maximum value of 
a{pp h tc) grows here in the naively expected way (^23)^ (Fig- 17.7b ). Above this 
value (523 ^ 10~^'^) the parameters have to be adjusted to provide and acceptable range 
for B{b — » S7) ^. In this region the SUSY mass Msusy grows fFig. l7.8b ). and decreases, 
but not so fast as in the bs channel case fFig. 17.4b ). At the same time At increases with 
increasing ^23 > 10~^'^. Most of the physical squark masses grow, but one of the up-type 
squarks (stop squark) can always have the minimum allowed mass - Eq. ()7.2j) . In this 
region the observables under study grow more slowly, since their original behavior 
is partially cancelled by the growing of MgusY with 623- In Fig. 17.91 we see, again, that 
a{pp ^ h ^ tc) grows with the SUSY mass scale (although in a way less pronounced 
than in the bs channel), due to the relaxation of the b ^ s'-f constraint for large Msusy- 
While a{pp h ^ qq') is augmented nearly three orders of magnitude in the range 
^susY = 200 — 1000 GeV fFig. 17. 5|) . the tc channel undergoes only an increase of roughly 



^The two lines appearing in this region for Bnib — > 57) mean that the maximum of cr(pp —f H'^ tc) 
is attained either by the maximum or the minimum allowed value of 5(5 ^ 57), our Monte-Carlo sampling 
procedure picks either choice with equal probability for each value of 52z- 
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Figure 7.7: Maximum SUSY-QCD contributions to a{pp ^ h ^ to), Eg. (|7.1|) . as a 
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Figure 7.8: Value of (a) up-type squark physical masses (M^), four of them are degenerate; 
(b) the parameters ()7.2j) from the maximization of the channel in Fig. 17.7b . 
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a factor 10 in the same parameter range. 

How do the two Higgs boson channels and compare as sources of tc events? 
The (y'^'-fusion mechanism is one of the leading processes for Higgs boson production at 
relatively small values of tan /5 > 1 - associated production with bb remaining still sizeable. 
Due to CP conservation, the squark contributions to gg —>■ cancel out at one-loop and 
only the quark contributions remain [136]. For large squark masses the production cross- 
sections for and A^ are similar, the latter being slightly larger. We see from Fig. l7.6l 
that, despite the similarity in production, the channel gives larger FCNC rates than 
the A^ one. The excess of FCNC events from the former can be explained mainly from 
the constraints that we have imposed from the very beginning on the squarks masses 
in relation to the Higgs boson masses (see Eq. ()7.2|) ). As we have noted in section 17^ 
we can have squarks of the same chirality-type in the qq* vertex, whereas they must 
necessarily be of opposite chirality-type in the A^qq* case. As a result, for small m^o the 
squark mass constraints expressed in ()7.2|) allow the FCNC cross-section maximization 
process to pick points near the saturation of the mass condition 2 Mg- > muo + 50 GeV, 
but not of Mq- + Mq. > m^o + 50 GeV for {i ^ j) because only one of the up-squarks 
can be light (see Fig. 17.8b ). This produces an enhancement of the branching ratio of 
tc and for this reason this FCNC channel dominates in the relatively small m^o 
region. However, as soon as m^o is sufficiently heavy the second mass constraint can also 
be satisfied and then the two curves in Fig. 17.61 tend to converge. 

From the combined analysis of Figs. I7.6II7.9| we conclude that in the case of the 
channel we expect a maximum of ~ 300 events/100 fb~^ decays into top quarks at the 
LHC. This maximum is achieved for a CP-odd Higgs boson mass of m^o ~ 300 GeV, and 
a moderately low tan /5 ~ 5. This rate can grow one order of magnitude by a lower value 
of tan/5 ~ 2, but decreases significantly with m^o. The maximum is obtained at the 
largest possible value of ^23, and a moderate SUSY mass scale Msusy ~ 600 — 800 GeV, 
but having one of the squarks light. While the number of events is significantly lower 
than the 6s-channel ones, the tc-channels offers a better opportunity for detection, due to 
the much lower background. 

7.5 Discussion and conclusions 

We have carried out a systematic study of the production rate of FCNC processes at the 
LHC mediated by the decay of neutral Higgs bosons of the MSSM: cr{pp h ^ q q') {h = 
h^, H^, A^) - see Eq. 1)7.11) . Specifically, we have concentrated on the FCNC production of 
the heavy quark pairs qq' = bs and tc, because they are the only ones that have a chance of 



124 Production and FCNC decay of MSSM Higgs bosons into heavy quarks in the LHC 



J2 



10" 





/ 



— ct(pp -> H -> t c) 

— - ct(pp -> a -> t c) 



,-5 



tanP = 5 

m,o = 300 GeV 



10" 



200 



400 



600 



800 



M, 



SUSY 



[GeV] 



Figure 7.9: Maximum SUSY-QCD contributions to a^pp ^ h ^ tc), Eq. (jTHj), as a 
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being detected. We have focused on the FCNC supersymmetric effects stemming from the 
strongly interacting sector of the MSSM, namely from the gluino-mediated flavor-changing 
interactions. We have performed a maximization of the event rates in the parameter space 
under a set of conditions that can be considered "irreducible", see Eq. fl7.2|l . i.e. we cannot 
further shorten this minimal set (e.g. by making additional assumptions on the relations 
among the parameters) without potentially jeopardizing the conclusions of this study. 
Even within this restricted parameter subspace the computer analysis has been rather 
demanding. The numerical scan has been performed using Monte Carlo techniques which 
we have partially cross-checked with more conventional methods. The maximization of the 
cross-sections fl7.1|) has been performed by simultaneously computing the corresponding 
MSSM quantum effects on the (relatively well-measured) low-energy FCNC decay b ^ s'y 
and requiring that the experimental limits on this observable are preserved. 

To summarize our results: the total number of FCNC heavy flavor events originating 
from supersymmetric Higgs boson interactions at the LHC can be large (of order 10^), but 
this does not mean that they can be easily disentangled from the underlying background 
of QCD jets where they are immersed. For example, it is well known that the simple two- 
body decay h ^ bb is impossible to isolate due to the huge irreducible QCD background 
from bb dijets - a result that holds for both the SM and the MSSM [13,14]. This led a 
long time ago to complement the search with many other channels, particularly /i — 7 7 
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which has been identified as an excellent signature in the appropriate range. Similarly, 
the FCNC Higgs boson decay channels may help to complement the general Higgs boson 
search strategies, mainly because the FCNC processes should be essentially free of QCD 
background. Notwithstanding other difficulties can appear, such as misidentification of 
jets. For instance, for the bs final states misidentification of 6-quarks as c-quarks in cs- 
production from charged currents may obscure the possibility that the 6s-events can be 
really attributed to Higgs boson FCNC decays. This also applies to the to final states, 
where misidentification of 6-quarks as c-quarks in e.g. tb production might be a source 
of background to the to events, although in this case the clear-cut top quark signature 
should be much more helpful (specially after an appropriate study of the distribution of 
the signal versus the background). However, to rate the actual impact of these disturbing 
effects one would need an additional study which is beyond the scope of the present work. 

An interesting (and counter-intuitive) result of our work is that for all the Higgs boson 
channels h = h^, H'^, A^, the FCNC cross-section a{pp —^h^qq') increases with growing 
SUSY mass scale Msusy- Due to this effective "non-decoupling" behavior (which is more 
pronounced for the bs channel) the FCNC rates are maximal when the overall squark and 
gluino mass scale is of order of Msusy ^ 1 TeV - with the only proviso that for the tc- 
channel a single squark should have a low mass (> 150 GeV). Moreover, we find that the 
two types of FCNC final states {bs and to) prefer different ranges of tan j3. The 6s-channel 
is most efficient at high tan f] > 30, whereas the tc-channel works better in the regime of 
low tan/3 < 10 (see below). As for the mixing parameter (which is the fundamental 
supersymmetric FCNC parameter of our analysis, see its definition in ()(i.5|) ) we remark 
that the maximum number of events is not always attained for the largest possible values 
of it (due to the infiuence of the b —>■ s'-f constraint): in the 6s-channel the maximum is 
achieved for moderate values in the range 0.05 < ^23 ^ 0.2, whereas the tc-channel prefers 
the maximum allowed values in our analysis (^23 ^ 0.8). We also remark that the naive 
expectation a{pp h ^ q q') ^ does not always apply. 

The number of FCNC events originating from the two channels bs and tc is not alike, 
and it also depends on the particular Higgs boson. For the bs final states we have found 
that the optimized value of (j{pp ^ h ^ qq') produced by our analysis is ~ 12 pb. This 
amounts to ~ 12 x 10^ events per J Cdt = 100 fb~^ of data at the LHC. The most favorable 
Higgs boson channel is the one corresponding to the lightest MSSM Higgs boson, h^. For 
this boson there are non-trivial correlations between the two factors in Eq. 1)7.11) . namely 
between the Higgs boson production cross-section and the FCNC branching ratio. These 
correlations permit an increase of the total number of FCNC events up to two orders of 
magnitude in certain cases as compared to the number of events produced by the heavy 
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Higgs bosons and A^, which are essentially free of these correlations. The latter stem 
from relevant quantum effects on the parameters of the Higgs boson sector at one-loop 
precisely in the regions of our interest. 

On the other hand, the maximum value of a{pp — ^ /i — * tc) is more moderate, to 
wit: 3 X 10~^pb, or ~ 300 events/100 fb~^. For the total integrated luminosity during 
the operative lifetime of the LHC, which amounts to some (300 — 400) fb~^, we estimate 
that a few thousand to events could be collected in the most optimistic conditions. This 
number is of course sensitive to many MSSM parameters, but most particularly to two: 
m^o and tan j3. The mass of the CP-odd Higgs boson should not be heavier than m^o ~ 
(400 — 500) GeV if one does not want to decrease the number of events below a few hundred 
per 100 fb~^. On the other hand the number of events is very much dependent on the 
particular range of tan (3. As we have said above, the lowest possible values are preferred 
for the to channel, but the sensitivity in this range is so high that the order of magnitude 
of a{pp —^h-^tc) may change for different (close) choices of tan/?. Throughout all the 
analysis of the to channel we have fixed tan (3 at an intermediate value, but the maximum 
number of events per 100 fb"^ would grow from ~ 300 (for our standard choice tan (3 = 5) 
up to ~ (500, 900, 2000) if we would have chosen tan P = (4, 3, 2) respectively. Such 
lower values of tan jS are usually avoided in some MSSM analyses in the literature, but 
as a matter of fact there are no fully water-tight experimental bounds on tan j3 excluding 
this lower range, apart from the more incontrovertible strict lowest limit tan/3 > 1. We 
recall that the lower limit on tan (3 is obtained indirectly from the LEP exclusion data 
on the light Higgs boson search, and is therefore very sensitive to the inputs used in the 
computation, specially the top quark mass [108]. Unlike the difficulties in the 6s-channel, 
and in spite of the substantially smaller number of events, we deem more feasible to 
extract the to signal at the LHC, due to the presence of the quark top, which carries a 
highly distinguishable signature. 

At this point a comparison with the previous chapter is in order. In that work we 
have studied in quite some detail the maximum FCNC production rates of Higgs bosons 
decaying into to final states within the general two-Higgs-doublet model. It was found that 
the maximal branching ratio in the 2HDM takes place in the type II 2HDM (or 2HDM 
II), and reads B^\h to) ~ 10~^, whereas in the 2HDM I it is comparatively negligible. 
After a detailed computation of the event rates (including also the particular restrictions 
of the b s'-f process, which are different in the 2HDM case as compared to the MSSM) 
the conclusion was that several hundred to events could be collected at the LHC under 
optimal conditions. We clearly identified which are the most relevant Higgs boson modes 
for this purpose and the domains of the 2HDM II parameter space where these events 
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could originate from. To make it short, our conclusion in chapter |S1 was that the is 
the most gifted decay and the ideal situation occurs when tan/5 and tana are both large, 
and also when the CP-odd state A'^ is much heavier than the CP-even ones 
Furthermore we found that in the general 2HDM the state never gives any appreciable 
FCNC rate into to. It is easy to see that the mode bs is not favored either (unless 

tan/5 is very small and unusually light, both situations rather unappealing). 

How it compares with the MSSM case under study? To start with, we note that here 
the A^ channel gives essentially the same bs rate as the one, and that both modes 
can be quite relevant. At the same time the A^ rate into tc is, though not dominant, not 
negligible at all. Some few hundred events of this nature per 100 fb~^ are possible. If this 
is not enough, in the MSSM the most relevant tan (3 region for the tc final states is not the 
highest one (as in the 2HDM II case) but just the opposite: the lowest allowed one. This 
is an important difference, and one that should help to discriminate between the 2HDM 
and MSSM models in case that some tc events would be unambiguously tagged at the 
LHC. After all there are many high precision observables that are highly sensitive to the 
preferred range of tan (3, so that the favorite value of tan /5 could already be fixed from 
other experiments by the time that some FCNC events could be detected. Apart from 
the different correlation of the parameters in the non-supersymmetric and supersymmetric 
model, in the latter case the maximal event rates for the tc mode are typically one order 
of magnitude higher than the maximal event rates in the former. 

The corresponding study for the 2IIDM branching ratios into bs final states was per- 
formed in [51], although in this work the cross-section and number of events were not 
computed. However, from the maximum size of the expected branching ratios compatible 
with b s'j (viz. B^\h bs) ~ 10"^ - 10"^ in the 2HDM II) it is already pretty 
obvious that the number of events can never be competitive with the supersymmetric 
case where B'^^^^^^{h bs) ~ lO""' — 10~^as seen in chapter IHl As for the 2HDM I 
(which is insensitive to the b —* s'-f bounds) the branching ratios can be at most of order 
B^h bs) ~ 10~^ — 10"'^ and only so for very small values of tan/5 = 0.1 — 0.5 which 
are actually excluded in the MSSM. 

Summing up and closing: the maximum number of FCNC events in the MSSM case 
is larger than the highest expected rates both in the 2HDM I and II , the two kind of 
signatures of physics beyond the SM being perfectly distinguishable because the relevant 
regions of the parameter space are completely different. If a sample of FCNC events of 
this kind could be collected, we should be able to ascertain which is its ultimate origin. At 
the end of the day if one single thing should be emphasized is that the FCNC event rate 
into bs or tc is so extremely tiny in the SM that if only a dozen events of this kind could be 
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captured under suitable experimental conditions it would be an undeniable signature of 
new physics. From what we have seen, the odds should be heavily in favor of attributing 
it to a supersymmetric origin. 



Chapter 8 
Conclusions 



In this PhD. Thesis we have presented a monographic study of some rare Flavour Chang- 
ing Neutral Currents (FCNC) induced in the Two-Higgs-Doblet-Models (of both kinds 
2HDM I and II) and in the Minimal Supersymmetric Standard Model (MSSM). In par- 
ticular we have studied the FCNC top quark decay in the 2HDM, t ch, Higgs boson 
production and decay in the 2HDM, pp ^ h + X ^ tc, h tc, and the Higgs boson 
production and decay in the MSSM, pp ^ h + X ^ tc,bs and h — * tc, bs, in the context 
of the LHC, where h = h^, H^, are the three neutral Higgs bosons in these extended 
Higgs models. 

The main results obtained are the following: 

• The non standard effects on the top quark FCNC decay into Higgs boson in the 
2HDM II can be very important, at most ten orders of magnitude larger then in the 
SM, bringing the branching ratios at the level of 10~^ or higher in some cases. Thus 
there is a real chance for seeing rare decays of that sort at the LHC. The effects of 
the 2HDM I, although sizeable, do not reach the visible level at the LHC. 

It is useful to compare the t he processes with the more conventional rare decay 
t — s> gc. The maximum rates for the leading FCNC processes ()4.H) and ()4.2|) in the 
2HDM II (resp. in the MSSM) satisfy the relations 

B{t ^gc)< 10"^(10~^) <B{t^hc)^ 10"\ (8.1) 

where it is understood that h is or iJ", but not both, in the 2HDM II; whereas h 
is most likely but it could also be and in the MSSM. These decay ratios 
are compatible with their observation at the LHC if we compare them with the 
estimations of the sensitivities for 100 fh~^ of integrated luminosity in the relevant 
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colliders [25,133,166]: 

LHC -.Bit cX) > 5 X 10"^ 
LC:B{t^cX)>5xlO-\ (8.2) 
TEV33 -.Bit cX) > 5 X 10"^ . 

The origin of the FCNC effects from the general 2HDM contributions is purely 
electroweak (namely it stems from enhanced trilinear Higgs boson couplings), in 
contrast to the MSSM case where they are mainly from the strong (gluino) sector, 
with the electroweak contributions one-two orders of magnitude below. This can 
be important because the MSSM strong couplings are very weakly restrained from 
experiment. 

Although the maximum ratio is similar there are different strategies to distinguish 
between the 2HDM and MSSM FCNC Higgs boson effects in the case that these 
decays would be detected. A possible strategy would be to look for different sig- 
natures for the process t ^ he ebb. In the favorite FCNC region ()4.7|) of the 
2HDM II, the combined decay t ^ h e ^ ebb is possible only for or for 
but not for both - Cf. Fig. I4.3fe i - whereas in the MSSM, together with if", are 
highlighted for llQGeV < rrij^o < rrit, with no preferred tan/? value. And similarly, 
t ^ e is also non-negligible for m^o < 120 GeV [43]. 

• The potential enhancement of the 2HDM contributions to the Higgs boson decay into 
top and charm quarks may reach up to ten billion times the SM value B{H^^^ — > 
tc) ~ 10~^^, thereby bringing the maximum value of the FCNC branching ratio 
B{h^ ^ tc) to the level of ~ 10~^. Adding to this the 2HDM contributions of the 
Higgs production [pp H+X) we find that the events pp H+X tc+X could 
be visible at the LHC, getting a few hundred events for the and for the H^, but no 
for the A^. Thus, theses events can be very useful to complement the Higgs boson 
detection strategies, like for example those based on the decay — > bb (affected 
with a huge QCD background) or on the decay — s> 77 (more promising). 

There are some strategies to distinguish between the 2HDM and the MSSM. An 
obvious one is the detection of the channel A^ — * tc that is highly suppressed in 
the 2HDM but it is not so in the MSSM. Another possibility is using the fact that 
the MSSM contribution to FCNC processes is mainly from the gluino, which is not 
particularly sensitive to tan f3, within a range of high and moderately values of this 
parameter. If, for example, a few of these events were observed and at the same 
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time the best MSSM fits to the electroweak precision data would favor moderate 
values of tan/3, say in the range 10 — 20, then it is clear that those events could 
originate from the FCNC gluino interactions but in no way within the context of 
the general 2HDM. 

• The SUSY-QCD contributions can enhance the maximum expectation for the FCNC 
decay rates B{h —* qq'), where qq' = tc or bs, reaching the level of 10"'^, in particular 
^max^/^o ^ ^^/) _ 3 X 10-3 and 5'^^^(if7A° ^ gg') ~ 9 X 10"^. This corresponds 
to an scenario where r(/i° bb) is suppressed by radiative corrections to the CP- 
even mixing angle a (the so-called small ae// scenario). When it is not the case the 
maximum ratios give the same result except for the light Higgs boson B^^^[h^ 
qq') ^ 1.3 X lO"'^. So we get branching rations three to nine orders of magnitude 
higher than in the SM case, depending on the qq' channel. 

The total number of FCNC heavy flavor events (tc, bs) originating from supersym- 
metric Higgs boson interactions at the LHC can be large (of order 10^ for bs), but 
this does not mean that they can be easily disentangled from the underlying back- 
ground of QCD jets where they are immersed. For the bs flnal state the maximum 
number of events is ~ 12 x 10^ for a 100 fb~^ of integrated luminosity at the LHC, 
being the most favorable Higgs boson channel the lightest MSSM Higgs boson, . 
There are some correlations in this channel that increase the number of events up 
to two orders of magnitude in certain cases as compared with the other channels 
, . On the other hand, for the case of tc the maximum number of events is 
~ 300 events/100 fb~^ for tan/3 = 5, which amounts to few thousand tc events dur- 
ing the lifetime of the LHC. This value is sensitive to the rest of the parameters, 
in particular m^o and tan/3. The mass of the CP-odd Higgs boson should not be 
heavier than m^o ~ (400 — 500) GeV if one does not want to decrease the number of 
events below a few hundred per 100 fb~^. On the other hand the number of events is 
very much dependent on the particular range of tan /3, with a preferred low values, 
and changing orders of magnitude with tiny changes in tan (3. For example for the 
values tan/3 = (5,4,3,2) we get a number of events ~ (300,500,900,2000). Al- 
though such lower values are usually avoided in the literature, they are not strictly 
excluded experimentally. 

We emphasize the primary differences with the 2HDM. In the MSSM, (1) the 
and channels give essentially the same bs rate, and it can be quite important; 
(2) the A^ rate into tc is not negligible; (3) the relevant range in tan (3 is the lowest 
allowed one. In general the maximum number of FCNC events in the MSSM case 



132 



Conclusions 



is larger than the highest expected rates both in the 2HDM I and II, the two kind 
of signatures of physics beyond the SM being perfectly distinguishable because the 
relevant regions of the parameter space are completely different. 

Although we have shown that the FCNC events considered in this Thesis are in princi- 
ple visible (namely, we have proven that the branching ratios and production cross section 
are sizeable enough as compared to other rare, but detected, events), it should be empha- 
sized that the final word can only be said after an appropriate signal versus background 
experimental study, which is out of the scope of this theoretical investigation. However, 
we consider highly encouraging that there is no theoretical a priori obstruction for these 
processes not to be visible in the experimental setup of the LHC collider. In particular the 
tc + ic events, with the top quark in the final state, should be very helpful and essentially 
free from background. 

From these positive theoretical results we can assert that the pathway to seeing new 
physics through FCNC decays of the top quark and Higgs boson is thus potentially open. 
It is now an experimental challenge to accomplish this program using the high luminosity 
super-colliders round the corner, in particular the Large Hadron Collider at CERN, which 
is scheduled to start operating within two years from now, i.e. in 2007. 



Appendix A 
Vertex functions 



In this appendix we briefly collect, for notational convenience, the basic vertex functions 
frequently referred to in the text. In practice we have performed the calculations using 
the algebraic and numerical programs FeynArts, FeynCalc and Loop Tools [128-130]. The 
given formulas are exact for arbitrary internal masses and external on-shell momenta. 

Most of them are an adaptation to the g^j^^ = {~\ } metric of the standard formulae 

of Refs. [167-169]. The basic one-, two- and three-point scalar functions are: 

Ao{m) = [ <rq ^ (A.l) 
J [g^ — m^] 

5o(p,mi,m2) = / JTTT"^ — ^2 2T ' (^'2) 

J [q^ - mf\ [{q + py - 



Co{p,k,mi,m2,m3) = I d^'q— ^ — — ^ ; (A.3) 

[q^ — m{\ [[q + pY — 7712] [\q + p + ky — mi^\ 



using the integration measure 



The two and three-point tensor functions needed for our calculation are the following 



\B„,B,,B,Mp„) = / (A.5) 



[CQ,C^,C^y\{p,k,mi,m2,m^) = 

[g^ (lt^, %(lu\ 



(Tq 



[q^ — mf] [(g + pY — ml] [{q + p + k^ — m| 



(A.6) 
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By Lorentz covariance, they can be decomposed in terms of the above basic scalar func- 
tions and the external momenta: 

Bo{p,mi,m2) = Ao{m2) + mlBo{p,mi,m2) , 
B^{p, nil, 1712) = p^Bi{p, nil, 1712) , 

Bf,^{p,mi,m2) = p^iPuB2i{p,mi,m2) + g^,uB22{p,mi,m2) , 
Co{p,k,mi,m2,m3,) = Bo{k,m2,ms) + mlCo{p, k,mi,m2,m3) , 
C^{p, k, mi, m2, mg) = p^Cu + k^,Ci2 , 

Cf,y{p, k, mi, 1712, ms) = p^PuC2i + k^ki,C22 + {Pf^k^ + k^pu)C23 + gf,uC24 , (A. 7) 
where we have defined the Lorentz invariant functions: 

Bi{p, nil, 1712) = ^[Ao{nii) - Ao{ni2) - fiBo{p,nii,ni2)], (A. 8) 

B2i{p,nii,ni2) = —-^ — —[{n-2)Ao{ni2)-2nilBo{p,nii,ni2) 
Ip'^yn — 1) 

- nfiBi{p,nii,ni2)], (A. 9) 

522b, mi, ma) = _ [Ao{m2) + 2mlBo{p, nii, ma) +/i^i(p, mi, ma)], (A. 10) 

j _ y I Bo{p+ k,nii,ni3) - Bo{k,ni2,nis) - fiCo , (All) 
C12 J \ Bo{p,nii,ni2) - Bo{p + k,nii,m3) - f2Co 

^'^^\=y( ^i(^' + ^'"^1'"^3) +5o(^,"^2,m3) -/iCii -2C24 , (A 12) 

C23 I \ Bi{p,mi,ni2) - Bi{p + k,mi,m3) - f2Cii 



C22 = TTTTT^ / 1 \9i {-pk[Bi{p + k, mi, ma) - Bi{k, ni2, m^) - fiC_ 



12 



2[p2A;2 - (pk) 

p^[-Bi{p+k,nii,m3) - /2C12 - 2C24]}, (A.13) 



C24 = —^-—[Bo{k, ni2, mg) + 2m2Co + fiCn + /2C12] , (A.14) 
2[n — 2) 

the factors /i^2 and the matrix Y, 

r 2,2 2 

Ji = p + m^ - m^. 



/a = k'^ + 2pk + m2 — m 



3 ' 



^ 2[p2p _ (p^)2] I p2 1 • (A-IS) 
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The UV divergences for n — 4 can be parametrized as 



e 

A 



n - 4 
2 



+ 7s - ln(47r) 



(A.16) 



being 7^; the Euler constant. In the end one is left with the evaluation of the scalar 
one-loop functions: 



Ao{m) 
Bo{p,mi,m2) 



2 

A + In ^ - 2 + ln[(xi - l)(x2 - 1)] 



(A.17) 



+Xi m h a;2 m 



Co{p, k, nil, 1712,1713) 



Xi — 1 
-i \ 1 1 



X2-I 



with 



2:1,2 = Xi^2{P,mi,17l2) 



1 22 
1 m| — m2 



(A.18) 
(A.19) 



2^2 ±^AV2(/,^2^^2)^ (A.20) 



K^,y,z) = [x - {y/y - y/zf][x - {y/y + ^/zY] , 

and where S is a bookkeeping device for the following alternate sum of twelve (complex) 
Spence functions: 
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The Spence function is defined as 

Sp{z) 

and we have set, on one hand: 



^ ln(l - zt) 



dt , 



(A.22) 



-1,2 

h,2 



^1,2 



xi,2{p, rn2,mi) , 
Xi^2{p + k,m3,mi] 
xi,2{,k, 1713, 1712) ; 



(A.23) 



136 



Vertex functions 



and on the other: 

yi = yo + ^, 2/2 

where 



1 - e ' e ' ^° 2 pk + p^^ 



g = —k"^ + — m| , = —p^ — 2pk — + m\ 
and ^ is a root (always real for external on-shell momenta) of 



+ 2pki + = {). 



(A.24) 



(A.25) 



(A.26) 



Derivatives of some 2-point functions are also needed in the calculation of self-energies, 
and we use the following notation: 

d_ 

dp"^ 



—B^{p,mi,m2) = B'^{p, nil, 1712). 



(A.27) 



We can obtain all the derivatives from the basic B'r,: 



1 



Xi{xi — 1) In 



p2 A^/2(p2, m^, m|) 
Xi — 1 



- X2{x2 - 1) In 



X2-I 



,(A.28) 



Xi / \ X2 

which has a threshold for \p\ = mi + m2 and a pseudo-threshold for \p\ = \mi — m2\- 



A.l Limit of heavy internal masses 

Next we will present the results when the external moments are much heavier than the 
internal masses {p^ -C m?). The basic equations is: 



1 



1 



{q + pY — m? — m? 

1 



p^ + 2p ■ q 
q2 _ 



g2 _ ^2 



1 - 



p^ + 2p ■ q f p^ + 2p ■ q 



q2 _ j^2 



q2 _ j^2 



(A.29) 



which we have to substitute in each of the denominators of the n-point functions. But, 
all the lineal terms in the expansion become zero under the integrations j d^q, because 
of the denominator is a spherically symmetric function. The others terms are integrals of 
functions of q^ and p^ /{q^ — im?Y which are negligible in this limit. 
Thus the effective substitution is: 

(A.30) 



(g -|- pY — ni? q^ — w? ' 
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in the cases where <^ m^. 

So in this hmit the ra-point functions do not depend on the moments. It allows 
us to define, without ambiguities, that the n-point functions without the moments are 
proportional to the complete ?7,-point functions with the moments equals to zero. As: 



Co(0,0,mi,m2,m3) 



167r2 



Co(mi,m2,m3) 



(A.31) 



with 



Co (mi, 7712,7713) 



777T — rrio 



rF(77ii, 7772, 7773) 



7-1/ \ 1 rrii lltin llbn 

F(777i, 7772, 7773) = log— H 2 2^^^~ 2 2 ^"^S 2 



777q 



777o 



7772 '^S ~ '^1 "^3 ~ "^2 "^2 

F(7772, 77^1, 77^3) . 



(A.32) 

(A.33) 



One can see that if, for example, 77i^ ^ m\^m\, then the function ()A.3H) only depends 
on this heavy mass as Co oc The same is true for 7773 ^ 777^, 777| and 7773 ^ 777^, 7773. 
The function C^,y(77ii, 7772, 7773) is logarithmically divergent and can be easily written 

as: 



C^^(77ii, 7772, 7773) = Qpiv 



1 1 

-5o(77i2,77i3) + -777^Co(77Zi, 7712,7773) 
77 4 



(A.34) 



But we can not substitute directly the value of Bq because is divergent and divided by n. 



77io 



Bo{m2,ms) = A - 1 + logm^ ^ 2 



7775 



mo — m 



mo 



lim —Bninio, rri'^) = lim - (l 



A + Bo(m2,m3) 



finite 



A - - + 5o(m2,m3) I finite 



(A.35) 
(A.36) 



where the finite part is defined as the finite part of the ()A.35|) . The final result is: 



C/,i.(mi,m2,m3) = -Qf,^ 



A + C(mi,m2,m3)] , 



(A.37) 



where the normalized 3-point function is: 



C (mi, 1712,1713,] 



-- + logmg + Eo{mi, 1712, 1713) 



(A.38) 



with 



Eo{mi, 1712,1713) 



rrin 



TTlo 1 1 La . , 

log H 2 2-^("^l' "^2, ^3] 
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1 1 

— m2 



2 ' 2 2 ^ 
m2 m| — m2 



(A.39) 
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Some frequent special cases are: 

£„(M,™,m) = ^;^j-^-(^;^j-^) log-, (A.40) 

9 

Eo{m,m,M)=\og— + Eo{M,m,m). (A.41) 
The function (70(^1, m2, 777,3) is very similar to the previous one, with the result: 

Co{mi,m2, ms) = 5o(m2, ms) + mlCo{mi, m2, ms) = A + C{mi,m2, m^) + ^ . (A.42) 
And the last one is Si (mi, 7712): 



Appendix B 

Diagonalizing a squared mass matrix 



In this appendix we will show a group of formulas useful to diagonalize a 2 x 2 squared 
mass matrix M(2). We will fix our convention for the orthogonal matrixes that diagonalize 
M^'^\ showing the analytic expressions for the eigenvalues and the rotations angles, and 
finally how parametrize the sector in terms of the eigenvalues. 
For the mass Lagrangian in the electroweak basis 

Cm = ^0'^M(2)0/ (B.l) 

where 

are the electroweak eigenstates, and 



M(2) ^ (B.3) 

is the symmetric squared mass matrix that mix electroweak states. We define R as the 
orthogonal matrix that diagonalize M*^^); 

(cos 09 sin 009 \ ,^ , 

^ (B.4) 
— sin (j)ip cos j 

RM^^^R^ = diag{m2, ml} (B.5) 
The new mass eigenstates will be: 

= R(j)' (B.6) 
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with the mass eigenstates the Lagrangian can be written as 

Cm = \m\<Pl + \nl<Pl (B.7) 

The mass eigenvalues can be written in terms of the matrix elements as: 

1 



1,2 2 



h^^{a- 6)2 + 4c2 (B.8) 



where mKm"^) is the lightest (heaviest) particle and in general we can said that is a 
monotonic decreasing (increasing) function with the mixing term |c|. The mixing angle 
that diagonalize ()B.3|) can be extracted from the relation: 

tan(2(^) = (B.9) 
a — b 

to know the correct quadrant of ip is more useful this pair of relations: 
A/(a - 6)2 + 4c2 ^/(a - 6)^ + 4c2 



From ()B.8|) we can extract two relations: 

mi+m2 = a + 6 (B.H) 



-m? ^m\ = V(a - 6)^ + 4c2 (B.12) 
and from the invariant of the determinant: 

m^mj = a6 — c^. (B.13) 

These relations can be very useful to reparametrize the mass sector. We have put 
everything in terms of the three degrees of freedom of the matrix, i.e. a, 6, c. We can 
choose other sets to represent the degrees of freedom and give the rest in terms of these. 
Next we will present some examples of these reparametrizations. 

• Parameters mi, m2, 

a = ^ \m\{\ -cos2(y9) + m2(l + cos2v9)] (B.14) 

6 = 1 [mi(l + cos2^)+m2(l-cos2v?)] (B.15) 

c = — - \{rn\ — rr?^ sin2(y9] (B.16) 
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Parameters mi, a, c 



1 



+ m 



a — mr 



m2 = ^ + a 



2 



2c 

tan2<y5 



9 

a — mi ^ 

i a— mi 



Parameters mi,m2,c 



a = - 



b 



ml + ml± J {ml - mlY - 4c^ 



2 I 2 

+ m-^ =1= 
2c 



ml — m^y — 4c^ 



(B.17) 
(B.18) 
(B.19) 

(B.20) 
(B.21) 
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